Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 867,998 A→G 5.4% intergenic (+93/+84) 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,9980AG 5.4% 279.9 / 4.3 96intergenic (+93/+84)PP_0750/mqo‑Ihypothetical protein/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (50/38);  new base G (3/2);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

GCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCA  >  NC_002947/867849‑868145
                                                                                                                                                     |                                                                                                                                                   
ggcGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCGCAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGGGCGGGCa                                                                                                                                                     <  2:33558/148‑1 (MQ=255)
gCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCa                                                                                                                                                     >  2:1149764/1‑150 (MQ=255)
 ttaagccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaacgccctcttcgcgggcagGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCAt                                                                                                                                                    >  2:1539096‑M1/89‑150 (MQ=255)
 ccGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCAt                                                                                                                                                    >  2:1481817/1‑150 (MQ=255)
  cGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATg                                                                                                                                                   <  1:1149764/150‑1 (MQ=255)
    agccggtctgaaacgcaaaaacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgccctcttcgcgggcgcGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGcc                                                                                                                                                 >  1:243913‑M1/86‑150 (MQ=255)
    aCCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGTGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGcc                                                                                                                                                 >  1:1594136/1‑150 (MQ=255)
      cAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCcg                                                                                                                                               >  1:717901/1‑150 (MQ=255)
            aaTGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCCGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCTGGCATTACCGAACCTGTGGGAGCGGGCATGGCCGCGAaga                                                                                                                                         >  1:40294/1‑150 (MQ=255)
              tGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAgg                                                                                                                                       >  1:1285952/1‑150 (MQ=255)
                ccccTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGc                                                                                                                                     >  1:482657/1‑150 (MQ=255)
                 cccTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCg                                                                                                                                    >  2:987433/1‑150 (MQ=255)
                 cccTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCg                                                                                                                                    <  1:1272591/150‑1 (MQ=255)
                 cccTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCg                                                                                                                                    >  2:467221/1‑150 (MQ=255)
                      aGTACATCCAGATCATAGATAAAGCGCGAGGCGGATCTGTAGCCAGTGCCGGCCTCTTCGCGCGCGCGCACGCTCCCCCAGGAGTCGCGCCGCACCTGCGGGCATTACCCAACCTGGGGGAGCAGCCATGCCCGCGAAGAGGGCGGTACa                                                                                                                               <  2:10194/150‑1 (MQ=255)
                       aacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgccctcttcgcgggcgcGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAagag                                                                                                                                        <  1:582370‑M1/74‑1 (MQ=255)
                       aacgcccctcggggcgttttttgtgcctggcgcttgatcagcctgtaccgccctcttcgcgggcgcGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAagag                                                                                                                                        >  2:582370‑M1/67‑140 (MQ=255)
                        tACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAgg                                                                                                                             >  2:1596744/1‑150 (MQ=255)
                          cATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGCACAGGCt                                                                                                                           <  2:219687/150‑1 (MQ=255)
                           aTCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCCATACCTAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTg                                                                                                                          >  1:514651/1‑150 (MQ=255)
                            tCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGa                                                                                                                         <  1:1481817/150‑1 (MQ=255)
                            tCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGa                                                                                                                         <  2:859252/150‑1 (MQ=255)
                            tCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGa                                                                                                                         >  2:796506/1‑150 (MQ=255)
                            tCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGa                                                                                                                         <  1:28421/150‑1 (MQ=255)
                            tCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGa                                                                                                                         >  2:28421/1‑150 (MQ=255)
                                  gggcgttttttgtgcctggcgcttgatcagactgtaccggcctcttcgcgggcgcGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAgc                                                                                                                   >  2:1552154‑M1/56‑150 (MQ=255)
                                    atcGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAgcgc                                                                                                                 <  1:987433/150‑1 (MQ=255)
                                    atcGATAAAGCCCGTGGCTGATCTGTAGCCAGGGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAgcgc                                                                                                                 >  2:112752/1‑150 (MQ=255)
                                     tcGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGACCCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAgg                                                                                                                             <  2:1056234/137‑1 (MQ=255)
                                     tcGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAgg                                                                                                                             >  1:1056234/1‑137 (MQ=255)
                                      cGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCa                                                                                                               <  1:796506/150‑1 (MQ=255)
                                       gcttaaGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAg                                                                                                              >  2:1264347/5‑150 (MQ=255)
                                       gATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAg                                                                                                              <  1:1596744/150‑1 (MQ=255)
                                       gATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAg                                                                                                              >  1:1148066/1‑150 (MQ=255)
                                         tAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGc                                                                                                            <  2:717901/150‑1 (MQ=255)
                                                cgtcggtGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACaaaaa                                                                                                     >  2:406678/7‑150 (MQ=255)
                                                 gtcgCTGATCTGTAGGCAGTGCCGGCCTCTTCGCGCGCGCGCCCGCTCACACAGGTGTCGCGCCGCACGTGCGGGCATTACCGAAACTGTGGGAGCGGGCACGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACaaaaaa                                                                                                    <  1:1301866/147‑1 (MQ=255)
                                                 gTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGAACaaaaaa                                                                                                    >  2:849451/1‑150 (MQ=255)
                                                  tGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAAc                                                                                                   <  1:960570/150‑1 (MQ=255)
                                                      tGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATc                                                                                                                       <  1:567055/126‑1 (MQ=255)
                                                      tGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATc                                                                                                                       >  2:567055/1‑126 (MQ=255)
                                                      tGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGccc                                                                                               >  2:7443/1‑150 (MQ=255)
                                                      gcttgatcagcctgtaccgccctcttcgcgggcatgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGCCGGTACAGTCTGATCAAGCGCCAGGCACAAAAAACGccc                                                                                               >  2:1552154‑M2/96‑150 (MQ=255)
                                                       gATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAg                                                                                                                              >  1:762147/1‑118 (MQ=255)
                                                       gATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAg                                                                                                                              <  2:762147/118‑1 (MQ=255)
                                                          cTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAg                                                                                           <  2:482657/150‑1 (MQ=255)
                                                            gTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAggg                                                                                         >  1:1085998/1‑150 (MQ=255)
                                                             tAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCAt                                                                                                                                                    >  1:524123/1‑90 (MQ=255)
                                                             tAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCAt                                                                                                                                                    <  2:524123/90‑1 (MQ=255)
                                                               gCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAggg                                                                                         <  2:105255/147‑1 (MQ=255)
                                                               gCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAggg                                                                                         >  1:105255/1‑147 (MQ=255)
                                                                ccAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGt                                                                                     >  1:1589174/1‑150 (MQ=255)
                                                                 cAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGtt                                                                                    >  1:939282/1‑150 (MQ=255)
                                                                   gTGCCGGCCTCTTCGCGGGCGCGCCCGCCCCCACAGTCGTTGCGCCGACCCTGCGGGCATGACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGtttt                                                                                  <  2:514651/150‑1 (MQ=38)
                                                                      ccGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGc                                                                               <  1:467221/150‑1 (MQ=255)
                                                                        ggCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCTCGGGGCGtttt                                                                                  >  2:715552/1‑145 (MQ=255)
                                                                        ggCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCTCGGGGCGtttt                                                                                  <  1:715552/145‑1 (MQ=255)
                                                                        ggCCTCGTCGCGGGCTCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGt                                                                             <  2:40294/150‑1 (MQ=255)
                                                                         gCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGtt                                                                            >  1:871840/1‑150 (MQ=255)
                                                                         gCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGtt                                                                            <  2:392763/150‑1 (MQ=255)
                                                                           ctcttcgcgggcatgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGtt                                                                                    >  2:582370‑M2/75‑140 (MQ=255)
                                                                           ctcttcgcgggcatgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGtt                                                                                    <  1:582370‑M2/66‑1 (MQ=255)
                                                                           ctctTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTc                                                                          <  2:1285952/150‑1 (MQ=255)
                                                                                    ggcatgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCt                                                                 >  1:243913‑M2/66‑150 (MQ=255)
                                                                                       atgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcCTGCCCGCGAAGAGGGCGTTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTaa                                                              >  2:1539096‑M2/63‑150 (MQ=255)
                                                                                         gcCCGCACCCACAGGAAGCGCAGCGCAGCTGCGGGCATTACCGCACCACTGCGCACCAGCTCGCCGGATTTTTGGTCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGc                                                            <  1:1264347/150‑1 (MQ=255)
                                                                                           ccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACaaaaaa                                                                                                    <  1:399864‑M2/50‑1 (MQ=255)
                                                                                           ccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACaaaaaa                                                                                                    >  2:399864‑M2/59‑108 (MQ=255)
                                                                                            CGCTCCCAGAGGAGTCGCGCCGAACCTGGGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAACAGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTa                                                         >  2:812059‑M1/1‑150 (MQ=255)
                                                                                            CGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTa                                                         >  2:1594136‑M1/1‑150 (MQ=255)
                                                                                                 ccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttttt                                                                                 <  2:1300575‑M2/69‑1 (MQ=255)
                                                                                                 ccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttttt                                                                                 >  1:1300575‑M2/53‑121 (MQ=255)
                                                                                                  CACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                   >  2:871840‑M1/1‑150 (MQ=255)
                                                                                                   acaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                  <  2:243913‑M2/100‑1 (MQ=255)
                                                                                                   ACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                  <  1:1292834‑M1/150‑1 (MQ=255)
                                                                                                   ACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                  >  1:7443‑M1/1‑150 (MQ=255)
                                                                                                    caggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAgg                                                                                                                                       >  1:1066527‑M2/50‑64 (MQ=255)
                                                                                                    caggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAgg                                                                                                                                       <  2:1066527‑M2/15‑1 (MQ=255)
                                                                                                    CAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttttt                                                                                 <  1:840732‑M1/118‑1 (MQ=255)
                                                                                                    CAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttttt                                                                                 >  2:840732‑M1/1‑118 (MQ=255)
                                                                                                        aGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGgccgc                                             <  2:1085998/150‑1 (MQ=255)
                                                                                                            gcgcCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCggg                                         >  1:622995/1‑150 (MQ=255)
                                                                                                            ccgccccACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCggg                                         <  1:674215/149‑1 (MQ=255)
                                                                                                              gccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGgct                                       >  1:616607/1‑150 (MQ=255)
                                                                                                              gccgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAACCGTAGTACTGTGGCGCCGGgct                                       >  2:298575/1‑150 (MQ=255)
                                                                                                                cgcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCTGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGgctgc                                     >  2:52242/1‑150 (MQ=255)
                                                                                                                    ccTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctc                                 >  1:1227546/1‑150 (MQ=255)
                                                                                                                      ggCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTg                               <  2:1589174/149‑1 (MQ=255)
                                                                                                                               tACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCa                      >  2:574380/1‑150 (MQ=255)
                                                                                                                                 ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                    >  2:450122/1‑150 (MQ=255)
                                                                                                                                 ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                    <  2:1292834/150‑1 (MQ=255)
                                                                                                                                 ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                    <  1:619759/150‑1 (MQ=255)
                                                                                                                                  cGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAg                           <  1:206466/142‑1 (MQ=255)
                                                                                                                                  cGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAg                           >  2:206466/1‑142 (MQ=255)
                                                                                                                                  cGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAggg                   >  2:1163686/1‑150 (MQ=255)
                                                                                                                                         gtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgtc            <  1:1596085‑M2/138‑1 (MQ=255)
                                                                                                                                                 gggcCTGCCCGCGAAGAGGGCGTTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCgg    <  1:1539096‑M2/146‑1 (MQ=255)
                                                                                                                                                   gcATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCa  >  2:1640469‑M2/3‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCTGATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCA  >  NC_002947/867849‑868145

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: