Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,659,737 A→T 5.3% T244S (ACT→TCT) ‡ metK → methionine adenosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,659,7370AT 5.3% 358.3 / 6.9 113T244S (ACT→TCT) ‡metKmethionine adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (60/47);  new base T (2/4);  total (62/51)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 113 total aligned reads displayed.

GTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGC  >  NC_002947/5659590‑5659885
                                                                                                                                                   |                                                                                                                                                    
ttttCTCAGAAATACCTGCAAGAAGCCGTGCTGGCGGTGATCGTCAAGTACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACt                                                                                                                                                    <  2:2202160/149‑1 (MQ=255)
       aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCa                                                                                                                                             <  1:557063/150‑1 (MQ=255)
       aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCa                                                                                                                                             <  1:2065767/150‑1 (MQ=255)
        gAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCaa                                                                                                                                            <  1:1433214/150‑1 (MQ=255)
        aaaaGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCag                                                                                                                                            <  1:1096751/149‑2 (MQ=255)
         aaaGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAg                                                                                                                                           >  2:937069/1‑150 (MQ=255)
         aaaGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAg                                                                                                                                           >  1:646630/1‑150 (MQ=255)
             aCCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatca                                                                                                                                       >  2:681090/1‑150 (MQ=255)
               cTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcctc                                                                                                                                     <  1:2153992/150‑1 (MQ=255)
               cTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                                                     >  1:1756905/1‑150 (MQ=255)
               cTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGACGCACTACTGAACAAAGGCACCCAGGAGCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                                                     >  1:296332/1‑150 (MQ=255)
                tGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcg                                                                                                                                    >  2:777846/1‑150 (MQ=255)
                  cAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgtc                                                                                                                                  <  2:754627/150‑1 (MQ=255)
                            tGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACg                                                                                                                        <  2:1497208/150‑1 (MQ=255)
                              aTGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGt                                                                                                                      <  1:937069/150‑1 (MQ=255)
                                 gAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCaa                                                                                                                                            <  1:3448/125‑1 (MQ=255)
                                 gAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCaa                                                                                                                                            >  2:3448/1‑125 (MQ=255)
                                     tGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATgg                                                                                                               <  2:296332/150‑1 (MQ=255)
                                        tCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCTCCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGccc                                                                                                            <  1:2200443/150‑1 (MQ=255)
                                             aatcacaCGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCa                                                                                                                                             <  1:127916/109‑1 (MQ=255)
                                             aaGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCAc                                                                                                       >  2:1735211/1‑150 (MQ=255)
                                             aaGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCAc                                                                                                       >  2:1336477/1‑150 (MQ=255)
                                              aGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACg                                                                                                      >  1:788212/1‑150 (MQ=255)
                                               gCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgg                                                                                                     <  2:343018/150‑1 (MQ=255)
                                                 acacGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCCCAAGATCATCGTCGACTCCTTCGGTGGCATGGCCCGCCACgggg                                                                                                   <  2:690108/150‑1 (MQ=255)
                                                   acGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGt                                                                                                 >  2:1906188/1‑150 (MQ=255)
                                                    cGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGGCTCCTACGGTGGCATGGCCCGTCACGGGGGTg                                                                                                <  2:1573881/150‑1 (MQ=255)
                                                      cTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgc                                                                                              <  2:2207324/150‑1 (MQ=255)
                                                              cgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACg                                                                                                                        <  1:243382/116‑1 (MQ=255)
                                                                 actactGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCgg                                                                                   >  2:335456/1‑150 (MQ=255)
                                                                     ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTg                                                                                                                     <  1:917525/112‑1 (MQ=255)
                                                                     ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTg                                                                                                                     >  2:917525/1‑112 (MQ=255)
                                                                     ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgggg                                                                                                   <  1:241119/130‑1 (MQ=255)
                                                                     ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgggg                                                                                                   >  2:241119/1‑130 (MQ=255)
                                                                      tGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCa                                                                                                                                             >  1:2230659/1‑87 (MQ=255)
                                                                      tGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCa                                                                                                                                             <  2:2230659/87‑1 (MQ=255)
                                                                      tGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAgg                                                                              >  1:466390/1‑150 (MQ=255)
                                                                      tGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGAGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAgg                                                                              >  1:62075/1‑150 (MQ=255)
                                                                        cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCa                                                                                                                  >  1:1670551/1‑112 (MQ=255)
                                                                        cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCa                                                                                                                  <  2:1670551/112‑1 (MQ=255)
                                                                        cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgcgc                                                                                            >  1:1751620/1‑134 (MQ=255)
                                                                        cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgcgc                                                                                            <  2:1751620/134‑1 (MQ=255)
                                                                         acaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGCGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAcc                                                                           >  2:1165638/1‑150 (MQ=255)
                                                                                cACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCac                                                                    <  2:646630/150‑2 (MQ=255)
                                                                                cACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCaa                                                                    <  1:2270685/150‑1 (MQ=255)
                                                                                 aCCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAg                                                                   <  1:777846/150‑1 (MQ=255)
                                                                                 aCCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAg                                                                   >  1:1427783/1‑150 (MQ=255)
                                                                                   ccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGt                                                                 <  2:1427783/150‑1 (MQ=255)
                                                                                   ccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGCGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAACGACCCGTCCAAGGt                                                                 >  1:2054135/1‑150 (MQ=255)
                                                                                    cAGTACCACATCAACCCGACCGGCAACTTCAACACCGGTGGCGCGGTGGGTGACGGCGGCGTGACTGGCCGCAAGATCATCGTCGTCTTCGACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTTTCCGTCAAGGACCCGTCCAAGGTc                                                                <  2:502148/150‑1 (MQ=255)
                                                                                     aGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCg                                                               >  2:903186/1‑150 (MQ=255)
                                                                                     aGTACCACATCAACCCGACCGGCAAATTCAGCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGGAAGAACATAGTCGACTCCTACGGCAGCATGGCCCGTCACGGGGGTGGCGCGTTTTCCGGCAAAGATCCCTGCAAgattg                                                               >  2:1163260/1‑146 (MQ=255)
                                                                                       tACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCACGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                             >  1:359480/1‑150 (MQ=255)
                                                                                       tACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTTGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                             <  2:1463138/150‑1 (MQ=255)
                                                                                          cacaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATGATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGGCCGt                                                          >  1:549495/1‑150 (MQ=255)
                                                                                          cacaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                             <  2:1704914/147‑1 (MQ=255)
                                                                                          cacaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                                                             >  1:1704914/1‑147 (MQ=255)
                                                                                            caTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAgg                                                                  <  2:2345628/140‑1 (MQ=255)
                                                                                            caTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAgg                                                                  >  1:2345628/1‑140 (MQ=255)
                                                                                            caTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTc                                                        >  1:1982760/1‑150 (MQ=255)
                                                                                             aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGt                                                                                                                      <  1:1904841/87‑1 (MQ=255)
                                                                                             aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGt                                                                                                                      >  2:1904933/1‑87 (MQ=255)
                                                                                             aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGt                                                                                                                      >  2:1904841/1‑87 (MQ=255)
                                                                                             aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAccc                                                                          <  2:1134832/131‑1 (MQ=255)
                                                                                             aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAccc                                                                          >  1:1134832/1‑131 (MQ=255)
                                                                                                aaCCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGc                                                                                                              >  1:777332/1‑92 (MQ=255)
                                                                                                aaCCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGc                                                                                                              <  2:777332/92‑1 (MQ=255)
                                                                                                 aCCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                            >  2:1679979/1‑125 (MQ=255)
                                                                                                 aCCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                            <  1:1679979/125‑1 (MQ=255)
                                                                                                 aCCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTg                                                   >  1:1938032/1‑150 (MQ=255)
                                                                                                  cccGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCAATAGGTCCCATGGCCAGTCACGGGGGTTGCGCGTTCTCAGGACCGGACGCGTCCAAGGTCGACCGTTCCGCTGc                                                  >  1:1252029/1‑150 (MQ=255)
                                                                                                      aCCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTAc                                              <  1:1906188/150‑1 (MQ=255)
                                                                                                        cGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgc                                            >  1:1923719/1‑150 (MQ=255)
                                                                                                           cAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGc                                                  <  1:2153432/141‑1 (MQ=255)
                                                                                                            aaCTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccggc                                        >  2:1206572/1‑150 (MQ=255)
                                                                                                              cTTCATCATCGTTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccggccg                                      >  2:601180/1‑150 (MQ=255)
                                                                                                               ttCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGc                                     >  1:1704847/1‑150 (MQ=255)
                                                                                                               ttCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGc                                     >  2:165364/1‑150 (MQ=255)
                                                                                                                  atcatcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCg                                                                         <  1:780534/111‑1 (MQ=255)
                                                                                                                  atcatcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCg                                                                         >  2:780534/1‑111 (MQ=255)
                                                                                                                   tcatcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACg                                 >  2:1937162/1‑150 (MQ=255)
                                                                                                                    catcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGt                                <  1:601180/150‑1 (MQ=255)
                                                                                                                      tcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCTGCTGCCTACGCCGGCCGCTACGTgg                              >  2:253831/1‑150 (MQ=255)
                                                                                                                          tGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGGAAGGACCCGCCCAATGCCGACCGTTCCTGTTCCTACGCCGGCCGATCCGTGGCaaa                          >  2:397998/1‑150 (MQ=255)
                                                                                                                          tGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCaa                          <  1:397998/150‑1 (MQ=255)
                                                                                                                          tGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCaa                          <  2:1756905/150‑1 (MQ=255)
                                                                                                                              ccGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAAc                      >  1:992211/1‑150 (MQ=255)
                                                                                                                                 gtgggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtc                   <  2:466390/150‑1 (MQ=255)
                                                                                                                                  tgggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtcg                  <  1:335456/150‑1 (MQ=255)
                                                                                                                                      tgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTcgc              >  2:451215/1‑150 (MQ=255)
                                                                                                                                      tgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTcgc              >  1:2220733/1‑150 (MQ=255)
                                                                                                                                       gACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTcgcc             >  1:1009486/1‑150 (MQ=255)
                                                                                                                                       gACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTcgcc             >  1:842652/1‑150 (MQ=255)
                                                                                                                                           gCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGccgccg         <  2:2220733/150‑1 (MQ=255)
                                                                                                                                             ggCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGt       >  2:2344646/1‑150 (MQ=255)
                                                                                                                                               ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCaa                                                                                <  2:576753/75‑1 (MQ=255)
                                                                                                                                                cTGACTGTCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTg    >  2:408111/1‑150 (MQ=255)
                                                                                                                                                cTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtc                   <  1:1687000/135‑1 (MQ=255)
                                                                                                                                                 tGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTgg   >  1:1412559/1‑150 (MQ=255)
                                                                                                                                                  gACTGGCCGCAAGATCATCATCGAATCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTCGGc  >  2:312463/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGC  >  NC_002947/5659590‑5659885

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