Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,660,125 A→C 5.5% D373A (GAC→GCC)  metK → methionine adenosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,660,1250AC 5.5% 250.3 / 5.1 91D373A (GAC→GCC) metKmethionine adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (46/39);  new base C (2/3);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GACGACAAGATCATCCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCA  >  NC_002947/5659977‑5660274
                                                                                                                                                    |                                                                                                                                                     
gacgacAAGATCATCCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGcgac                                                                                                                                                      <  2:1827588/150‑1 (MQ=255)
    acAAGATCATCCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGaca                                                                                                                                                  <  1:1411011/150‑1 (MQ=255)
      aaGATCATCCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcacc                                                                                                                                                >  1:634640/1‑150 (MQ=255)
        gATCATCCAGCTGGTGCGCGAGTGCTTCGACCTGCGCCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcacctt                                                                                                                                              <  1:1350299/150‑1 (MQ=255)
          tcatcCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGTCGCCGacagtcttc                                                                                                                                            <  2:637284/150‑7 (MQ=255)
          tcatcCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcaccttcc                                                                                                                                            >  2:898696/1‑149 (MQ=255)
          tcatcCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcaccttca                                                                                                                                            >  1:489693/1‑150 (MQ=255)
           catcCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcaccttcac                                                                                                                                           <  2:1695689/150‑1 (MQ=255)
           catcCAGCTGGTGCGCGAGTGCTTCAACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCAACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcaccttcac                                                                                                                                           >  2:170455/1‑150 (MQ=255)
             tcCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcaccttcacca                                                                                                                                         <  2:655062/150‑1 (MQ=255)
                aGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcaccttcaccacct                                                                                                                                      >  1:82733/1‑150 (MQ=255)
                      tGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGTCGCCGacagtctt                                                                                                                                             <  1:1172989/137‑6 (MQ=255)
                      tGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGTCGCCGacagtctt                                                                                                                                             >  2:1172989/1‑132 (MQ=255)
                      tGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGTCGCCGacagtctt                                                                                                                                             >  2:1172923/1‑132 (MQ=255)
                      tGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGTCGCCGacagtctt                                                                                                                                             <  1:1172923/137‑6 (MQ=255)
                        ggcgAGTGCGTCGACCGGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCCCCTTCACCTgg                                                                                                                              <  2:91569/149‑1 (MQ=255)
                        cgcgAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTgg                                                                                                                              <  1:500446/150‑1 (MQ=255)
                         gcgAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTggg                                                                                                                             >  1:844200/1‑150 (MQ=255)
                         gcgAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTggg                                                                                                                             >  1:582960/1‑150 (MQ=255)
                           gAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAg                                                                                                                           >  2:481191/1‑150 (MQ=255)
                                   cGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCaccga                                                                                                                   >  1:980974/1‑150 (MQ=255)
                                           gTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGccc                                                                                                           >  2:329548/1‑150 (MQ=255)
                                               gTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTgg                                                                                                                              >  2:669798/1‑127 (MQ=255)
                                               gTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTgg                                                                                                                              <  1:669798/127‑1 (MQ=255)
                                                 aCGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGt                                                                                                     <  2:582960/150‑1 (MQ=255)
                                                   gCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAAACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTg                                                                                                   >  1:1445041/1‑150 (MQ=255)
                                                    ccATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGc                                                                                                  <  1:737679/150‑1 (MQ=255)
                                                     cATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCg                                                                                                 >  1:1340291/1‑150 (MQ=255)
                                                      aTCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGACCAGTCGTTGCGt                                                                                                >  1:1787047/1‑150 (MQ=255)
                                                       tCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGTAAACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTg                                                                                               >  1:1861421/1‑150 (MQ=255)
                                                       tCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGACACTTCGGCCGTACTCCGCAGCGGAAGACTGTCAGCGAAGACACCTTTCGGAAGGGCACCTGGGAGCGCACCGACCGCGCCCATTCGTTGCGTg                                                                                               >  1:292781/1‑150 (MQ=255)
                                                         accaccATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGAc                                                                                             <  2:82733/150‑1 (MQ=255)
                                                                 gCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCgg                                                                                     >  1:915483/1‑150 (MQ=255)
                                                                 gCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCgg                                                                                     >  1:920885/1‑150 (MQ=255)
                                                                    cGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCt                                                                                  <  2:920885/150‑1 (MQ=255)
                                                                    cGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCt                                                                                  <  2:915483/150‑1 (MQ=255)
                                                                     gACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGGCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTg                                                                                 >  1:986304/1‑150 (MQ=255)
                                                                       ccTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTa                                                                               <  2:1329225/150‑1 (MQ=255)
                                                                        ctgctgCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGGACCGACGTGCTGCACCTTCACCTTGGAGTGACCCGAACGCACCAAATCAATGTGTGACGCTGCCGGCCTGTaa                                                                              >  1:1469168/1‑150 (MQ=255)
                                                                        ctgctgCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTaa                                                                              <  2:489693/150‑1 (MQ=255)
                                                                        ctgctgCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTaa                                                                              >  2:84802/1‑150 (MQ=255)
                                                                         tgctgcACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAg                                                                             <  2:55249/150‑1 (MQ=255)
                                                                             gcACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTc                                                                         <  2:446179/150‑1 (MQ=255)
                                                                                 gtgATGTATCAGGAAGCCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGc                                                                     >  2:472571/3‑150 (MQ=255)
                                                                                 ccGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGc                                                                     >  1:660817/1‑150 (MQ=255)
                                                                                  cGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCa                                                                    >  1:1909466/1‑150 (MQ=255)
                                                                                    aTGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAgc                                                                  >  2:229945/1‑150 (MQ=255)
                                                                                        aTCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCa                                                              >  1:1526087/1‑150 (MQ=255)
                                                                                        aTCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCa                                                              >  1:551320/1‑150 (MQ=255)
                                                                                        aTCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAACACTCTCCGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGGGCCa                                                              >  2:420627/1‑150 (MQ=255)
                                                                                         tCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAg                                                             >  2:533977/1‑150 (MQ=255)
                                                                                         tCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAg                                                             >  2:570493/1‑150 (MQ=255)
                                                                                            ggAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTg                                                                      >  2:665593/1‑138 (MQ=255)
                                                                                            ggAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTg                                                                      <  1:665593/138‑1 (MQ=255)
                                                                                            ggAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACg                                                          >  2:1122590/1‑150 (MQ=255)
                                                                                             gAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGa                                                         >  1:380252/1‑150 (MQ=255)
                                                                                               aaCCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGaaa                                                       <  1:1487033/150‑1 (MQ=255)
                                                                                                 ccGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGc                                                     <  2:980974/150‑1 (MQ=255)
                                                                                                   gCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGAAGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGccc                                                   >  1:347242/1‑150 (MQ=255)
                                                                                                    cTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAAACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGcccc                                                  <  2:1445041/150‑1 (MQ=255)
                                                                                                         tACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCAGCACCTTCACCTGGGAGCGCACCGACcgc                                                                                                               <  2:995907/84‑1 (MQ=255)
                                                                                                         tACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACcgc                                                                                                               >  1:995907/1‑84 (MQ=255)
                                                                                                          gcGGCCACTTCGGCCGTCCTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGGTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGa                                            <  2:308956/149‑1 (MQ=255)
                                                                                                           cGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAg                                           <  1:481191/150‑1 (MQ=255)
                                                                                                              ccACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcacc                                                                                                                                                >  2:1431558/1‑46 (MQ=255)
                                                                                                              ccACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGAcacc                                                                                                                                                <  1:1431558/46‑1 (MQ=255)
                                                                                                                 cTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCg                                     <  1:84802/150‑1 (MQ=255)
                                                                                                                  ttCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGa                                            <  2:1118714/142‑1 (MQ=255)
                                                                                                                  ttCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGa                                            >  1:1118714/1‑142 (MQ=255)
                                                                                                                  ttCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCgg                                    >  1:996353/1‑150 (MQ=255)
                                                                                                                     ggCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                 <  1:1122590/150‑1 (MQ=255)
                                                                                                                      gCCGTGCTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGATCGCTCGGCAgg                                <  2:292781/150‑1 (MQ=255)
                                                                                                                      gCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                <  2:844200/150‑1 (MQ=255)
                                                                                                                       ccGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAggg                               <  1:472571/150‑1 (MQ=255)
                                                                                                                        cGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgggg                              <  1:555730/150‑1 (MQ=255)
                                                                                                                             tCCGCAGCAGAAGACTGTCGGCGACGACaccttcaccaccttcacc                                                                                                                                 <  2:1271255/46‑1 (MQ=255)
                                                                                                                             tCCGCAGCAGAAGACTGTCGGCGACGACaccttcaccaccttcacc                                                                                                                                 >  1:1271255/1‑46 (MQ=255)
                                                                                                                             tCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCt                                                                                  >  2:1508903/1‑93 (MQ=255)
                                                                                                                             tCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCt                                                                                  <  1:1508903/93‑1 (MQ=255)
                                                                                                                             tCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCGTCACCTGGGAGCGGACCGACAACACCCAGTCGTTGCGTCACGCTGCCGGCCTTTGAGCTTCCGGCAGCGCCAGGCGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttt                         >  2:688251/1‑150 (MQ=255)
                                                                                                                               cGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                       >  2:542694/1‑150 (MQ=255)
                                                                                                                                gcagcaGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGACCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtg                      >  2:1111528/1‑150 (MQ=255)
                                                                                                                                 cagcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                     <  2:996353/150‑1 (MQ=255)
                                                                                                                                   gcagAAGACTGTCGGCGACGACACCTTCAGCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGGGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                   <  2:705689/150‑1 (MQ=255)
                                                                                                                                      gaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCt                                                                                  >  2:1678423/1‑84 (MQ=255)
                                                                                                                                      gaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCt                                                                                  <  1:1678423/84‑1 (MQ=255)
                                                                                                                                         gaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTa             <  1:542694/150‑1 (MQ=255)
                                                                                                                                           cTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACg           <  1:824105/150‑1 (MQ=255)
                                                                                                                                            tGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTcc                                                                        <  2:1667804/88‑1 (MQ=255)
                                                                                                                                            tGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTcc                                                                        >  1:1667804/1‑88 (MQ=255)
                                                                                                                                            tGTCGGCGAACACACCTCCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACgg          >  2:1019553/1‑150 (MQ=255)
                                                                                                                                                 gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAttt     >  2:853182/1‑150 (MQ=255)
                                                                                                                                                  cgacgaCACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAtttt    >  2:1276955/1‑150 (MQ=255)
                                                                                                                                                    acgacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCa  <  1:853182/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACGACAAGATCATCCAGCTGGTGCGCGAGTGCTTCGACCTGCGTCCGTACGCCATCACCACCATGCTCGACCTGCTGCACCCGATGTATCAGGAAACCGCTGCCTACGGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCA  >  NC_002947/5659977‑5660274

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