Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 585,863 T→C 6.1% intergenic (‑230/‑194) ompT ← / → pauD DLP12 prophage; outer membrane protease VII; outer membrane protein 3b/tRNA‑OTHER

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913585,8630TC 6.1% 345.0 / 9.5 114intergenic (‑230/‑194)ompT/pauDDLP12 prophage; outer membrane protease VII; outer membrane protein 3b/tRNA‑OTHER
Reads supporting (aligned to +/- strand):  ref base T (54/53);  new base C (5/2);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 116 total aligned reads displayed.

CACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCA  >  NC_000913/585714‑586012
                                                                                                                                                     |                                                                                                                                                     
cACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCtt                                                                                                                                                       >  2:583002/1‑150 (MQ=255)
   tGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTag                                                                                                                                                    <  1:1344226/150‑1 (MQ=255)
    gTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGATTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTaga                                                                                                                                                   >  2:602974/1‑150 (MQ=255)
      ttATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGGGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAgg                                                                                                                                                 <  2:525156/150‑1 (MQ=255)
       tataTTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATCTTTCGGTTTTCCTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGc                                                                                                                                                <  2:935008/150‑1 (MQ=255)
       tataTTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGc                                                                                                                                                <  1:566252/150‑1 (MQ=255)
          atTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCata                                                                                                                                             >  2:1428656/1‑150 (MQ=255)
          atTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCata                                                                                                                                             >  2:1145153/1‑150 (MQ=255)
             ttGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTCTAAAGGCatggca                                                                                                                                          <  1:1039424/150‑5 (MQ=255)
              tGTTTAGTATAGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTATTATGCTTTTGTGGCGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAg                                                                                                                                         <  2:1471351/150‑1 (MQ=255)
                      aTTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa                                                                                                                                 >  2:366392/1‑150 (MQ=255)
                      aTTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATa                                                                                                                                 <  1:1644570/150‑1 (MQ=255)
                         gactGATTATCACAATGCGCTATCTGTTTTTGGTTTAAATATCTGTTAATGTTTCATATTTCGGGTGTACTGTGTGGGTTTTGTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCa                                                                                                                              <  1:110261/147‑1 (MQ=255)
                         gAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCAAATAAGCGAAAATAGCa                                                                                                                              >  2:763981/1‑150 (MQ=255)
                         gAAAGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTCTCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCa                                                                                                                              <  2:1080884/150‑1 (MQ=255)
                             gATTATCACAATGCGCTATCTGTTTATGGTTTAATTATCTGTTATTGTTTCATATTACGGTATTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCATTTAGAGGCATATAAGCGAAAATAGCATGAg                                                                                                                          >  1:890897/1‑150 (MQ=255)
                              aTTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTGGAGGCATATAAGCGAAAATAGCATGAgg                                                                                                                         <  1:1912486/150‑1 (MQ=255)
                                 aTCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGttt                                                                                                                      >  2:1110391/1‑150 (MQ=255)
                                 aTCACAATGCGCTATCTGTTTTTGCTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGAATATAAGCTAAAATAGAATGAGGttt                                                                                                                      >  1:287039/1‑150 (MQ=255)
                                  tCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTa                                                                                                                     >  2:253796/1‑150 (MQ=255)
                                  tCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTa                                                                                                                     <  2:1027075/150‑1 (MQ=255)
                                   cacaATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTAt                                                                                                                    >  2:253055/1‑150 (MQ=255)
                                     caATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATcc                                                                                                                  <  2:277795/150‑1 (MQ=255)
                                       aTGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTc                                                                                                                <  1:1308509/150‑1 (MQ=255)
                                        tGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCa                                                                                                               >  1:658407/1‑150 (MQ=255)
                                            cTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCATCAATTa                                                                                                           >  1:977052/1‑150 (MQ=255)
                                                cTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATTAGGTTTATCCTCAATTACTAt                                                                                                       >  2:42930/1‑150 (MQ=255)
                                                    ttttGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCTTGAGGTTTATCCTCAATTACTATGttt                                                                                                   >  1:1753043/1‑150 (MQ=255)
                                                       tGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGtttttt                                                                                                <  2:857274/150‑1 (MQ=255)
                                                        ggTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTa                                                                                               <  1:925623/150‑1 (MQ=255)
                                                        ggTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGAATATAAGCGAAAATAGAATGAGGTTTATCCTCAATTACTATGTTTTTTa                                                                                               >  1:1516561/1‑150 (MQ=255)
                                                         gTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAg                                                                                              >  1:1781955/1‑150 (MQ=255)
                                                         gTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAg                                                                                              <  2:647007/150‑1 (MQ=255)
                                                         gTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAg                                                                                              <  1:433302/150‑1 (MQ=255)
                                                          tttAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGt                                                                                             <  2:1359293/150‑1 (MQ=255)
                                                          tttAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGt                                                                                             >  2:888395/1‑150 (MQ=255)
                                                          tttAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATCTATTTTAGt                                                                                             >  2:880388/1‑150 (MQ=255)
                                                           ttAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCATCAATTACTATGTTTTTTAGTa                                                                                            >  1:1179196/1‑150 (MQ=255)
                                                              aTTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaa                                                                                         <  1:735122/150‑1 (MQ=255)
                                                               ttATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaa                                                                                        <  1:52349/150‑1 (MQ=255)
                                                                 aTCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaaa                                                                                      <  1:1890731/150‑1 (MQ=255)
                                                                   cTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACaaaaaat                                                                                    >  2:1754314/1‑149 (MQ=255)
                                                                    tGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAaga                                                                                   >  2:238770/1‑150 (MQ=255)
                                                                     tttATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAagag                                                                                  <  1:1741361/149‑1 (MQ=255)
                                                                         ttGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGAc                                                                              <  2:1665854/150‑1 (MQ=255)
                                                                              tCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGACTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTATTACAAAAAAGAGGGACAAAAc                                                                         <  2:1800724/150‑1 (MQ=255)
                                                                               cATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACt                                                                        >  2:1494845/1‑150 (MQ=255)
                                                                                atatTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTg                                                                       <  2:1564377/150‑1 (MQ=255)
                                                                                  atTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgag                                                                     >  2:228827/1‑150 (MQ=255)
                                                                                  atTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgag                                                                     >  1:1312347/1‑150 (MQ=255)
                                                                                   tttCTGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTgaga                                                                    <  2:1118165/150‑1 (MQ=255)
                                                                                   tttCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACATAAAAGAGGGACAAAACTgagc                                                                    >  2:1861913/1‑149 (MQ=255)
                                                                                       ggTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGAcaca                                                                <  2:387430/150‑1 (MQ=255)
                                                                                       ggTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGAcaca                                                                >  2:1288002/1‑150 (MQ=255)
                                                                                       ggTTTTACTGTGTGGTTTTTTTATACTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGAcaca                                                                >  2:175884/1‑150 (MQ=255)
                                                                                                 tgttgtTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCg                                                      <  1:1262649/150‑1 (MQ=255)
                                                                                                  gtgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGc                                                     >  2:607098/1‑150 (MQ=255)
                                                                                                   tgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACATATAAGGCCTCGCa                                                    >  1:538206/1‑150 (MQ=255)
                                                                                                   tgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCATCGCa                                                    >  1:1436547/1‑150 (MQ=255)
                                                                                                   cgatctttttATGCTTTTGTGTTGCTTTTATCTATTTAAGTGCCATTCGTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATTTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCa                                                    <  2:1548664/145‑1 (MQ=255)
                                                                                                    ggTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                   >  2:1839039/1‑150 (MQ=255)
                                                                                                    ggTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                   >  1:1322678/1‑150 (MQ=255)
                                                                                                    ggTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaa                                                   >  1:1843624/1‑150 (MQ=255)
                                                                                                     gTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAAt                                                  >  1:1075965/1‑150 (MQ=255)
                                                                                                      tttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATg                                                 <  1:473815/150‑1 (MQ=255)
                                                                                                       ttttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATgg                                                <  2:1593009/150‑1 (MQ=255)
                                                                                                        tttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGc                                               >  1:1758157/1‑150 (MQ=255)
                                                                                                         ttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCaagatcg                                              >  1:1548664/1‑145 (MQ=255)
                                                                                                         ttttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCt                                              >  1:697254/1‑150 (MQ=255)
                                                                                                           ttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTg                                            <  1:763377/150‑1 (MQ=255)
                                                                                                              tGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCaa                                         <  1:1858133/150‑1 (MQ=255)
                                                                                                               gCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAg                                        <  1:1114849/150‑1 (MQ=255)
                                                                                                               gCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAg                                        <  1:1486222/150‑1 (MQ=255)
                                                                                                               gCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATACTCTATTACTATGTTTTTTTGTACAAAAAAGAGGGACAAATCTGAGACACATAAGGCGTCGCAATGGCTTgaatg                                        >  1:1195820/1‑146 (MQ=255)
                                                                                                                cTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAgg                                       >  1:1521859/1‑150 (MQ=255)
                                                                                                                cTTTAGTGGGGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTGTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAACGCCTCGCAATGGCTTGCAAgg                                       <  2:213279/150‑1 (MQ=255)
                                                                                                                   ttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCtt                                    <  1:440469/150‑1 (MQ=255)
                                                                                                                   ttGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCtt                                    <  1:44258/150‑1 (MQ=255)
                                                                                                                     gtggtgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATACAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTa                                  <  1:5912/150‑1 (MQ=255)
                                                                                                                       ggtgCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACa                                >  1:509135/1‑150 (MQ=255)
                                                                                                                          gCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGtt                            >  1:1185949/1‑150 (MQ=255)
                                                                                                                               tATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAATAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGATTGCAAGGATTTACATGTTTTGa                        >  1:255450/1‑150 (MQ=255)
                                                                                                                               tATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGttttta                        >  1:1353767/1‑148 (MQ=255)
                                                                                                                               tATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                        <  2:658407/150‑1 (MQ=255)
                                                                                                                               tATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                        <  1:697326/150‑1 (MQ=255)
                                                                                                                               tATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGATTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGa                        <  1:280538/150‑1 (MQ=255)
                                                                                                                                aTCTAATTAAGTGCCATGCCTTTAGAGGCATATAACCGAAAATAGCATGAGGTGTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAg                       <  2:289828/150‑1 (MQ=255)
                                                                                                                                 tCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAgg                      >  2:1831163/1‑150 (MQ=255)
                                                                                                                                   tATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTa                    >  2:47966/1‑150 (MQ=255)
                                                                                                                                       tAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTgg                >  1:533311/1‑150 (MQ=255)
                                                                                                                                        cgcttataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTggg               >  2:929308/8‑150 (MQ=255)
                                                                                                                                         gcttataTGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGa              >  2:1800653/7‑150 (MQ=255)
                                                                                                                                         aGTGCCATTCCTTTAGATGCATATAAGCGAAAATTGCATGAGTTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGa              <  2:1268358/150‑1 (MQ=255)
                                                                                                                                          gTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGAc             >  1:1728030/1‑150 (MQ=255)
                                                                                                                                           tGCGATGCATTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCGTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACg            <  2:890897/150‑1 (MQ=255)
                                                                                                                                            gCCATGCATTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgt           >  2:1408291/1‑150 (MQ=255)
                                                                                                                                             ccATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgtg          >  1:885792/1‑150 (MQ=255)
                                                                                                                                               attccTTTAGAGGCATATAAGCGAAAATAGCTTGAGGTTTATCCTCAATTACTATGTTTTTTAGTAGAAACTAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGggtg        <  1:306726/147‑1 (MQ=255)
                                                                                                                                                tGCCTCTAAAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGa       >  1:318774/1‑150 (MQ=255)
                                                                                                                                                     ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCa  <  1:1460173/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACTGTTTATATTTTGTTTAGTATTGAATGATTATCACAATGCGCTATCTGTTTTTGGTTTAATTATCTGTTATTGTTTCATATTTCGGTTTTACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCA  >  NC_000913/585714‑586012

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