Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,659,739 T→A 5.9% T244T (ACT→ACA) ‡ metK → methionine adenosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,659,7390TA 5.9% 256.1 / 6.0 85T244T (ACT→ACA) ‡metKmethionine adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (43/36);  new base A (2/3);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

TTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGA  >  NC_002947/5659591‑5659888
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                                              >  2:1451099/1‑120 (MQ=255)
                                    tGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCAGTCAGGCCGCagt                                                                                                                                              <  1:1451099/122‑3 (MQ=255)
                                     gATCGTCAAGCACACGCTGCCGGCCGAACTACTCCACAAAGGCACCCAGTACACCATCACCGAGAGCGGCAACATCAACATCGGGTTCGTGGTGGTGGACGTGGGCCTGGCTGTCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGc                                                                                                                 >  2:1654259/1‑150 (MQ=255)
                                       tCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGccc                                                                                                               <  1:441453/150‑1 (MQ=255)
                                          tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                                                        <  2:522302/122‑1 (MQ=255)
                                          tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                                                        >  1:522302/1‑122 (MQ=255)
                                          tCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTc                                                                                                            >  1:1143386/1‑150 (MQ=255)
                                            aaGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGat                                                                                                                                            <  2:1871872/116‑1 (MQ=255)
                                            aaGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGat                                                                                                                                            >  1:1871872/1‑116 (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                                                                                                                                                    <  1:1309433/104‑1 (MQ=255)
                                                acacGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGc                                                                                                                                                    >  2:1309433/1‑104 (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                                                                                                       <  2:1857212/111‑1 (MQ=255)
                                                                                      tACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACggg                                                                                                       >  1:1857212/1‑111 (MQ=255)
                                                                                         cacaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAccc                                                                             >  2:558076/1‑134 (MQ=255)
                                                                                         cacaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAccc                                                                             <  1:558076/134‑1 (MQ=255)
                                                                                          acaTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGtt                                                            <  1:68939/150‑1 (MQ=255)
                                                                                            aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                               >  2:216672/1‑129 (MQ=255)
                                                                                            aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                                                               <  1:216672/129‑1 (MQ=255)
                                                                                            aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTcc                                                          >  2:1365220/1‑150 (MQ=255)
                                                                                            aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTcc                                                          <  1:1365220/150‑1 (MQ=255)
                                                                                            aTCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTcc                                                          >  1:1593184/1‑150 (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)
                                                                                                                                                    tGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGa  >  2:1564484/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTACTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGA  >  NC_002947/5659591‑5659888

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