Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,659,737 A→T 6.0% 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 6.0% 244.3 / 5.2 86T244S (ACT→TCT) ‡metKmethionine adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base T (2/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

AGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGC  >  NC_002947/5659588‑5659885
                                                                                                                                                     |                                                                                                                                                    
aGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGa                                                                                                                                                      <  1:1004290/150‑1 (MQ=255)
 ggTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGCACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGAc                                                                                                                                                     <  2:1271504/150‑1 (MQ=255)
   tttCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTg                                                                                                                                                   >  2:695725/1‑150 (MQ=255)
     tCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTAAATCTGTGGGTGACTGCGGCCTGACTGGc                                                                                                                                                 <  1:1483438/150‑1 (MQ=255)
         aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCa                                                                                                                                             >  1:1685977/1‑150 (MQ=255)
         aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCa                                                                                                                                             <  1:1845158/150‑1 (MQ=255)
         aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCa                                                                                                                                             <  1:1849574/150‑1 (MQ=255)
          gAAAGACCTGCAAGAAGCCGTGATGGGGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCaa                                                                                                                                            <  2:1464430/150‑1 (MQ=255)
          gAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCaa                                                                                                                                            <  2:1731938/150‑1 (MQ=255)
           aaaGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAg                                                                                                                                           <  1:1164314/150‑1 (MQ=255)
            aaGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGa                                                                                                                                          <  2:1137078/150‑1 (MQ=255)
                  tGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGGGACTGCGGCCTGACTGGCCGCAAGATCAtcg                                                                                                                                    >  1:772389/1‑150 (MQ=255)
                       gaagCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGAc                                                                                                                               >  2:1479289/1‑150 (MQ=255)
                           ccGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCt                                                                                                                           >  2:1670546/1‑150 (MQ=255)
                            ggTGAGGGAGGTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACATGATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTa                                                                                                                          <  2:162723/149‑1 (MQ=255)
                               gATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACgg                                                                                                                       <  1:589542/150‑1 (MQ=255)
                                aTGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGt                                                                                                                      <  1:1127730/150‑1 (MQ=255)
                                 tGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTg                                                                                                                     <  1:218259/150‑1 (MQ=255)
                                       tGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCagt                                                                                                                                           <  1:1451099/122‑3 (MQ=255)
                                       tGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCagt                                                                                                                                           >  2:1451099/1‑120 (MQ=255)
                                        gATCGTCAAGCACACGCTGCCGGCCGAACTACTCCACAAAGGCACCCAGTACACCATCACCGAGAGCGGCAACATCAACATCGGGTTCGTGGTGGTGGACGTGGGCCTGGCTGTCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGc                                                                                                              >  2:1654259/1‑150 (MQ=255)
                                          tCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGccc                                                                                                            <  1:441453/150‑1 (MQ=255)
                                             tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                                                     >  1:522302/1‑122 (MQ=255)
                                             tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                                                     <  2:522302/122‑1 (MQ=255)
                                             tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTc                                                                                                         >  1:1143386/1‑150 (MQ=255)
                                               aaGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGat                                                                                                                                         >  1:1871872/1‑116 (MQ=255)
                                               aaGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGat                                                                                                                                         <  2:1871872/116‑1 (MQ=255)
                                                aGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCACCCCGACCGTCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGCCTGGCCGCAAGATCATCGTCCAATCCTACGGTGGCATGGCTCGTCACg                                                                                                      <  2:342232/150‑1 (MQ=255)
                                                 gCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgg                                                                                                     >  1:1579820/1‑150 (MQ=255)
                                                 gCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgg                                                                                                     >  2:916523/1‑150 (MQ=255)
                                                 gCACACGCGGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGGCTCCTACGGTGGCATGGCCCGTCACgg                                                                                                     <  2:1668473/150‑1 (MQ=255)
                                                   acacGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGc                                                                                                                                                 >  2:1309433/1‑104 (MQ=255)
                                                   acacGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGc                                                                                                                                                 <  1:1309433/104‑1 (MQ=255)
                                                   acacGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCAAATTGACTGGCCGCAAGATCATCGTCGGCTCCTACGGTGGCATGGCCCGTCACgggg                                                                                                   <  1:807734/150‑1 (MQ=255)
                                                     acGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGa                                                                                                                                          >  2:457551/1‑109 (MQ=255)
                                                     acGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGa                                                                                                                                          <  1:457551/109‑1 (MQ=255)
                                                      cGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACCTCAACCCGACCGGCAACTTCATCATCGGTCGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTg                                                                                                <  2:183476/150‑1 (MQ=255)
                                                        cTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGATGCCCGGTTGTGTGACTGCGTTCTGACGGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgc                                                                                              <  1:1670546/150‑1 (MQ=255)
                                                          gccggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgcgc                                                                                            <  2:1148715/150‑1 (MQ=255)
                                                             ggccgAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGtt                                                                                         <  2:1914588/150‑1 (MQ=255)
                                                                 gAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTcc                                                                                     >  1:1308275/1‑150 (MQ=255)
                                                                    ctactGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGc                                                                                  >  1:1117681/1‑150 (MQ=255)
                                                                     tactGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGcc                                                                                 >  2:190387/1‑149 (MQ=255)
                                                                        tGCACAAAGGCACCCAGTACTACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAgg                                                                              >  2:68939/1‑150 (MQ=255)
                                                                        tGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAAATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGGGGCGCGTTCTCCGGCAAgg                                                                              >  1:881382/1‑150 (MQ=255)
                                                                         gCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGa                                                                             <  1:695725/150‑1 (MQ=255)
                                                                         gCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGATCTGCTGCaatta                                                                             >  2:1078970/1‑147 (MQ=255)
                                                                          cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATGGGAGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                            >  1:987938/1‑150 (MQ=255)
                                                                          cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTCGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                            >  2:326433/1‑150 (MQ=255)
                                                                           acaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAcc                                                                           >  1:1287621/1‑150 (MQ=255)
                                                                                     ccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGt                                                                 >  1:1329225/1‑150 (MQ=255)
                                                                                       aGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCg                                                               >  1:1748457/1‑150 (MQ=255)
                                                                                       aGTACCACAACAAGGCGAGGTGGAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTATGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCg                                                               <  2:987938/150‑1 (MQ=255)
                                                                                         tACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACggg                                                                                                    >  1:1857212/1‑111 (MQ=255)
                                                                                         tACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACggg                                                                                                    <  2:1857212/111‑1 (MQ=255)
                                                                                            cacaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAccc                                                                          <  1:558076/134‑1 (MQ=255)
                                                                                            cacaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAccc                                                                          >  2:558076/1‑134 (MQ=255)
                                                                                             acaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGtt                                                         <  1:68939/150‑1 (MQ=255)
                                                                                               aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                            >  2:216672/1‑129 (MQ=255)
                                                                                               aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                            <  1:216672/129‑1 (MQ=255)
                                                                                               aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTcc                                                       >  1:1593184/1‑150 (MQ=255)
                                                                                               aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTcc                                                       >  2:1365220/1‑150 (MQ=255)
                                                                                               aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTcc                                                       <  1:1365220/150‑1 (MQ=255)
                                                                                                 cAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGc                                                     <  2:1593184/150‑1 (MQ=255)
                                                                                                    cccGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGc                                                  >  2:631590/1‑150 (MQ=255)
                                                                                                     ccGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGcc                                                 >  2:163071/1‑150 (MQ=255)
                                                                                                     ccGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGGGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGcc                                                 >  1:835677/1‑150 (MQ=255)
                                                                                                           ggCAACTTCATCATCGGTGGCCCGGTGGGTGCCTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgcc                                           <  2:749762/150‑1 (MQ=255)
                                                                                                              aaCTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgccggc                                        <  1:163071/150‑1 (MQ=255)
                                                                                                                 ttCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGCGCAACGTTGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCAGCAACGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGc                                     >  2:737679/1‑150 (MQ=255)
                                                                                                                    atcatcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTAc                                  >  1:1287896/1‑150 (MQ=255)
                                                                                                                       atcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTg                               >  2:1838612/1‑150 (MQ=255)
                                                                                                                       atcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTg                               >  1:913274/1‑150 (MQ=255)
                                                                                                                       atcGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTg                               >  1:1125358/1‑150 (MQ=255)
                                                                                                                         cGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGc                             <  2:1548399/150‑1 (MQ=255)
                                                                                                                                ccGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGTACCCGTCCAAGGTGGACCGGTCCGCTGCCGACGCCGGCCGCTACGTGGCCAAGacc                      >  2:1350299/1‑148 (MQ=255)
                                                                                                                                    tgggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCATGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtcg                  >  2:999004/1‑150 (MQ=255)
                                                                                                                                    tgggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtcg                  <  1:1479289/150‑1 (MQ=255)
                                                                                                                                      ggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCTACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtcgtc                >  1:1097496/1‑150 (MQ=255)
                                                                                                                                        atcttGCGGCCAGTCAGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTcgc              >  1:655062/5‑150 (MQ=255)
                                                                                                                                            tGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTc                                                                >  2:221096/1‑96 (MQ=255)
                                                                                                                                            tGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTc                                                                <  1:221096/96‑1 (MQ=255)
                                                                                                                                            tGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCAGCCATTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCgccgcc          <  1:1078970/150‑1 (MQ=255)
                                                                                                                                            tGCGGCCTGACTGGACGCAAGATCATCGTCGACTCCTACGGTGGCATGGACCGTCACGGGGGTGGCGCGTTGTCAGGTAGGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCgccgcc          >  1:704388/1‑150 (MQ=255)
                                                                                                                                                 ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAg                                                                   <  2:689259/88‑1 (MQ=255)
                                                                                                                                                 ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAg                                                                   >  1:689259/1‑88 (MQ=255)
                                                                                                                                                    gACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGAGTGGCCAGTACTTCTACAAGGATACGTCTGAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGc  <  1:134370/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGC  >  NC_002947/5659588‑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: