Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,147,304 G→C 5.9% R56P (CGG→CCG)  frwD → putative PTS enzyme IIB component FrwD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,147,3040GC 5.9% 458.3 / 13.1 136R56P (CGG→CCG) frwDputative PTS enzyme IIB component FrwD
Reads supporting (aligned to +/- strand):  ref base G (59/69);  new base C (2/6);  total (61/75)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 144 total aligned reads displayed.

AGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGC  >  NC_000913/4147155‑4147453
                                                                                                                                                     |                                                                                                                                                     
aGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATccgccg                                                                                                                                                       <  2:557150/150‑1 (MQ=255)
aGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATccgccg                                                                                                                                                       >  1:1721004/1‑150 (MQ=255)
  tAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCggg                                                                                                                                                     <  1:251764/150‑1 (MQ=255)
    aCCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCg                                                                                                                                                   <  2:1773174/150‑1 (MQ=255)
      cGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                                                                 <  1:1924312/150‑1 (MQ=255)
          tGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGtt                                                                                                                                             >  1:661030/1‑150 (MQ=255)
          tGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatc                                                                                                                                             <  2:272652/150‑4 (MQ=255)
           gCGTAAGGGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCTGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcc                                                                                                                                            <  2:2102690/150‑5 (MQ=255)
           gCGTAAGCTTCGTGGCTCATACTTATATGGCGTCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCATCATTGAAACTCATGTCGAGCTGGGAACTGAGAATCGTTTAGCGGATGAGGATATCCGCCCGTCGGAtatac                                                                                                                                            <  2:1676853/150‑5 (MQ=255)
           gCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcc                                                                                                                                            <  2:1977120/150‑5 (MQ=255)
           gCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGGCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGACATCCGCCCGGCGGAtatcc                                                                                                                                            <  1:1433691/150‑5 (MQ=255)
            cGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGc                                                                                                                                           >  2:1987834/1‑150 (MQ=255)
            cGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGc                                                                                                                                           <  1:234848/150‑1 (MQ=255)
            cGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGct                                                                                                                                          <  1:1185062/150‑1 (MQ=255)
              tAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCTGATGTTGctc                                                                                                                                         <  2:1765170/150‑1 (MQ=255)
              tAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGctc                                                                                                                                         <  1:280059/150‑1 (MQ=255)
              tAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGctc                                                                                                                                         >  2:1406215/1‑150 (MQ=255)
                aGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTg                                                                                                                                       <  1:2251745/150‑1 (MQ=255)
                  cggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                                                                     <  1:646430/150‑1 (MQ=255)
                   ggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTg                                                                                                                                    <  1:753694/150‑1 (MQ=255)
                    gcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGa                                                                                                                                   <  1:1539335/150‑1 (MQ=255)
                         gcgcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACg                                                                                                                              <  2:342205/150‑1 (MQ=255)
                          cgcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACgg                                                                                                                             >  1:450734/1‑150 (MQ=255)
                          cgcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACgg                                                                                                                             >  1:1823043/1‑150 (MQ=255)
                           gcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCTTTGAAACTCAGGGGGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCTCCGGGCGGATGTTGCTCTGTTGATTACGGa                                                                                                                            <  1:510757/150‑1 (MQ=255)
                           gcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGa                                                                                                                            <  1:776654/150‑1 (MQ=255)
                              tACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGatat                                                                                                                         >  2:329088/1‑150 (MQ=255)
                                cTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGAGGTTGCTCTGTTGATTACGGATATCg                                                                                                                       >  1:1714130/1‑150 (MQ=255)
                                 ttATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGa                                                                                                                      >  1:1854432/1‑150 (MQ=255)
                                      tGGCGGCGGAACGCCTGGAAAAGTTGTGACTCTTAGAGAAGTGGGGAGTCAGCATTTAAAATCAGGGCGCTCTCCGAACGGAGAATCGTTTAGCGGATGAGGATATACGCCGGGCGGATCTTCATCTCATCATAAAGGATATAGAGCTTg                                                                                                                 >  2:1530845/1‑150 (MQ=255)
                                        gcggcgGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGcc                                                                                                               <  1:1747236/150‑1 (MQ=255)
                                         cggcggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCg                                                                                                              >  1:1479361/1‑150 (MQ=255)
                                          ggcggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTTCTCTGTTGATTACGGATATCGAGCTTGCCgg                                                                                                             <  2:157960/150‑1 (MQ=255)
                                            cggAACGGCTGCAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACCGAAAATCGTTTAGCGGATGAGGATATCAGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTg                                                                                                           >  2:1064646/1‑150 (MQ=255)
                                              gAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGcc                                                                                                         <  1:2395691/150‑1 (MQ=255)
                                                 cGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAg                                                                                                      <  1:1369962/150‑1 (MQ=255)
                                                  ggCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGAGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGc                                                                                                     <  2:1540558/150‑1 (MQ=255)
                                                   gCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCg                                                                                                    <  2:878202/150‑1 (MQ=255)
                                                       gAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAttt                                                                                                <  1:1941828/150‑1 (MQ=255)
                                                           aGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAAc                                                                                            <  2:165754/150‑1 (MQ=255)
                                                             ttGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACAt                                                                                          <  1:329088/150‑1 (MQ=255)
                                                              tgtgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACAtt                                                                                         <  2:1353942/150‑1 (MQ=255)
                                                               gtgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTg                                                                                        <  1:933200/150‑1 (MQ=255)
                                                                 gCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGcc                                                                                      <  2:1494182/150‑1 (MQ=255)
                                                                   cTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGc                                                                                    <  2:2416842/150‑1 (MQ=255)
                                                                    tGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCt                                                                                   <  1:139148/150‑1 (MQ=255)
                                                                     gTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTa                                                                                  <  2:302151/150‑1 (MQ=255)
                                                                      ttAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAt                                                                                 <  2:1908379/150‑1 (MQ=255)
                                                                       tAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtg                                                                                <  1:713993/150‑1 (MQ=255)
                                                                        agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgt                                                                               <  1:92308/150‑1 (MQ=255)
                                                                        agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgt                                                                               >  2:1248302/1‑150 (MQ=255)
                                                                        agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCGCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgt                                                                               >  2:1571138/1‑150 (MQ=255)
                                                                         gagaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgtg                                                                              <  1:1987834/150‑1 (MQ=255)
                                                                           gaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCa                                                                            >  2:2481092/1‑150 (MQ=255)
                                                                                ggggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTTCTCTGTTGATTAGGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                       <  2:1721004/150‑1 (MQ=255)
                                                                                ggggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                       >  2:405616/1‑150 (MQ=255)
                                                                                ggggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                       <  1:2290368/150‑1 (MQ=255)
                                                                                 gggAGGCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgcag                                                                      <  1:1582315/150‑1 (MQ=255)
                                                                                   gAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgcagca                                                                    <  1:793133/150‑1 (MQ=255)
                                                                                     gTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATc                                                                  >  2:2165290/1‑150 (MQ=255)
                                                                                     gTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATc                                                                  >  1:1107166/1‑150 (MQ=255)
                                                                                       cAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAGTTGAACATTGCCGCTATGTGCAATGCAGCATCTa                                                                >  2:1897143/1‑150 (MQ=255)
                                                                                         gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACg                                                              <  2:1522589/150‑1 (MQ=255)
                                                                                         gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACg                                                              <  2:1714130/150‑1 (MQ=255)
                                                                                           aTTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCa                                                            >  1:1090802/1‑150 (MQ=255)
                                                                                              gAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTc                                                         >  2:2423476/1‑150 (MQ=255)
                                                                                              gAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTc                                                         >  1:2121918/1‑150 (MQ=255)
                                                                                               aaaCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTcc                                                        <  2:1854432/150‑1 (MQ=255)
                                                                                                 aCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTg                                                      <  2:2314450/150‑1 (MQ=255)
                                                                                                 aCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTg                                                      >  2:1144703/1‑150 (MQ=255)
                                                                                                 aCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTg                                                      <  2:794348/150‑1 (MQ=255)
                                                                                                  cTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGc                                                     >  2:1841203/1‑150 (MQ=255)
                                                                                                      gggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGa                                                 >  1:2210855/1‑150 (MQ=255)
                                                                                                       ggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAg                                                >  2:1902719/1‑150 (MQ=255)
                                                                                                       ggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAg                                                <  1:1790343/150‑1 (MQ=255)
                                                                                                        gcgcgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGc                                               <  2:656783/150‑1 (MQ=255)
                                                                                                            gcTGGGAACGGAGAATCGTATAGAGGATGAGGATATCCACCGGTCGGATGTTCCTCTGTTGATTACGGATATAGAGCTTGCAGGTCCCGAGCGATATGAACATTGACGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAgaagca                                           >  2:2388857/1‑145 (MQ=255)
                                                                                                                 gAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCggg                                      <  1:1522871/150‑1 (MQ=255)
                                                                                                                      gagaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATg                                 <  1:1696912/150‑1 (MQ=255)
                                                                                                                             gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                          <  1:2102861/150‑1 (MQ=255)
                                                                                                                             gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                          <  1:964308/150‑1 (MQ=255)
                                                                                                                               ttAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGCTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTg                        <  1:2454966/150‑1 (MQ=255)
                                                                                                                                 aGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgcg                      >  2:1676020/1‑150 (MQ=255)
                                                                                                                                 aGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgcg                      >  1:2487900/1‑150 (MQ=255)
                                                                                                                                   cGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGTCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCa                    >  2:290666/1‑150 (MQ=255)
                                                                                                                                   cGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCa                    >  1:335944/1‑150 (MQ=255)
                                                                                                                                    ggATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCaa                   >  2:1199895/1‑150 (MQ=255)
                                                                                                                                       tGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGt                >  1:403669/1‑150 (MQ=255)
                                                                                                                                       gagcaacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGt                >  2:646377/8‑150 (MQ=255)
                                                                                                                                         aGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGc              >  1:874002/1‑150 (MQ=255)
                                                                                                                                          ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCAt                                                           <  1:1100943/104‑1 (MQ=255)
                                                                                                                                          ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCAt                                                           >  2:1100943/1‑104 (MQ=255)
                                                                                                                                          ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCt             >  1:1156944/1‑150 (MQ=255)
                                                                                                                                          ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCt             >  2:54565/1‑150 (MQ=255)
                                                                                                                                           gATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCtt            <  1:91261/150‑1 (MQ=255)
                                                                                                                                              atCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCt         >  1:1989004/1‑150 (MQ=255)
                                                                                                                                               tCCGCCGTGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTTAACATTGCCGCTATGTTCAATGCATCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTg        <  2:2287102/150‑1 (MQ=255)
                                                                                                                                                 cgccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTgcg      >  2:1506448/1‑150 (MQ=255)
                                                                                                                                                   ccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGcgcc    >  1:2224579/1‑150 (MQ=255)
                                                                                                                                                   ccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGcgcc    >  1:362792/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGC  >  NC_000913/4147155‑4147453

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