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.

GGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAAT  >  NC_000913/3757802‑3758096
                                                                                                                                                  |                                                                                                                                                     
ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACaa                                                                                                                                                   >  2:595958/1‑150 (MQ=255)
ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACaa                                                                                                                                                   >  1:1580324/1‑150 (MQ=255)
ggAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACaa                                                                                                                                                   >  1:40738/1‑150 (MQ=255)
ggAATAAAGAACGTTGAAGCTGCCATAACACACTACTTAAAGTGTAAAAAAGTGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACaa                                                                                                                                                   <  1:1742484/150‑1 (MQ=255)
 gAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAAt                                                                                                                                                  >  1:82298/1‑150 (MQ=255)
 gAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATACATTAAACAAt                                                                                                                                                  >  2:155295/1‑150 (MQ=255)
  aaTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGAACATTGTATTCATTAAACAAtg                                                                                                                                                 >  1:1167157/1‑150 (MQ=255)
       aGAAAGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTAATGAGGATTTTCATTCGAAAATATGATCTGAATTAATGATTATCAAAATTAATCTAATAAAAAGCACATTTTATTCATTAAACAATTTGCtt                                                                                                                                            >  2:128536/1‑150 (MQ=255)
            gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttttt                                                                                                                                       >  1:1077179/1‑150 (MQ=255)
            gTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCttttttt                                                                                                                                       >  2:581816/1‑150 (MQ=255)
             ttGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTa                                                                                                                                      >  2:1209783/1‑150 (MQ=255)
                  gCTGCCATAACACACTCCTGAAAGTGAAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGATTTTTTTAAAtta                                                                                                                                 >  2:275608/1‑149 (MQ=255)
                   cTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCa                                                                                                                                <  1:581656/150‑1 (MQ=255)
                     gCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATa                                                                                                                              >  2:1465545/1‑150 (MQ=255)
                      ccATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATaa                                                                                                                             >  1:911229/1‑150 (MQ=255)
                          aaCACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAAt                                                                                                                         >  2:479108/1‑150 (MQ=255)
                          aaCACAATCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGTATTTACATTCGATATTATGATCAGTATTATTGATTATAAAAATTAATCTAATAAAAAGATAATTGTATAAATTATTAAATCTGATTTTTTTAAATTCATAAGatt                                                                                                                         >  1:299429/1‑148 (MQ=255)
                           acacacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTAATTAAACAATGTGCTTTTTTTAAATTCATAAGAAtt                                                                                                                        >  2:273718/1‑150 (MQ=255)
                              cacTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTg                                                                                                                     <  2:624309/150‑1 (MQ=255)
                               aatacttaaaGTGTAAAAAATGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGt                                                                                                                    <  2:509471/143‑1 (MQ=255)
                                cTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTa                                                                                                                   <  2:712946/150‑1 (MQ=255)
                                  ccTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTAtt                                                                                                                 >  1:546117/1‑150 (MQ=255)
                                   cTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTa                                                                                                                >  1:839822/1‑150 (MQ=255)
                                           gtAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAatttattt                                                                                                        >  2:1708982/1‑150 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       <  2:771099/150‑1 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       >  1:1298104/1‑150 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       >  1:1019391/1‑150 (MQ=255)
                                            tAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTATTAAATTCATAAGAATTTTGTATTAATTTAtttt                                                                                                       >  2:725442/1‑150 (MQ=255)
                                             aaaaaaGGTCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTa                                                                                                      <  1:1606439/150‑1 (MQ=255)
                                              aaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                                                                     <  1:378562/150‑1 (MQ=255)
                                              aaaaaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTaa                                                                                                     >  2:790175/1‑150 (MQ=255)
                                                 aaGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAAc                                                                                                  <  2:940749/150‑1 (MQ=255)
                                                   gggcgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATGTAATAAAAAGCACATTGTATTAATTAAACTATTTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACtt                                                                                                <  2:391590/150‑1 (MQ=255)
                                                    ggcgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACtta                                                                                               >  1:705947/1‑150 (MQ=255)
                                                     gcgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCAAAAGAATTTTTTATTAATTTATTTTAAAACttaa                                                                                              >  1:1385848/1‑150 (MQ=255)
                                                      cgggcATCTTACTGATGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaat                                                                                             <  1:791945/150‑1 (MQ=255)
                                                      cgggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCATTATTATTGATTATCAAAATTAATCAAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaat                                                                                             <  2:138878/150‑1 (MQ=255)
                                                        ggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaatta                                                                                           >  2:1251526/1‑150 (MQ=255)
                                                        ggcAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACttaatta                                                                                           <  1:1503161/150‑1 (MQ=255)
                                                           aGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTtaattaata                                                                                        >  2:1594394/1‑149 (MQ=255)
                                                            gCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCa                                                                                       >  2:1617024/1‑150 (MQ=255)
                                                                  tGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCt                                                                                 <  1:1242414/150‑1 (MQ=255)
                                                                       aTTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTAttt                                                                            <  1:1129233/150‑1 (MQ=255)
                                                                           tCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTcc                                                                        <  2:342387/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  2:209663/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  1:49640/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  2:777098/150‑1 (MQ=255)
                                                                            cATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCg                                                                       <  2:82298/150‑1 (MQ=255)
                                                                             aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                      <  1:626389/150‑1 (MQ=255)
                                                                             aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCgg                                                                      <  2:1500542/150‑1 (MQ=255)
                                                                             aTTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTcagt                                                                      >  1:488945/1‑147 (MQ=255)
                                                                              ttCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGa                                                                     >  2:779253/1‑150 (MQ=255)
                                                                               tCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGaa                                                                    >  2:528606/1‑150 (MQ=255)
                                                                               tCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGGGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGaa                                                                    <  2:1139340/150‑1 (MQ=255)
                                                                               acGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTACGGaa                                                                    >  1:760729/2‑150 (MQ=255)
                                                                                cGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCTGaaa                                                                   >  1:1540373/1‑150 (MQ=255)
                                                                                 gAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCTGAaat                                                                  >  2:1124952/1‑150 (MQ=255)
                                                                                 gAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaat                                                                  >  1:695051/1‑150 (MQ=255)
                                                                                   aaaTATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataat                                                               <  1:1475167/150‑1 (MQ=255)
                                                                                       atGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAaataataat                                                            <  1:75085/150‑1 (MQ=255)
                                                                                                attattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCgcgc                                                   <  1:622794/150‑1 (MQ=255)
                                                                                                 ttattGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCa                                                  <  2:1738939/150‑1 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCGATaaa                                             >  2:1732768/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaaa                                             >  2:1670832/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaaa                                             >  1:1618526/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATaaa                                             >  2:1748428/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGAGCAATaaa                                             >  1:903372/1‑150 (MQ=255)
                                                                                                      gATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTACTTATTTTTCAGGAAATAATAATGCTTCGAGCAATaaa                                             >  1:37243/1‑150 (MQ=255)
                                                                                                      aaTTATCAAAATTAATCTAATAAAAAGCACATTGTAATAATTAAACAATGTGCTTTATATAAATAAATAAGAATTTTGTATTAATTTAAATTAAAAATTAATTAATCATTTCCTTATTTTTCCGGAAAAAATAATGCGACGCGCAATaaa                                             >  1:903370/2‑150 (MQ=255)
                                                                                                         tATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGa                                          >  2:1140588/1‑150 (MQ=255)
                                                                                                          aTCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGAt                                         <  2:1580324/150‑1 (MQ=255)
                                                                                                            aaaaaTTAATCTAATAAAAAGCACTTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                       <  2:850020/149‑1 (MQ=255)
                                                                                                            aaaaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                       <  1:523702/149‑1 (MQ=255)
                                                                                                             aaaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATAAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAt                                      >  2:372857/1‑150 (MQ=255)
                                                                                                                 ttAATCTAATAAAAATCACATTTTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTcc                                  <  1:1366522/150‑1 (MQ=255)
                                                                                                                   aaTCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCac                                <  2:1298104/150‑1 (MQ=255)
                                                                                                                    aTCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacg                               >  1:656700/1‑150 (MQ=255)
                                                                                                                    aTCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacg                               <  1:1029805/150‑1 (MQ=255)
                                                                                                                      cTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                             >  1:617387/1‑150 (MQ=255)
                                                                                                                      ataataataaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgac                             <  1:1650369/149‑1 (MQ=255)
                                                                                                                        taaaaaaaGCACATTGTTTAATGAATACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCATTATTTTTCCGGAAATAATAATGAGTCGCGCAATAAATGATCATTTCAacgacga                           >  2:1057648/4‑150 (MQ=38)
                                                                                                                          taaAAAGCACATTTTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGTATTTTGTATTAATTTAATTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCacgacgacg                         <  1:316744/150‑1 (MQ=255)
                                                                                                                           aaaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                        >  2:1042248/1‑150 (MQ=255)
                                                                                                                                 cacaTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAg                  >  2:719440/1‑150 (MQ=255)
                                                                                                                                 cacaTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAg                  >  2:318329/1‑150 (MQ=255)
                                                                                                                                  acaTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGc                 >  2:406457/1‑150 (MQ=255)
                                                                                                                                         aTTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCg          >  2:207653/1‑150 (MQ=255)
                                                                                                                                         aTTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCg          <  1:581816/150‑1 (MQ=255)
                                                                                                                                          ttCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCgg         <  2:1392433/150‑1 (MQ=255)
                                                                                                                                            cATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTATTTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc       >  2:131378/1‑150 (MQ=255)
                                                                                                                                            cATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc       >  2:548015/1‑150 (MQ=255)
                                                                                                                                            cATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc       >  2:723822/1‑150 (MQ=255)
                                                                                                                                             aTTAAACATTGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa      >  2:221441/1‑150 (MQ=255)
                                                                                                                                             aTTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa      <  1:199567/150‑1 (MQ=255)
                                                                                                                                             aTTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa      <  1:1253277/150‑1 (MQ=255)
                                                                                                                                               tAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCaaa    <  2:705947/150‑1 (MQ=255)
                                                                                                                                                 aaCAATGTGCTTTTTTTAAATTCATAAGATTTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  <  1:1482221/150‑1 (MQ=255)
                                                                                                                                                 aaCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  >  2:1820159/1‑150 (MQ=255)
                                                                                                                                                 aaCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  >  1:519714/1‑150 (MQ=255)
                                                                                                                                                 aaCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat  >  1:741373/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAAT  >  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: