Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 578,333 A→C K76T (AAA→ACA) 
K3Q (AAG→CAG) 
rzpD →
rzoD →
DLP12 prophage; putative prophage endopeptidase RzpD
DLP12 prophage; putative prophage lysis lipoprotein RzoD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,3330AC99.2% 314.3 / ‑3.6 119coding (227/462 nt)
coding (7/183 nt)
rzpD
rzoD
DLP12 prophage; putative prophage endopeptidase RzpD
DLP12 prophage; putative prophage lysis lipoprotein RzoD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (59/59);  minor base . (1/0);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

Only 100 of 121 total aligned reads displayed.

ATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGA  >  NC_000913/578197‑578468
                                                                                                                                        |                                                                                                                                       
aTCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGc                                                                                                                                      >  1:1171973/1‑140 (MQ=255)
 tCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCt                                                                                                                                     >  1:246631/1‑140 (MQ=255)
 tCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCt                                                                                                                                     <  1:1666200/140‑1 (MQ=255)
 tCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCt                                                                                                                                     >  2:200600/1‑140 (MQ=255)
   gCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGa                                                                                                                                   <  2:795588/140‑1 (MQ=255)
   gCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGa                                                                                                                                   <  1:1407880/140‑1 (MQ=255)
     cTACAAAGTCCAGCGCGACGAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGaaa                                                                                                                                 >  2:1889318/1‑140 (MQ=255)
      tACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGaaaa                                                                                                                                >  1:1080513/1‑140 (MQ=255)
       aCAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAat                                                                                                                               >  2:1694266/1‑140 (MQ=255)
          aaGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAatgat                                                                                                                            <  1:200600/140‑1 (MQ=255)
          aaGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAatgat                                                                                                                            <  2:1090275/140‑1 (MQ=255)
           aGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAatgatg                                                                                                                           >  1:1626395/1‑140 (MQ=255)
             tCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGct                                                                                                                         <  2:405081/140‑1 (MQ=37)
             tCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGct                                                                                                                         <  2:1597271/140‑1 (MQ=37)
             tCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGct                                                                                                                         >  1:2118389/1‑140 (MQ=37)
             tCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGct                                                                                                                         <  2:190898/140‑1 (MQ=37)
              ccAGCGCGACAAAAATGCCAGAGAACTCAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGAGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGctct                                                                                                                       <  2:524275/140‑1 (MQ=255)
               cAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGa                                                                                                                                   >  2:1412718/1‑128 (MQ=25)
               cAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGa                                                                                                                                   <  1:1412718/128‑1 (MQ=25)
               cAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGctct                                                                                                                       >  1:930037/1‑140 (MQ=25)
                aGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTg                                                                                                                      <  2:1929415/140‑1 (MQ=21)
                  cgcgACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCg                                                                                                                    >  1:87364/1‑140 (MQ=25)
                       cAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGtgatg                                                                                                               >  1:1730301/1‑140 (MQ=25)
                       cAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGtgatg                                                                                                               >  1:2211545/1‑140 (MQ=25)
                       cAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGtgatg                                                                                                               >  2:1966188/1‑140 (MQ=25)
                                cAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCg                                                                                                      <  1:211659/140‑1 (MQ=21)
                                cAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCg                                                                                                      <  2:660974/140‑1 (MQ=21)
                                 agagAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGc                                                                                                     <  1:1694266/140‑1 (MQ=21)
                                    gaACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTgg                                                                                                  >  1:542396/1‑140 (MQ=37)
                                       cTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcg                                                                                               <  1:2096231/140‑1 (MQ=21)
                                          aaGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcg                                                                                            >  1:1402816/1‑140 (MQ=37)
                                           aGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcgt                                                                                           <  2:1928378/140‑1 (MQ=21)
                                             cTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGtcgtcgtcg                                                                                         >  2:622758/1‑140 (MQ=25)
                                             cTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGtcgtcgtcg                                                                                         >  2:1074784/1‑140 (MQ=21)
                                               aGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGt                                                                                       >  2:53374/1‑140 (MQ=18)
                                               aGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGt                                                                                       >  1:983219/1‑140 (MQ=18)
                                               aGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGt                                                                                       >  2:1305498/1‑140 (MQ=18)
                                               aGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGt                                                                                       >  1:412947/1‑140 (MQ=18)
                                                 cGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg                                                                                     <  2:443921/140‑1 (MQ=17)
                                                 cGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg                                                                                     >  2:1598163/1‑140 (MQ=18)
                                                 cGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg                                                                                     >  1:2034733/1‑140 (MQ=17)
                                                 cGAACGCGGCAATTACTGACAGGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg                                                                                     <  2:394168/140‑1 (MQ=17)
                                                   aaCGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGca                                                                                   <  1:323725/140‑1 (MQ=17)
                                                   aaCGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGca                                                                                   <  2:2060771/140‑1 (MQ=17)
                                                    aCGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGcac                                                                                  >  1:1858694/1‑140 (MQ=18)
                                                     cgcgGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGcaca                                                                                 >  2:1458724/1‑140 (MQ=18)
                                                          cAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCaaa                                                                            >  2:1920356/1‑140 (MQ=18)
                                                           aaTTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAg                                                                           <  1:1966188/140‑1 (MQ=17)
                                                           aaTTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAg                                                                           <  2:845935/140‑1 (MQ=17)
                                                             ttACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCa                                                                         <  2:2211545/140‑1 (MQ=17)
                                                             ttACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCa                                                                         <  2:740506/140‑1 (MQ=17)
                                                               aCTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAg                                                                           <  1:1523151/136‑1 (MQ=17)
                                                               aCTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAg                                                                           >  2:1523151/1‑136 (MQ=18)
                                                               aCTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGt                                                                       <  1:450438/140‑1 (MQ=17)
                                                               aCTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGt                                                                       <  1:177620/140‑1 (MQ=17)
                                                               aCTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGt                                                                       <  1:143223/140‑1 (MQ=17)
                                                                cTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTc                                                                      <  2:412947/140‑1 (MQ=17)
                                                                                               tgcGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCTCATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCt                                       <  2:198544/140‑1 (MQ=17)
                                                                                               tgcGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCt                                       <  2:1585510/140‑1 (MQ=17)
                                                                                               tgcGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAACTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCt                                       <  1:2158608/140‑1 (MQ=17)
                                                                                                 cgcTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTcc                                     >  2:618446/1‑140 (MQ=17)
                                                                                                   cTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCgg                                   >  1:198751/1‑140 (MQ=17)
                                                                                                   cTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCgg                                   >  1:2018837/1‑140 (MQ=17)
                                                                                                   cTCGATGCAAAATACACGAAGGAGTTAGCTGATGCG‑‑ACCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCg                                 >  2:46642/1‑140 (MQ=255)
                                                                                                    tagATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGc                                  <  1:2149306/138‑1 (MQ=17)
                                                                                                    tCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGc                                  <  2:1887252/140‑1 (MQ=18)
                                                                                                    tCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGc                                  <  1:1074784/140‑1 (MQ=18)
                                                                                                    tCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGc                                  <  2:1626395/140‑1 (MQ=17)
                                                                                                    tCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGc                                  <  2:1483984/140‑1 (MQ=18)
                                                                                                    tCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGc                                  <  1:1601482/140‑1 (MQ=18)
                                                                                                    tCGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGc                                  <  1:622758/140‑1 (MQ=18)
                                                                                                     cGATGCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCg                                 >  1:2040181/1‑140 (MQ=17)
                                                                                                         gCAAAATACACGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGa                             <  2:2118389/140‑1 (MQ=18)
                                                                                                                  acGAAGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGc                    >  1:2127938/1‑140 (MQ=17)
                                                                                                                     gaGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTc                 >  2:2074995/2‑140 (MQ=17)
                                                                                                                     aaGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTc                 >  2:541581/1‑140 (MQ=17)
                                                                                                                      aGGAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcc                <  1:1086266/140‑1 (MQ=17)
                                                                                                                       ggAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccc               >  2:873352/1‑140 (MQ=17)
                                                                                                                       ggAGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccc               >  1:2154763/1‑140 (MQ=17)
                                                                                                                         aGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccccc             >  2:2189222/1‑140 (MQ=17)
                                                                                                                         aGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccccc             >  2:38496/1‑140 (MQ=17)
                                                                                                                         aGTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccccc             >  2:1476516/1‑140 (MQ=17)
                                                                                                                          gTTGGCTGTTGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCGCCGGCGTGGATAATGCAGCCTcccccc            <  2:131365/140‑1 (MQ=255)
                                                                                                                          gTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccc            <  2:2034733/140‑1 (MQ=17)
                                                                                                                          gTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccc            <  1:1476516/140‑1 (MQ=17)
                                                                                                                          gTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccc            <  1:1888012/140‑1 (MQ=17)
                                                                                                                          gTTAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccc            <  1:2155419/140‑1 (MQ=17)
                                                                                                                           ttAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATg                        >  2:2217616/1‑127 (MQ=17)
                                                                                                                           ttAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATg                        <  1:2217616/127‑1 (MQ=17)
                                                                                                                            tAGCTGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGa          <  2:1352776/140‑1 (MQ=17)
                                                                                                                                tGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTgg      >  2:2120421/1‑140 (MQ=17)
                                                                                                                                tGATGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTgg      >  2:743715/1‑140 (MQ=17)
                                                                                                                                  aTGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCa    >  2:2186915/1‑140 (MQ=17)
                                                                                                                                  aTGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCa    >  1:665193/1‑140 (MQ=17)
                                                                                                                                  aTGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCa    >  1:923455/1‑140 (MQ=17)
                                                                                                                                  aTGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCa    >  2:1403125/1‑140 (MQ=17)
                                                                                                                                   tGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAg   >  2:788488/1‑140 (MQ=17)
                                                                                                                                   tGCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAg   >  1:127463/1‑140 (MQ=17)
                                                                                                                                    gCGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGa  <  2:542396/140‑1 (MQ=17)
                                                                                                                                     cGACAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGa  <  1:1825165/139‑1 (MQ=17)
                                                                                                                                        |                                                                                                                                       
ATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGA  >  NC_000913/578197‑578468

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