Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 867,936 G→A 9.3% intergenic (+31/+146) PP_0750 → / ← mqo‑I hypothetical protein/malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947867,9360GA 9.3% 631.3 / 40.3 221intergenic (+31/+146)PP_0750/mqo‑Ihypothetical protein/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (92/106);  new base A (15/5);  total (107/114)
Fisher's exact test for biased strand distribution p-value = 1.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 229 total aligned reads displayed.

GAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAA  >  NC_002947/867789‑868085
                                                                                                                                                   |                                                                                                                                                      
 aaGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGC                                                                                                                                                     <  1:638672‑M1/150‑1 (MQ=255)
 aaGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCACTTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGC                                                                                                                                                     <  1:1272662‑M1/150‑1 (MQ=255)
  aGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCC                                                                                                                                                    >  2:5190678‑M1/1‑150 (MQ=255)
     tGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTAGCGGGCGCGCCCcc                                                                                                                                                 >  1:778661‑M1/1‑148 (MQ=255)
        cctggacctggaCGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCccacc                                                                                                                                              >  2:272730‑M1/1‑145 (MQ=255)
        cctggacctggaCGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCgaa                                                                                                                                              <  2:647139‑M1/150‑4 (MQ=255)
        cctggacctggaCGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCC                                                                                                                                              <  2:4128034‑M1/150‑1 (MQ=255)
        cctggacctggaCGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCGCC                                                                                                                                              >  1:5552488‑M1/1‑150 (MQ=255)
         ctggacctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCgaag                                                                                                                                             <  1:1451159‑M1/150‑5 (MQ=255)
         ctggacctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCC                                                                                                                                             <  1:3515630‑M1/150‑1 (MQ=255)
          tggacctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCTCA                                                                                                                                            <  1:2544772‑M1/150‑1 (MQ=255)
          tggacctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCA                                                                                                                                            >  2:470107‑M1/1‑150 (MQ=255)
            gacctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACA                                                                                                                                          <  1:2934325‑M1/150‑1 (MQ=255)
            gacctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACA                                                                                                                                          <  1:564389‑M1/150‑1 (MQ=255)
            cggccctggacctggcagagagcagcccggcggccaagcactacgattaagccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatctgtagccagtgccggCCTCTTCGCGGGCGCGCCCGCTCCcaca                                                                                                                                          >  1:5393437‑M1/123‑150 (MQ=255)
             acctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAG                                                                                                                                         <  1:3549653‑M1/150‑1 (MQ=255)
               ctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGa                                                                                                                                       <  2:3854110/150‑1 (MQ=255)
               ctggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGa                                                                                                                                       <  1:221421/150‑1 (MQ=255)
                tggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCCCAGGgg                                                                                                                                      >  2:1957416‑M1/1‑148 (MQ=255)
                tggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTAGCGGGCGCGCCCCCGCCCCCAGGAg                                                                                                                                      >  2:622055/1‑150 (MQ=255)
                 ggacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATATGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGt                                                                                                                                     <  1:4994286/150‑1 (MQ=255)
                  gacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCA                                                                                                                                            >  2:3885576‑M1/1‑142 (MQ=255)
                  gacGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCA                                                                                                                                            <  1:3885576‑M1/142‑1 (MQ=255)
                   acGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTcg                                                                                                                                   <  2:5034311/150‑1 (MQ=255)
                    cGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTcgc                                                                                                                                  <  2:1806896/150‑1 (MQ=255)
                      gCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCgcgc                                                                                                                                <  2:342896/150‑1 (MQ=255)
                      gCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGAGTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCGGATCTGTAGCCAGTGCCGGCCTCTGCGCGGGCGCGCCCGCGCCCACAGGAGTCgcgc                                                                                                                                >  2:1647973/1‑150 (MQ=255)
                        aaGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGcgccg                                                                                                                              >  2:3625961/1‑150 (MQ=255)
                         aGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGcgccgc                                                                                                                             <  2:3585103/150‑1 (MQ=255)
                           gACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGTGCCCGCTCCCACAGGttcggtaatgccc                                                                                                                           >  1:279523‑M1/1‑137 (MQ=255)
                               aGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGc                                                                                                                       <  1:3339349/150‑1 (MQ=255)
                                gCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCg                                                                                                                      <  2:1178937/150‑1 (MQ=255)
                                 cACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCCCAGGAGTCGCGCCGCACCTGCgg                                                                                                                     <  1:2922603/150‑1 (MQ=255)
                                      cAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCGCCCCCAGGAGTCGCGCCGCACCAGCGGGCaat                                                                                                               >  2:1042990/1‑148 (MQ=255)
                                        cccggcggccaagcactacgcttaagccggtccgcaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgcCCACTTCGCGGGCAGGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTAc                                                                                                             >  1:1805957‑M1/95‑150 (MQ=255)
                                        aGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTAc                                                                                                             <  2:478665/150‑1 (MQ=255)
                                            cTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGaa                                                                                                         <  2:820622/150‑1 (MQ=255)
                                            cTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTC                                                                                                                                               >  2:3453109‑M1/1‑113 (MQ=255)
                                              tcgtcgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCa                                                                                                                            >  2:5055784/1‑130 (MQ=255)
                                                 tcgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgt                                                                                                    >  2:391394‑M2/1‑145 (MQ=255)
                                                 tcgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgt                                                                                                    >  2:3757443‑M2/1‑145 (MQ=255)
                                                 tcgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGttcggtaatgcccgcaggtgcggcgcgactcctgt                                                                                                     <  2:1534000‑M1/150‑36 (MQ=255)
                                                 caagcactacgcttaagccggtctgaaacgcataaacgcccctcggggcgtttattgtgcctggcgcttgatcagcctgtaccgcCCTCTTCGCGGGCATGCCCGCTCCCACAGGAGTCGCGCCGCACCTTCGGGCATTACCGAACCtgt                                                                                                    >  1:3210082‑M1/86‑150 (MQ=255)
                                                  cgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCTGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgtg                                                                                                   >  1:4544690‑M2/1‑144 (MQ=255)
                                                  cgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGTGCCCGCTCCCACAGGttcggtaatgcccgcaggtgcggcgcgactcctgtg                                                                                                    <  2:279523‑M1/150‑37 (MQ=255)
                                                  cgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGttcggtaatgcccgcaggtgcggcgcgactcctgtg                                                                                                    <  1:5276817‑M1/150‑37 (MQ=255)
                                                  cgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgtg                                                                                                   >  1:3650059‑M2/1‑144 (MQ=255)
                                                  cgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgtg                                                                                                   >  2:5199807‑M2/1‑144 (MQ=255)
                                                  cgAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGttcggtaatgcccgcaggtgcggcgcgactcctgtg                                                                                                    <  2:1384585‑M1/150‑37 (MQ=255)
                                                  aagcactacgcttaagccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgcCCTCTTCGCGGGCACGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCtgtg                                                                                                   >  1:3423178‑M1/85‑150 (MQ=255)
                                                   gAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgtgg                                                                                                  >  2:4647441‑M2/1‑143 (MQ=255)
                                                   gAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCACCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCCAAccggtgg                                                                                                  <  2:1237217‑M2/150‑8 (MQ=255)
                                                   agcactacgcttaagccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgcCCTCTTCGCGGGCATGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTgg                                                                                                  >  1:1497836‑M1/84‑150 (MQ=255)
                                                     gCAGATTGCCGACCCGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGttcggtaatgcccgcaggtgcggcgcgactcctgtggga                                                                                                 <  1:3659683‑M1/150‑40 (MQ=255)
                                                     cactacgcttaagccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagccagtgccggCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGa                                                                                                >  2:4319582‑M1/82‑150 (MQ=255)
                                                       ctccgcttacgccggtctgaaacgcaacaacgcccctcggggcgttttttgggcctggcgcttgatcagcctgtaccgcCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGc                                                                                              >  1:2353570‑M1/80‑150 (MQ=255)
                                                       aGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgtgggagc                                                                                              >  2:3523044‑M2/1‑139 (MQ=255)
                                                        gacgcttaagccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgcCCTCTTCGCGGGCATGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCg                                                                                             >  2:4463219‑M1/79‑150 (MQ=255)
                                                         aTTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgtgggagcgg                                                                                            >  2:2180426‑M2/1‑137 (MQ=255)
                                                          ttGCCGACCAGTGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgtgggagcggg                                                                                           >  1:2715176‑M2/1‑136 (MQ=255)
                                                              cGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCTGCTCCCCCAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAcctgggggagcgggcctg                                                                                       <  2:5359878‑M2/150‑19 (MQ=255)
                                                               aagccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagccagtgccggCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGc                                                                                      >  2:2943629‑M1/72‑150 (MQ=255)
                                                                      gcAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGaa                                                                               <  2:693869/150‑1 (MQ=255)
                                                                      gcAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGaa                                                                               <  1:470107/150‑1 (MQ=255)
                                                                      gcAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCCCCGCACCAGCCGGCATTACCGAAcctgtgggagcgggccagcccgcgaa                                                                               <  2:5318947‑M2/150‑27 (MQ=255)
                                                                      gcAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCCCTCCCCCCCGAGTCGCGCCGCACCTGGGGGCCGTACCGAACCCGTGGGAGCGGGCAAGCCcgcgga                                                                               >  2:5318944/1‑148 (MQ=255)
                                                                        aaTGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCcgctccca                                                                             <  1:2484093/150‑6 (MQ=255)
                                                                           gCCCCTGAGTACATCTAGATAATGGATAAAGCCCGTGGGTGATCTTTAGCCAGTGCCGGCATCTTCGGGGGCGCGCCCGCTCCCACCGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAggg                                                                          <  1:1588855/150‑1 (MQ=255)
                                                                             cgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgcCCTCTTCGCGGGCACGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCg                                                                        >  2:2452178‑M1/58‑150 (MQ=255)
                                                                             cccTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACCGGAGTCGCGCCGCACCGGCGGGCATTACCGAACCTGTGGGGGCGGGCAAGCCCCCGAAGAGGGCg                                                                        >  2:2754255/1‑150 (MQ=255)
                                                                              ccTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTAGCGGGCGCGCCCGCTCCCACACGACTAGCGCCGCCCCAGCCGGCAATACCCAACCGGTGGGAGCGGGCAAGCCCCCGAAGAGGGCgg                                                                       >  2:5277939/1‑150 (MQ=255)
                                                                               cTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCGTCGCGGGCGCGCCCCCACCCACACGGGTCGCGCCGCACCAGCGGGCATTACCGAACCGGTGGGAGCGGGCAAGCCCCCGACCAGGGCggg                                                                      >  2:4305716/1‑149 (MQ=255)
                                                                                  aGTACAGCGAGATCATCGATACAGCCCGGGGCTGCTCTGTAGCCCGGGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCCCCGGAGTCGCGCCGCACCTGCGGGCATGACCGACCCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACa                                                                   <  2:5380704/150‑1 (MQ=255)
                                                                                   gTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAg                                                                  <  2:3154980/150‑1 (MQ=255)
                                                                                       aTCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGGGGGAGCGGGCCTGCCCGCGAAGAGGGCGGTACAGGCTg                                                              >  2:4558877/1‑150 (MQ=255)
                                                                                                 tcGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGcc                                                    <  2:3237088/150‑1 (MQ=255)
                                                                                                          cccGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACaaa                                           <  1:5023494/150‑1 (MQ=255)
                                                                                                             gTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTAc                                                                                                             >  2:1454787/1‑81 (MQ=255)
                                                                                                                  ttatgtGTAGCCTGTGCGGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGAATTACCTAACCTGTGGGAGGGGGCTTGCCCGGGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGccc                                   <  1:309231/145‑1 (MQ=38)
                                                                                                                  tGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCCCGCCCCCAGGCGGCGCGCCGCACCAGCGGGCCATAACCGACCGGTGGGGGCCGGCCAGCCCCCGCAGAGGGCGGGGCAGGCGGATCAACCCCCCCGCACAAAAAAAGccc                                   >  2:4764496/1‑150 (MQ=255)
                                                                                                                     tCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGa                                <  2:2891148/150‑1 (MQ=255)
                                                                                                                        gTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGTTCAAGCGCCAGGCACAAAAAACGcccctcggg                             <  2:4926591/150‑6 (MQ=255)
                                                                                                                        gTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCCCGCCCACAGGAGTCGCCCCCCACCAGCGGGCAATACCCAAACGGGGGGGGCGGGCCAGCCCCCCGAGAGGGCGGTACCGGCGGATGAAACCCCACGCACAAAAAAACccccccggg                             >  2:3025691/1‑145 (MQ=255)
                                                                                                                        gTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCCCGCCCACACGGGTGGCCCCCCCCCCGCGGGCCATAACCCAACCGGGGGGGCGGGCCAGCCCGCCCCGGGGGCGGGACAGGCCGTAGCAGCCCCCGGCCCAAAAAACCCCCCCAggg                             >  2:2360805/1‑150 (MQ=255)
                                                                                                                        gTAGCCAGTGCCGGCCTCCTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCCCCGCACCAGCGGGCAATAACGAACCAGGGGGCGCGGGCATGCCCCCCAAGAGGGCGGTACAGGCGGAAAAAGCGCCCGGCACAAAAAACGCCCCCAggg                             >  2:3994105/1‑150 (MQ=255)
                                                                                                                              tgTGCGGGCCTGTGCGGGGGCGCGCCCGCTCCCCCCGGCGGCGCGCCGCACCTGCGGGCATTCCCGACCCGGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttt                       <  1:1214430/149‑1 (MQ=255)
                                                                                                                                ggCCGGCCTCTTCGCGGGCGCGCCCGCTCCCCCCGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttttt                     <  1:882415/149‑1 (MQ=255)
                                                                                                                                  ccGGCCTCGGCGCGGGCGCGCCCGCGCCCACCGGAGTCGCGCCGCACCTGCGGGCAGGACCGACCCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAACGCCCCTCGGGGCGttttttgt                  <  2:2841837/150‑4 (MQ=255)
                                                                                                                                      CCTCTGCGCGGCCGCGCCCGCCCCCACAGGAGTCGCGCCGCACCTGCGGGAATTACCGAACCTGTGGGCGCGGGCAGGCCGGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGttt               >  1:3932607‑M1/1‑150 (MQ=255)
                                                                                                                                       CTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTc              >  2:2500512‑M1/1‑150 (MQ=255)
                                                                                                                                       CTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTc              >  2:778661‑M1/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAAccggtgggagcgggcacgcccccgaagaggccggcactggctacagatcagccacgggctttatcgatgatctcgatgtactcaggggcatt             <  2:3250032‑M2/150‑93 (MQ=255)
                                                                                                                                          ttCGCGGGCGCGCCCGCGCCCACACGAGTCGCGCCGCACCAGCGGGCCATACCGAACCGGGGGGAGCGGGCATGCCCGCGCCGAGGGCGGTACAGGCTTATAAAGCGCCCGGCACAAAAAACGCCCCGAGGGGCGGTGTTGCGTTTgaga           >  2:1392622/1‑150 (MQ=255)
                                                                                                                                           TCGCGGGCGCGCCCGCTCCCACCGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGGGGGCGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGAc          >  2:2440664‑M1/1‑150 (MQ=255)
                                                                                                                                              gGGGCGCGCCCGCTCCCCCCGGAGTCGCGCCGCACCTGCGGGCATTACCGACCCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCgg       >  1:4171759‑M1/2‑150 (MQ=255)
                                                                                                                                              cgGGCCCGCCCCCTCCCCCCGCCGTCGCGCCGCCCCTGCGGGCATTACCGAACCGTTGGCCGCGCGAAGCCCCGCCAAGAGGGCGGTACAGGCTGATCACGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCgg       <  1:1957416/150‑1 (MQ=255)
                                                                                                                                               GGGCGCGCCCGCTCCCCCCGGAGGCGCGCCGCACCTGCGGGCATTCCCGAACCTGTGGGCGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGc      >  1:272730‑M1/1‑150 (MQ=255)
                                                                                                                                                   gcgcCCGCTCCCACCGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTaa  <  2:2116503/150‑1 (MQ=255)
                                                                                                                                                   gcgcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTaa  >  1:3484542/1‑150 (MQ=255)
                                                                                                                                                   gcgcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTaa  <  2:5552488/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAA  >  NC_002947/867789‑868085

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