Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,684,016 A→G 8.2% S379G (AGC→GGC)  rstB → sensory histidine kinase of RstAB two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,0160AG 8.2% 328.7 / 16.0 123S379G (AGC→GGC) rstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base A (58/50);  new base G (7/3);  total (68/55)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 125 total aligned reads displayed.

ATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCG  >  NC_000913/1683867‑1684165
                                                                                                                                                     |                                                                                                                                                     
aTAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCa                                                                                                                                                       >  2:1857621/1‑150 (MQ=255)
 tAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAg                                                                                                                                                      >  1:1383487/1‑150 (MQ=255)
    cGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCg                                                                                                                                                   <  2:18220/150‑1 (MQ=255)
       ccTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGa                                                                                                                                                <  1:757316/150‑1 (MQ=255)
        cTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGTATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                                                                               <  2:468322/150‑1 (MQ=255)
        cTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                                                                               <  1:160477/150‑1 (MQ=255)
            ggTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCt                                                                                                                                           <  1:1023555/150‑1 (MQ=255)
              tACTGCCATTCCACGGTTGACACCAGCCTTCTACTGGCGGGGAATAGAGCGACATTAATTTTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggc                                                                                                                                         <  1:1241954/150‑5 (MQ=255)
              tACTGCCATTCAAAGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGATCGACATTAATTGTTGAGGATGATGGCCCAGGGATTCCCCCAGAAAACCGCGAACATCTCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggc                                                                                                                                         <  1:1807654/150‑5 (MQ=255)
                cTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGAAGGCCCAGGGATTGCCCCAGAACACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAc                                                                                                                                       >  1:273528/1‑150 (MQ=255)
                 tGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                                                                                      >  2:838520/1‑150 (MQ=255)
                 tGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                                                                                      >  1:1451354/1‑150 (MQ=255)
                 tGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGCATCGCTCAAcc                                                                                                                                      >  1:1613455/1‑150 (MQ=255)
                  gCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACcg                                                                                                                                     >  1:988605/1‑150 (MQ=255)
                     aTTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACcggcg                                                                                                                                  <  2:1608736/150‑1 (MQ=255)
                           cGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCCGCCGGGATCGCTCAACCGGCGGCTGCg                                                                                                                            >  1:1990077/1‑150 (MQ=255)
                               tGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggct                                                                                                                        >  2:425645/1‑150 (MQ=255)
                                gAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                                                                                       <  2:525882/150‑1 (MQ=255)
                                   aCCAGCCTGCTACTGTCGGGGAAAAGAGCGACATTAATTGTTGAGGAGGATGGCCCAGGGATTGCCCCAGAAACCCGCGAACATCTCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgggg                                                                                                                    <  1:1079771/150‑1 (MQ=255)
                                       gCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGggctggggctgg                                                                                                                <  2:489173/150‑1 (MQ=255)
                                        ccTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCCGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGc                                                                                                               >  1:1831411/1‑150 (MQ=255)
                                        ccTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGc                                                                                                               <  1:323672/150‑1 (MQ=255)
                                         cTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCa                                                                                                              >  2:1733771/1‑150 (MQ=255)
                                          tGCTACGGTCGGGGAATAGAGCTACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCaa                                                                                                             <  1:1807933/150‑1 (MQ=255)
                                           gCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAt                                                                                                            <  2:1454129/150‑1 (MQ=255)
                                             tACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTg                                                                                                          >  1:779174/1‑150 (MQ=255)
                                              aCTGTCGGGGAATAGAGCGACATTAAATGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACAGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCACACGGGAACGCTCAACCGGCGGCTGCGGGCTGGGGCAGGCAATTGt                                                                                                         >  1:847587/1‑150 (MQ=255)
                                               cTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCCGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTc                                                                                                        >  2:1849172/1‑150 (MQ=255)
                                                  tCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAc                                                                                                     <  1:1956206/150‑1 (MQ=255)
                                                  tCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAc                                                                                                     <  1:885151/150‑1 (MQ=255)
                                                  tCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAc                                                                                                     >  2:1494838/1‑150 (MQ=255)
                                                  tCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAc                                                                                                     <  1:425645/150‑1 (MQ=255)
                                                  tCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCGCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAc                                                                                                     <  2:1287789/150‑1 (MQ=255)
                                                        aaTAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtata                                                                                               <  2:244047/150‑1 (MQ=255)
                                                           agagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCa                                                                                            <  1:871547/150‑1 (MQ=255)
                                                           agagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTATGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCa                                                                                            >  2:1765924/1‑150 (MQ=255)
                                                            gagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGTTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCAc                                                                                           >  1:1632774/1‑150 (MQ=255)
                                                              gCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATTGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTg                                                                                         >  1:1590662/1‑150 (MQ=255)
                                                              ccGACATTAATTGGTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTg                                                                                         <  2:714951/149‑1 (MQ=255)
                                                                  cATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGATCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCaa                                                                                     <  2:225435/150‑1 (MQ=255)
                                                                  cATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCaa                                                                                     <  2:141047/150‑1 (MQ=255)
                                                                    ttAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATg                                                                                   <  2:1108488/150‑1 (MQ=255)
                                                                    ttAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATg                                                                                   >  1:1356609/1‑150 (MQ=255)
                                                                       aTTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGTCAATTGTCCACTCTATAGCACTGGCAATGGGc                                                                                <  2:1613455/150‑1 (MQ=255)
                                                                       aTTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGc                                                                                <  2:1717719/150‑1 (MQ=255)
                                                                       aTTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGc                                                                                <  1:685185/150‑1 (MQ=255)
                                                                        ttgttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCg                                                                               >  2:768046/1‑150 (MQ=255)
                                                                          gttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGt                                                                             >  2:1818096/1‑150 (MQ=255)
                                                                          gttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGt                                                                             >  1:1715406/1‑150 (MQ=255)
                                                                              aGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACgg                                                                         <  1:1573019/150‑1 (MQ=255)
                                                                               ggATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCACCCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt                                                                        >  2:1442111/1‑150 (MQ=255)
                                                                                  tgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTaa                                                                     >  1:645121/1‑150 (MQ=255)
                                                                                  tgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTaa                                                                     >  2:1584313/1‑150 (MQ=255)
                                                                                     tgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGGACCTTTTGTTCGCCTCGATCACAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCCGTCTATAGCACTGGCAATGGGCGGTACGGTTATTtg                                                                  >  2:979406/1‑150 (MQ=255)
                                                                                      ggCCCAGGGATTGCACCAGAAAACCGCGAACATATATTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCTGGCTGGGGCTGGAAATTGTCCACTCTATTGCACTGGCAATGGGCGGTACGGTTAATtgt                                                                 >  1:1255615/1‑150 (MQ=255)
                                                                                        cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGATCGCCTCGCTACCAGCCGGGATCGCTCACCCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCCCTCTATAGCTCTGGCAATGGGCGGGTCGGTTAATTGTGa                                                               >  1:108394/1‑150 (MQ=255)
                                                                                              gATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAACCGGGATCGCTCAACCGGCGGCTGCGGTCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTTTGACACCAg                                                         >  2:751194/1‑150 (MQ=255)
                                                                                                ttGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAGCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCg                                                       <  2:845372/150‑1 (MQ=255)
                                                                                                tggCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCg                                                       <  2:837973/148‑1 (MQ=255)
                                                                                                   ccccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAAc                                                    >  1:1874633/1‑150 (MQ=255)
                                                                                                     ccAGAAAACCGCGAACATATCTTTGCCCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTg                                                  >  1:803148/1‑150 (MQ=255)
                                                                                                      cAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTgg                                                 <  1:711826/150‑1 (MQ=255)
                                                                                                        gAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCAACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGggt                                               >  1:946376/1‑150 (MQ=255)
                                                                                                                cgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCg                                       <  2:1134225/150‑1 (MQ=255)
                                                                                                                 gAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGc                                      <  2:273528/150‑1 (MQ=255)
                                                                                                                    cATATCTTTGAACCTTTTGTGCGCCTCGATCCCAGCCGGGATCGCTCCACCGGCGGCTGCGGGCGGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtc                                   <  1:1816547/150‑1 (MQ=255)
                                                                                                                      tatCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtctc                                 >  2:1925217/1‑150 (MQ=255)
                                                                                                                        tCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGt                               >  2:1718291/1‑150 (MQ=255)
                                                                                                                         cTTTGACCCGGTAGTGCGCCTCGCCACCCGCCCGGATCGCTCAACCGGCGGCTGCGGCCGGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                              <  1:1246948/150‑1 (MQ=255)
                                                                                                                         cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                              >  2:1136028/1‑150 (MQ=255)
                                                                                                                         cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                              <  2:2018596/150‑1 (MQ=255)
                                                                                                                           ttGAACCTTTTGTTCGCCTCGATCCCCGCCGGGATCGCTCAACCGGCGGCGGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTa                            <  2:751887/150‑1 (MQ=255)
                                                                                                                           ttGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTa                            <  1:593011/150‑1 (MQ=255)
                                                                                                                             gAAACTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAATATATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAgg                          >  2:1028142/1‑149 (MQ=255)
                                                                                                                                 cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGc                      >  1:542666/1‑150 (MQ=255)
                                                                                                                                  ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGcc                     >  2:1875669/1‑150 (MQ=255)
                                                                                                                                  ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGcc                     >  1:673035/1‑150 (MQ=255)
                                                                                                                                     tGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACTCCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGtt                  >  1:32798/1‑150 (MQ=255)
                                                                                                                                     tGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGTCCGCTTCTCGTTTAGCTGGCCGtt                  >  1:863082/1‑150 (MQ=255)
                                                                                                                                     tGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGTCCGCTTCTCGTTTAGCTGGCCGtt                  >  1:1377370/1‑150 (MQ=255)
                                                                                                                                      gTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGGTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCATCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTa                 >  2:962279/1‑150 (MQ=255)
                                                                                                                                      gTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGAAATGGGCGGTACGGTTAATTGTGACACCAGCGATCTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTa                 >  2:1450035/1‑150 (MQ=255)
                                                                                                                                       ttCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt                >  1:11991/1‑150 (MQ=255)
                                                                                                                                       ttCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGGTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt                >  2:720655/1‑150 (MQ=255)
                                                                                                                                       ggttgagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt                >  2:431263/8‑150 (MQ=255)
                                                                                                                                          tgagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTATTTAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGCTTAGCTGGCCGTTATGGc             >  2:1051876/5‑150 (MQ=255)
                                                                                                                                          tgagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCAGGCAATTGTCCACTCTATAGCACTGGCAATCGGCGCTACGGTTAATTGTGACACAAGCGAACTGGGTGGTGACAGCTTCTCGTTTAGCTGGAAGTTATGGc             >  1:1038375/5‑150 (MQ=255)
                                                                                                                                          gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACGGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGc             <  1:267889/150‑1 (MQ=255)
                                                                                                                                           gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGGAATGGGCGGTACGGTTACTTGTGACACCAGCGAAGTGGGTGGTGCCCGCTTCTCGTTTAGCGGGCCGTTATGGCa            >  1:201295/4‑150 (MQ=255)
                                                                                                                                           gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCa            >  1:483831/4‑150 (MQ=255)
                                                                                                                                            agcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCAt           >  1:1445078/3‑150 (MQ=255)
                                                                                                                                            agcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCAt           >  1:400180/3‑150 (MQ=255)
                                                                                                                                              cGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTTTGACACCAGCGATCTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATaa         >  2:965393/1‑150 (MQ=255)
                                                                                                                                               gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGTGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc        >  2:747901/1‑150 (MQ=255)
                                                                                                                                               gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc        >  2:753080/1‑150 (MQ=255)
                                                                                                                                               gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc        >  1:1289766/1‑150 (MQ=255)
                                                                                                                                                 tCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACAt      <  2:704463/150‑1 (MQ=255)
                                                                                                                                                   ccAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATcc    >  2:947581/1‑150 (MQ=255)
                                                                                                                                                    cAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATccc   >  2:1869895/1‑150 (MQ=255)
                                                                                                                                                     aGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCg  <  2:1831411/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCG  >  NC_000913/1683867‑1684165

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