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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,147,3040GC 6.2% 401.3 / 15.8 129R56P (CGG→CCG) frwDputative enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (72/48);  new base C (2/6);  total (74/55)
Fisher's exact test for biased strand distribution p-value = 6.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

Only 100 of 135 total aligned reads displayed.

AGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCG  >  NC_000913/4147155‑4147452
                                                                                                                                                     |                                                                                                                                                    
aGTAACCGCCTGCGTACGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTGGTGCCTGCTAGAGACGTGGGGAGGCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATccgccg                                                                                                                                                      <  1:1044753/150‑1 (MQ=255)
   aaCCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGc                                                                                                                                                   <  2:1657057/150‑1 (MQ=255)
     ccGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCgg                                                                                                                                                 >  1:1792109/1‑150 (MQ=255)
          tGCGTAAGCGGCGTGGCGCATACTTATATGGCGGGGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCGGGGAACGGAGAAGCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatc                                                                                                                                            <  2:1358592/150‑4 (MQ=255)
           gCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAATGTTGGGCCTGTGAGAGAATTGGGGAGTCAGCATTTCAACTCAGGGCGAGCTGGGAACGGAGAATGGTTTATCGGTTGAGTATATCAGCCGGGCGGGTGTTg                                                                                                                                           >  2:917296/1‑150 (MQ=255)
           gCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcc                                                                                                                                           <  2:1743176/150‑5 (MQ=255)
           gCGTAAGCGGCGGGGCGCATACTTATATGGCGTCGGAACGGCTGGATAAGTTGTGCCTGTTAGAGAAGCGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcc                                                                                                                                           <  2:58561/150‑5 (MQ=255)
            cGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGc                                                                                                                                          <  1:99395/150‑1 (MQ=255)
            cGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcct                                                                                                                                          <  1:363657/150‑6 (MQ=255)
            cGTAAGCGGCGGGGCGCATACTTATATGGCGGCGGAACGGCGGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcct                                                                                                                                          <  1:1760229/150‑6 (MQ=255)
               aaGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGctct                                                                                                                                       >  1:1017646/1‑150 (MQ=255)
                aGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTg                                                                                                                                      <  2:1256769/150‑1 (MQ=255)
                aGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTg                                                                                                                                      >  2:701894/1‑150 (MQ=255)
                 gcggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGt                                                                                                                                     >  1:1043714/1‑150 (MQ=255)
                   ggcgTGGCGCATACTTATATGGCGGTGGAACGGCTGGAAAAGTTGTGCGTGTTAGAGAAGTGGGTAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTg                                                                                                                                   <  2:1865110/150‑1 (MQ=255)
                     cgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGAt                                                                                                                                 <  1:1316925/150‑1 (MQ=255)
                          cgcATACTTATATGGCGGCGGAACGGCTGGAAAAGGTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACgg                                                                                                                            <  1:1136340/150‑1 (MQ=255)
                            cATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGat                                                                                                                          <  2:1691121/150‑1 (MQ=255)
                               aCTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATc                                                                                                                       <  2:847131/150‑1 (MQ=255)
                                cTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCg                                                                                                                      <  2:1123530/150‑1 (MQ=255)
                                   ttatGGCGGCGGAACGGCTGGAACAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGc                                                                                                                   <  2:277138/149‑1 (MQ=255)
                                       ggcggcGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCGGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGc                                                                                                               <  2:247858/150‑1 (MQ=255)
                                             ggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGc                                                                                                         <  2:1565815/150‑1 (MQ=255)
                                                tcgtcttgAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTGAGGTCGCGCTGGGAACGGAGAATCGTTTTGCGGATGAGGTTATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATAGAGCTTGCCGGTGCCGa                                                                                                      <  2:859151/143‑1 (MQ=255)
                                                aCGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGa                                                                                                      <  1:902718/150‑1 (MQ=255)
                                                 cGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAg                                                                                                     <  1:1846945/150‑1 (MQ=255)
                                                  ggCTGGTACAGTGGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGc                                                                                                    <  1:759694/150‑1 (MQ=255)
                                                  ggCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGc                                                                                                    <  2:20388/150‑1 (MQ=255)
                                                   gCTGGAAAAGTGGTGCCTGTTAGAGAAGTGGGGAGTCAGTATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCg                                                                                                   <  1:569787/150‑1 (MQ=255)
                                                    cTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGa                                                                                                  >  1:402298/1‑150 (MQ=255)
                                                    cTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGa                                                                                                  >  1:1078202/1‑150 (MQ=255)
                                                     tGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGATGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAt                                                                                                 >  2:1268298/1‑150 (MQ=255)
                                                      ggAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAtt                                                                                                >  2:673618/1‑150 (MQ=255)
                                                      ggAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAtt                                                                                                >  2:1373732/1‑150 (MQ=255)
                                                        aaaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTg                                                                                              <  2:1402910/150‑1 (MQ=255)
                                                         aaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATTTCGAGCTTGCCGGTGCCGAGCGATTTGa                                                                                             >  2:891540/1‑150 (MQ=255)
                                                         aaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGa                                                                                             <  1:666484/150‑1 (MQ=255)
                                                          atgTTTTGCCTGTTTGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCTGGCGGATGTTGCTCTGTTGATTACGTATATCGAGCTTGCCGTTGCCGAGCGATTTGaa                                                                                            <  2:1195646/148‑1 (MQ=255)
                                                          aaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGaa                                                                                            <  1:516212/150‑1 (MQ=255)
                                                          aaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGaa                                                                                            <  1:893384/150‑1 (MQ=255)
                                                          aaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGaa                                                                                            <  1:1806473/150‑1 (MQ=255)
                                                            gTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACa                                                                                          >  2:1757928/1‑150 (MQ=255)
                                                            gTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACa                                                                                          >  2:342290/1‑150 (MQ=255)
                                                               gtgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGTTGTTGCTCTGGTGATTACGGATATCGAGCTTGCCGGTGCCGCGCGATTGGAACATTg                                                                                       >  2:538059/1‑150 (MQ=255)
                                                               gtgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTg                                                                                       >  2:1377780/1‑150 (MQ=255)
                                                               gtgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTATGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTg                                                                                       >  2:169977/1‑150 (MQ=255)
                                                                tgCCGGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGATTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGc                                                                                      <  1:142752/150‑1 (MQ=255)
                                                                 gCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGcc                                                                                     >  2:1766991/1‑150 (MQ=255)
                                                                 gCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGcc                                                                                     >  1:1503524/1‑150 (MQ=255)
                                                                   cTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGc                                                                                   <  2:967010/150‑1 (MQ=255)
                                                                    tGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTTCGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCt                                                                                  >  2:1619117/1‑150 (MQ=255)
                                                                              gTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATg                                                                        >  1:445198/1‑150 (MQ=255)
                                                                               tGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgc                                                                       >  2:887707/1‑150 (MQ=255)
                                                                               tGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgc                                                                       >  2:879674/1‑150 (MQ=255)
                                                                                ggggAGTCAGCATTGATTCTGAGGGCGCGCAGGTAACGGAGTATCGTATAGCGGAGGAGGATATGTGCAGGGCGGATGTTGCTCTGTTGATTACGGAAATCGTGCTTGCCGTTGACGAGCGATATGAACATTGCCTCTATGTGCAATgca                                                                      <  2:1601333/150‑1 (MQ=255)
                                                                                ggggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGTGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                      <  1:640231/150‑1 (MQ=255)
                                                                                  ggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgcagc                                                                    <  2:1329373/150‑1 (MQ=255)
                                                                                  ggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgcagc                                                                    >  2:485899/1‑150 (MQ=255)
                                                                                    aGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCAt                                                                  >  1:575816/1‑150 (MQ=255)
                                                                                        aGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTAc                                                              <  2:1830024/150‑1 (MQ=255)
                                                                                         gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACg                                                             <  2:704068/150‑1 (MQ=255)
                                                                                         gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACg                                                             >  1:20757/1‑150 (MQ=255)
                                                                                         ccaatGAAACGCAGGGCGCGCTCGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGACGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACg                                                             <  1:49083/146‑1 (MQ=255)
                                                                                          cATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCCGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGAGCAATGCAGCAACTACGc                                                            >  2:462107/1‑150 (MQ=255)
                                                                                              gAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTc                                                        >  2:1312999/1‑150 (MQ=255)
                                                                                                   tCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCg                                                   <  1:1654033/150‑1 (MQ=255)
                                                                                                     aGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg                                                 <  1:1081177/150‑1 (MQ=255)
                                                                                                      gggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGa                                                >  1:21603/1‑150 (MQ=255)
                                                                                                        gcgcgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGc                                              >  2:1704710/1‑150 (MQ=255)
                                                                                                         cgcgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGcc                                             <  1:253311/150‑1 (MQ=255)
                                                                                                           cgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGc                                           >  2:838175/1‑150 (MQ=255)
                                                                                                            gcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCa                                          >  1:119429/1‑150 (MQ=255)
                                                                                                             cTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAg                                         >  2:1905337/1‑150 (MQ=255)
                                                                                                                  aaCGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGt                                    >  2:1257666/1‑150 (MQ=255)
                                                                                                                   aCGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTa                                   >  1:1243468/1‑150 (MQ=255)
                                                                                                                        gaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAg                              >  2:740029/1‑150 (MQ=255)
                                                                                                                        gaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAg                              >  2:631222/1‑150 (MQ=255)
                                                                                                                             gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                         >  1:1728800/1‑150 (MQ=255)
                                                                                                                             gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCgg                         >  1:589736/1‑150 (MQ=255)
                                                                                                                               ttAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTg                       >  2:340883/1‑150 (MQ=255)
                                                                                                                                tAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATCTCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgc                      <  2:1040597/150‑1 (MQ=255)
                                                                                                                                tAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGCGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgc                      >  2:1495930/1‑150 (MQ=255)
                                                                                                                                tAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgc                      >  1:1304589/1‑150 (MQ=255)
                                                                                                                                tAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgc                      <  2:668055/150‑1 (MQ=255)
                                                                                                                                  gCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTAGACCATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCCCGGTgcgc                    >  2:1213060/1‑150 (MQ=255)
                                                                                                                                   cGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCa                   >  1:6225/1‑150 (MQ=255)
                                                                                                                                    ggATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCaa                  >  1:787491/1‑150 (MQ=255)
                                                                                                                                       cagcaacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGAGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGt               >  1:259762/8‑150 (MQ=255)
                                                                                                                                         aGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGc             >  2:425365/1‑150 (MQ=255)
                                                                                                                                          ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCt            >  2:575986/1‑150 (MQ=255)
                                                                                                                                          ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCt            >  1:1291456/1‑150 (MQ=255)
                                                                                                                                           aacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCtt           >  1:1202809/4‑150 (MQ=255)
                                                                                                                                            atatCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCttt          >  2:491697/1‑150 (MQ=255)
                                                                                                                                             tatCCGCCGGGCGGATGTTGCTCTGTTGATTACGGCTCTCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTc         <  2:1017646/150‑1 (MQ=255)
                                                                                                                                             tatCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTc         >  1:948279/1‑150 (MQ=255)
                                                                                                                                             tatCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTc         >  2:1852290/1‑150 (MQ=255)
                                                                                                                                               tCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTg       >  2:419191/1‑150 (MQ=255)
                                                                                                                                               tCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTACTg       >  1:1462733/1‑150 (MQ=255)
                                                                                                                                                ccgccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTgc      >  1:339122/1‑150 (MQ=255)
                                                                                                                                                    cgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGcgccg  <  1:1721630/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCG  >  NC_000913/4147155‑4147452

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: