Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,572,928 C→G 6.6% A759P (GCA→CCA)  pqqL ← putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,9280CG 6.6% 311.8 / 10.4 103A759P (GCA→CCA) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base C (41/54);  new base G (5/2);  total (47/56)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 105 total aligned reads displayed.

CTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAC  >  NC_000913/1572789‑1573077
                                                                                                                                           |                                                                                                                                                     
cTGTATGCTCCAGATGCCTGTTCACGACTATTAACACGTAGATCTTTTGCCAGTGCGACGGTAAAAGCAGCGAGCTCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCa                                                                                                                                             <  2:1197617/150‑1 (MQ=255)
  gTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                                                                           <  1:1112754/150‑5 (MQ=255)
   tATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcaca                                                                                                                                          <  1:1564483/150‑6 (MQ=255)
   tATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGGCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAttt                                                                                                                                          <  2:1434609/150‑1 (MQ=255)
   tATGCTCCAGAAGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAACCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAttt                                                                                                                                          <  1:159268/150‑1 (MQ=255)
    aTGCTCCAGATGCCTGTTCACGATTATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                         <  1:159443/150‑1 (MQ=255)
    aTGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                         <  2:1164472/150‑1 (MQ=255)
    aTGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCAGCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                                                                         <  1:1031899/150‑1 (MQ=255)
      gCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGt                                                                                                                                       <  1:834500/150‑1 (MQ=255)
          cAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtt                                                                                                                                   >  2:1713538/1‑150 (MQ=255)
           aGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCttt                                                                                                                                  <  2:234206/150‑1 (MQ=255)
              tGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                                                                               >  1:771346/1‑150 (MQ=255)
              tGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                                                                               <  1:391159/150‑1 (MQ=255)
              tGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGGCGGCAGATTAACAGGTGGCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                                                                               <  1:951160/150‑1 (MQ=255)
                    ttCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAAc                                                                                                                         >  1:575468/1‑150 (MQ=255)
                       aCGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACTCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTCTGCCACAGGTTCATTTTGTTCTTTTACAGTAACAGa                                                                                                                      >  1:267569/1‑150 (MQ=255)
                       aCGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                                                                      >  1:1573235/1‑150 (MQ=255)
                        cGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGAt                                                                                                                     >  2:1503151/1‑150 (MQ=255)
                          aaTATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGc                                                                                                                   <  1:984754/150‑1 (MQ=255)
                              ttAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTg                                                                                                               <  2:664599/150‑1 (MQ=255)
                                 acacGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTcc                                                                                                            >  2:923212/1‑150 (MQ=255)
                                  cacGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCg                                                                                                           >  2:1064312/1‑150 (MQ=255)
                                     gTAGATCTTTTTCCAGTGCGACGTTAAAAGCATCTAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcg                                                                                                        <  2:1045935/150‑1 (MQ=255)
                                           gttttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGTTTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg                                                                                                  <  2:907264/149‑1 (MQ=255)
                                           cTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg                                                                                                  <  2:1733884/150‑1 (MQ=255)
                                            ttttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGt                                                                                                 >  1:1247066/1‑150 (MQ=255)
                                              ttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTa                                                                                               <  2:722940/150‑1 (MQ=255)
                                               tGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                              >  2:439627/1‑150 (MQ=255)
                                               tGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAGTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa                                                                                              >  2:1311697/1‑150 (MQ=255)
                                                 ccAGTGCGACGTTAAAAGCATCGTGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTCACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                                                                            <  2:124372/150‑1 (MQ=255)
                                                 ccAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                                                                            >  2:721335/1‑150 (MQ=255)
                                                   aGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGt                                                                                          <  1:1713538/150‑1 (MQ=255)
                                                     tGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGttt                                                                                        <  2:72325/150‑1 (MQ=255)
                                                        gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                                                                     <  1:382591/150‑1 (MQ=255)
                                                        gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                                                                     <  1:1664293/150‑1 (MQ=255)
                                                           gTTAAAAGCATCGAGCGCCATACGCGTCGGCAGAGTAACAGGTGTCCGGGATTCATAACGCTTCCACTGTGAAACCTGTGCCCCAGGTTCATGTTGTTCTTTTACAGTTACCGATGCGTTGTCCGTCGCGCGAGTGTATGGTTTTCCTGc                                                                                  >  2:1135637/1‑150 (MQ=255)
                                                            ttAAAAGCATCGAGCGCCATACGCGTCGGCTGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTTTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGAAAATGGTTTACCTGCg                                                                                 <  1:152902/150‑1 (MQ=255)
                                                            ttAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACCGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCTCAGGCTCATTTTGTTCTTTTCCAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCg                                                                                 <  1:828095/150‑1 (MQ=255)
                                                            ttAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCg                                                                                 <  1:649718/150‑1 (MQ=255)
                                                            ttAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCg                                                                                 <  1:1268033/150‑1 (MQ=255)
                                                              aaaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGc                                                                               >  2:1308580/1‑150 (MQ=255)
                                                              aaaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGc                                                                               <  2:863883/150‑1 (MQ=255)
                                                               aaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCt                                                                              <  1:145262/150‑1 (MQ=255)
                                                                 aGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTaa                                                                            >  1:786072/1‑150 (MQ=255)
                                                                  gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAAt                                                                           >  2:1240011/1‑150 (MQ=255)
                                                                   cATCGAGCGCCATAGGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATg                                                                          <  2:168178/150‑1 (MQ=255)
                                                                   cATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATg                                                                          >  2:216330/1‑150 (MQ=255)
                                                                    aTCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATgg                                                                         <  1:1064312/150‑1 (MQ=255)
                                                                     tCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGc                                                                        <  2:484154/150‑1 (MQ=255)
                                                                      cGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCg                                                                       >  2:628093/1‑150 (MQ=255)
                                                                         gcgcCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTAGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGaaa                                                                    >  1:1298921/1‑149 (MQ=255)
                                                                            ccATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag                                                                 <  2:1447869/150‑1 (MQ=255)
                                                                                aCGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgtg                                                             <  2:1660409/150‑1 (MQ=255)
                                                                                aCGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgtg                                                             <  1:874865/150‑1 (MQ=255)
                                                                                  gcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGtt                                                           <  2:467661/150‑1 (MQ=255)
                                                                                   cgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGttt                                                          <  2:1060798/150‑1 (MQ=255)
                                                                                      cGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgat                                                       <  1:1503151/150‑1 (MQ=255)
                                                                                       ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGTGGCTAATGGCGAATCAGAGTGTttgatt                                                      <  2:505324/150‑1 (MQ=255)
                                                                                       ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgatt                                                      >  1:1105713/1‑150 (MQ=255)
                                                                                       ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgatt                                                      <  2:1305565/150‑1 (MQ=255)
                                                                                       ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGGTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgatt                                                      >  2:1056181/1‑150 (MQ=255)
                                                                                        gCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCAATGTGAAACCTGTGACACAGGTTCATTTTGTTCATTTACAGAAACCGATCAGTTGTCCGTCCCGCGAGGTTATGGTTTACCTGCGGCTAATGGCGAATCATCGTGTttgattg                                                     >  1:839859/1‑150 (MQ=255)
                                                                                         cAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGGCGCGCGAGTTAATGGTTTCCCTGCGGCTAAGGGCGAATCCGAGTGTttgattga                                                    >  2:246053/1‑150 (MQ=255)
                                                                                           taTTAACAGGTGTCCGTGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATc                                                  <  2:267569/149‑1 (MQ=255)
                                                                                           gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATc                                                  <  1:1536617/150‑1 (MQ=255)
                                                                                              tcactggTGTCCGGGAATCATAACGCTTCCACTGTGAAACGTGTGCCAGAGGTTCATTTTGTTCTTTTACAGTATCCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTa                                               <  1:254075/145‑1 (MQ=255)
                                                                                               aaCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCAGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTaa                                              >  2:754094/1‑150 (MQ=255)
                                                                                                 cAGGTGTCCGGGATTCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGt                                            <  1:1732689/150‑1 (MQ=255)
                                                                                                 cAGGTGTCCGGGAATCATAACTCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACTGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTAACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGt                                            <  1:258539/150‑1 (MQ=255)
                                                                                                      gtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAcgc                                       >  1:600235/1‑150 (MQ=255)
                                                                                                       tCCGGGAATCATAACGCATCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAcgcg                                      <  2:1287071/150‑1 (MQ=255)
                                                                                                         cGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGAGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTa                                    >  2:1088642/1‑150 (MQ=255)
                                                                                                          gggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTaa                                   >  1:187714/1‑150 (MQ=255)
                                                                                                           ggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTGGATTGATACTAAGTAACGCGTAAt                                  >  2:512442/1‑150 (MQ=255)
                                                                                                                cATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGAGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACg                             >  1:331425/1‑150 (MQ=255)
                                                                                                                 aTAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGc                            <  1:216330/150‑1 (MQ=255)
                                                                                                                 aTAACGCGTCCACTTGGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTTCCTGCGGCTAATGGCGATTCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGc                            <  2:889414/150‑1 (MQ=255)
                                                                                                                     cGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACg                        <  1:1437235/150‑1 (MQ=255)
                                                                                                                           cACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTg                  >  2:1358902/1‑150 (MQ=255)
                                                                                                                           cACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTg                  <  1:628093/150‑1 (MQ=255)
                                                                                                                            caaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGt                 >  1:1266560/8‑150 (MQ=255)
                                                                                                                            aCTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAcgcg                                      >  2:333614/1‑129 (MQ=255)
                                                                                                                            aCTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAcgcg                                      <  1:333614/129‑1 (MQ=255)
                                                                                                                             cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGTGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtg                >  1:1402100/1‑149 (MQ=255)
                                                                                                                             cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtg                >  2:164553/1‑149 (MQ=255)
                                                                                                                             cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtc                >  2:771266/1‑150 (MQ=255)
                                                                                                                             aaaatgaaCCTGTGGCACAGGTTCATTTTGTTATTTTACAGTAACCGATGCGTTGTCCGTCTCGCGAGTTAATGGTTTACCTGCGTCTAATGGCGAATCAGATTGTTTGATTGATCCTAAGTAACGCGTAAATAACGCCACGAGTTTGta                >  2:1657273/7‑149 (MQ=255)
                                                                                                                              tgtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtct               <  1:371514/150‑1 (MQ=255)
                                                                                                                              aaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtct               >  1:1000609/6‑150 (MQ=255)
                                                                                                                              aaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGATAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACTAGTTTGtct               >  2:1369267/6‑150 (MQ=255)
                                                                                                                               aatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtctt              >  2:686002/5‑150 (MQ=255)
                                                                                                                                tgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtcttc             >  2:1343262/1‑150 (MQ=255)
                                                                                                                                tgAAACCTGTGCCACAGGTTCATTTTGTACTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTTATGGTTTACCTGCGGCTAATCGCGAATCAGAGTGATTGATAGAGAATAATTAACGAGTAATTAAAGCAACGAGTTTGtcttc             >  2:1056900/1‑150 (MQ=255)
                                                                                                                                    aCCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCg         <  1:1435767/150‑1 (MQ=255)
                                                                                                                                    aCCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCg         <  2:786072/150‑1 (MQ=255)
                                                                                                                                          gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTa   >  1:528830/1‑150 (MQ=255)
                                                                                                                                          gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTa   >  1:1669819/1‑150 (MQ=255)
                                                                                                                                          gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCAAATCAGAGTGTTTGTTTGTTCCTAAGTGACGCGTAATTAACGCCACGAGATTGTCTTCTGCGAGATTa   >  1:282454/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTTCATTTTGTTCTTTTACGGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGAGCCTAAGTACCGCCTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAc  >  2:115669/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAc  >  2:1365848/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
CTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAC  >  NC_000913/1572789‑1573077

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: