Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,757,947:1 +T 5.1% intergenic (‑120/+70) yiaY ← / ← selB L‑threonine dehydrogenase/selenocysteinyl‑tRNA‑specific translation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,757,9471.T 5.1% 550.4 / 24.0 136intergenic (‑120/+70)yiaY/selBL‑threonine dehydrogenase/selenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base . (76/53);  new base T (5/2);  total (81/55)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 141 total aligned reads displayed.

GGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAAT  >  NC_000913/3757802‑3758096
                                                                                                                                                  |                                                                                                                                                     
ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACaa                                                                                                                                                   >  2:595958/1‑150 (MQ=255)
ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACaa                                                                                                                                                   >  1:1580324/1‑150 (MQ=255)
ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACaa                                                                                                                                                   >  1:40738/1‑150 (MQ=255)
ggAATAAAGAACGTTGAAGCTGCCATAACACACTACTTAAAGTGTAAAAAAGTGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACaa                                                                                                                                                   <  1:1742484/150‑1 (MQ=255)
 gAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAAt                                                                                                                                                  >  1:82298/1‑150 (MQ=255)
 gAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATACATTAA‑ACAAt                                                                                                                                                  >  2:155295/1‑150 (MQ=255)
  aaTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGAACATTGTATTCATTAA‑ACAAtg                                                                                                                                                 >  1:1167157/1‑150 (MQ=255)
       aGAAAGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTAATGAGGATTTTCATTCGAAAATATGATCTGAATTAATGATTATCAAAATTAATCTAATAAAAAGCACATTTTATTCATTAA‑ACAATTTGCtt                                                                                                                                            >  2:128536/1‑150 (MQ=255)
            gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCttttttt                                                                                                                                       >  1:1077179/1‑150 (MQ=255)
            gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCttttttt                                                                                                                                       >  2:581816/1‑150 (MQ=255)
             ttGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTa                                                                                                                                      >  2:1209783/1‑150 (MQ=255)
                  gCTGCCATAACACACTCCTGAAAGTGAAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGATTTTTTTAAAtta                                                                                                                                 >  2:275608/1‑149 (MQ=255)
                   cTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCa                                                                                                                                <  1:581656/150‑1 (MQ=255)
                     gCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATa                                                                                                                              >  2:1465545/1‑150 (MQ=255)
                      ccATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATaa                                                                                                                             >  1:911229/1‑150 (MQ=255)
                          aaCACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAAt                                                                                                                         >  2:479108/1‑150 (MQ=255)
                          aaCACAATCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGTATTTACATTCGATATTATGATCAGTATTATTGATTATAAAAATTAATCTAATAAAAAGATAATTGTATAAATTAT‑TAAATCTGATTTTTTTAAATTCATAAGatt                                                                                                                         >  1:299429/1‑148 (MQ=255)
                           acacacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTAATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAAtt                                                                                                                        >  2:273718/1‑150 (MQ=255)
                              cacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTg                                                                                                                     <  2:624309/150‑1 (MQ=255)
                               aatacttaaaGTGTAAAAAATGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGt                                                                                                                    <  2:509471/143‑1 (MQ=255)
                                cTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTa                                                                                                                   <  2:712946/150‑1 (MQ=255)
                                  ccTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTAtt                                                                                                                 >  1:546117/1‑150 (MQ=255)
                                   cTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTa                                                                                                                >  1:839822/1‑150 (MQ=255)
                                           gtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatttattt                                                                                                        >  2:1708982/1‑150 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       <  2:771099/150‑1 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       >  1:1298104/1‑150 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       >  1:1019391/1‑150 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTATTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       >  2:725442/1‑150 (MQ=255)
                                             aaaaaaGGTCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTa                                                                                                      <  1:1606439/150‑1 (MQ=255)
                                              aaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                                                                     <  1:378562/150‑1 (MQ=255)
                                              aaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                                                                     >  2:790175/1‑150 (MQ=255)
                                                 aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                                                                  <  2:940749/150‑1 (MQ=255)
                                                   gggcgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATGTAATAAAAAGCACATTGTATTAATTAA‑ACTATTTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACtt                                                                                                <  2:391590/150‑1 (MQ=255)
                                                    ggcgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACtta                                                                                               >  1:705947/1‑150 (MQ=255)
                                                     gcgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCAAAAGAATTTTTTATTAATTTATTTTAAAACttaa                                                                                              >  1:1385848/1‑150 (MQ=255)
                                                      cgggcATCTTACTGATGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaat                                                                                             <  1:791945/150‑1 (MQ=255)
                                                      cgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCATTATTATTGATTATCAAAATTAATCAAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaat                                                                                             <  2:138878/150‑1 (MQ=255)
                                                        ggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaatta                                                                                           >  2:1251526/1‑150 (MQ=255)
                                                        ggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaatta                                                                                           <  1:1503161/150‑1 (MQ=255)
                                                           aGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTtaattaata                                                                                        >  2:1594394/1‑149 (MQ=255)
                                                            gCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCa                                                                                       >  2:1617024/1‑150 (MQ=255)
                                                                  tGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCt                                                                                 <  1:1242414/150‑1 (MQ=255)
                                                                       aTTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAttt                                                                            <  1:1129233/150‑1 (MQ=255)
                                                                           tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTcc                                                                        <  2:342387/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  2:209663/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  1:49640/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  2:777098/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  2:82298/150‑1 (MQ=255)
                                                                             aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                      <  1:626389/150‑1 (MQ=255)
                                                                             aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                      <  2:1500542/150‑1 (MQ=255)
                                                                             aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTcagt                                                                      >  1:488945/1‑147 (MQ=255)
                                                                              ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGa                                                                     >  2:779253/1‑150 (MQ=255)
                                                                               tCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGaa                                                                    >  2:528606/1‑150 (MQ=255)
                                                                               tCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGGGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGaa                                                                    <  2:1139340/150‑1 (MQ=255)
                                                                               acGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTACGGaa                                                                    >  1:760729/2‑150 (MQ=255)
                                                                                cGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCTGaaa                                                                   >  1:1540373/1‑150 (MQ=255)
                                                                                 gAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCTGAaat                                                                  >  2:1124952/1‑150 (MQ=255)
                                                                                 gAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaat                                                                  >  1:695051/1‑150 (MQ=255)
                                                                                   aaaTATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGC‑TTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataat                                                               <  1:1475167/150‑1 (MQ=255)
                                                                                       atGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataataat                                                            <  1:75085/150‑1 (MQ=255)
                                                                                                attattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCgcgc                                                   <  1:622794/150‑1 (MQ=255)
                                                                                                 ttattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCa                                                  <  2:1738939/150‑1 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCGATaaa                                             >  2:1732768/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaaa                                             >  2:1670832/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaaa                                             >  1:1618526/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaaa                                             >  2:1748428/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGAGCAATaaa                                             >  1:903372/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTACTTATTTTTCAGGAAATAATAATGCTTCGAGCAATaaa                                             >  1:37243/1‑150 (MQ=255)
                                                                                                      aaTTATCAAAATTAATCTAATAAAAAGCACATTGTAATAATTAA‑ACAATGTGCTTTATATAAATAAATAAGAATTTTGTATTAATTTAAATTAAAAATTAATTAATCATTTCCTTATTTTTCCGGAAAAAATAATGCGACGCGCAATaaa                                             >  1:903370/2‑150 (MQ=255)
                                                                                                         tATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGa                                          >  2:1140588/1‑150 (MQ=255)
                                                                                                          aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                         <  2:1580324/150‑1 (MQ=255)
                                                                                                            aaaaaTTAATCTAATAAAAAGCACTTTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                       <  2:850020/149‑1 (MQ=255)
                                                                                                            aaaaaTTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                       <  1:523702/149‑1 (MQ=255)
                                                                                                             aaaaTTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATAAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAt                                      >  2:372857/1‑150 (MQ=255)
                                                                                                                 ttAATCTAATAAAAATCACATTTTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTcc                                  <  1:1366522/150‑1 (MQ=255)
                                                                                                                   aaTCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCac                                <  2:1298104/150‑1 (MQ=255)
                                                                                                                    aTCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacg                               >  1:656700/1‑150 (MQ=255)
                                                                                                                    aTCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacg                               <  1:1029805/150‑1 (MQ=255)
                                                                                                                      cTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                             >  1:617387/1‑150 (MQ=255)
                                                                                                                      ataataataaGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                             <  1:1650369/149‑1 (MQ=255)
                                                                                                                        taaaaaaaGCACATTGT‑TTAATGAATACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCGGAAATAATAATGAGTCGCGCAATAAATGATCATTTCAacgacga                           >  2:1057648/4‑150 (MQ=38)
                                                                                                                          taaAAAGCACATTTTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGTATTTTGTATTAATTTAATTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacgacg                         <  1:316744/150‑1 (MQ=255)
                                                                                                                           aaaaaGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                        >  2:1042248/1‑150 (MQ=255)
                                                                                                                                 cacaTTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAg                  >  2:719440/1‑150 (MQ=255)
                                                                                                                                 cacaTTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAg                  >  2:318329/1‑150 (MQ=255)
                                                                                                                                  acaTTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGc                 >  2:406457/1‑150 (MQ=255)
                                                                                                                                         aTTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCg          >  2:207653/1‑150 (MQ=255)
                                                                                                                                         aTTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCg          <  1:581816/150‑1 (MQ=255)
                                                                                                                                          ttCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCgg         <  2:1392433/150‑1 (MQ=255)
                                                                                                                                            cATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTATTTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc       >  2:131378/1‑150 (MQ=255)
                                                                                                                                            cATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc       >  2:548015/1‑150 (MQ=255)
                                                                                                                                            cATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc       >  2:723822/1‑150 (MQ=255)
                                                                                                                                             aTTAA‑ACATTGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa      >  2:221441/1‑150 (MQ=255)
                                                                                                                                             aTTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa      <  1:199567/150‑1 (MQ=255)
                                                                                                                                             aTTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa      <  1:1253277/150‑1 (MQ=255)
                                                                                                                                               tAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCaaa    <  2:705947/150‑1 (MQ=255)
                                                                                                                                                 a‑aCAATGTGCTTTTTTTAAATTCATAAGATTTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  <  1:1482221/150‑1 (MQ=255)
                                                                                                                                                 a‑aCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  >  2:1820159/1‑150 (MQ=255)
                                                                                                                                                 a‑aCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  >  1:519714/1‑150 (MQ=255)
                                                                                                                                                 a‑aCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  >  1:741373/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAA‑ACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAAT  >  NC_000913/3757802‑3758096

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑