Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 109,959 A→C 5.4% intergenic (+22/‑22) PP_0103 → / → ctaD cytochrome c oxidase subunit 2/cytochrome c oxidase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947109,9590AC 5.4% 217.9 / 2.6 73intergenic (+22/‑22)PP_0103/ctaDcytochrome c oxidase subunit 2/cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base C (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CAAGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGTT  >  NC_002947/109811‑110107
                                                                                                                                                    |                                                                                                                                                    
cAAGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaa                                                                                                                                                     <  1:1061365/150‑1 (MQ=255)
 aaGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaat                                                                                                                                                    >  1:979374/1‑150 (MQ=255)
         gTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCa                                                                                                                                            <  1:1426604/150‑1 (MQ=255)
            ggAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGga                                                                                                                                         <  2:1070765/150‑1 (MQ=255)
                 tCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAg                                                                                                                                    >  1:186490/1‑150 (MQ=255)
                   gACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGa                                                                                                                                  >  2:907461/1‑150 (MQ=255)
                    aCATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGAc                                                                                                                                 >  2:1446011/1‑150 (MQ=255)
                     cATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACa                                                                                                                                <  2:898715/150‑1 (MQ=255)
                      aTCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACCCGCGGCAATCAATGGATTGCAGGAGACAGGACAt                                                                                                                               >  2:755302/1‑150 (MQ=255)
                          ctgctgTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGt                                                                                                                           <  1:1446011/150‑1 (MQ=255)
                           tgctgTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACTTGGGGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTg                                                                                                                          >  1:1117736/1‑150 (MQ=255)
                                tAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGAGACATGGTGACGCCGAAGGACGTGCTGGAGAACTACCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTg                                                                                                                     >  2:874967/1‑150 (MQ=255)
                                    cACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcg                                                                                                                 >  1:1055315/1‑150 (MQ=255)
                                     aCCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcga                                                                                                                <  1:1476931/150‑1 (MQ=255)
                                     aCCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcga                                                                                                                >  1:511838/1‑150 (MQ=255)
                                       cTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgacg                                                                                                              >  1:610184/1‑150 (MQ=255)
                                                 aaTGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCa                                                                                                    <  2:610184/150‑1 (MQ=255)
                                                 aaTGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCa                                                                                                    >  1:559456/1‑150 (MQ=255)
                                                 aaTGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCa                                                                                                    >  1:1543226/1‑150 (MQ=255)
                                                  aTGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCAc                                                                                                   >  2:294310/1‑150 (MQ=255)
                                                           gCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGAc                                                                                          >  1:1472080/1‑150 (MQ=255)
                                                           gCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGAc                                                                                          >  2:347639/1‑150 (MQ=255)
                                                           gCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGAc                                                                                          >  1:833138/1‑150 (MQ=255)
                                                           gCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGAc                                                                                          >  2:1446397/1‑150 (MQ=255)
                                                           gCAACAACAAAGGCGACATGGTGAAGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGAc                                                                                          >  1:1041689/1‑150 (MQ=255)
                                                             aacaacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCa                                                                                        >  1:12075/1‑150 (MQ=255)
                                                             aacaacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCa                                                                                        >  1:743147/1‑150 (MQ=255)
                                                                aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATgc                                                                                     >  1:1371833/1‑150 (MQ=255)
                                                                  caaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATgcgc                                                                                   >  1:1411714/1‑150 (MQ=255)
                                                                       gCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAg                                                                                                                       <  1:1318669/109‑1 (MQ=255)
                                                                       gCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAg                                                                                                                       >  2:1318669/1‑109 (MQ=255)
                                                                               gTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCaa                                                                      <  1:105888/150‑1 (MQ=255)
                                                                                    gCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgacga                                                                                                             <  2:564254/106‑1 (MQ=255)
                                                                                    gCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgacga                                                                                                             >  1:564254/1‑106 (MQ=255)
                                                                                      cGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaat                                                                                                                                                    <  1:1294839/65‑1 (MQ=255)
                                                                                      cGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaat                                                                                                                                                    >  2:1294839/1‑65 (MQ=255)
                                                                                       gAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACAt                                                                                                                               <  2:1507865/85‑1 (MQ=255)
                                                                                       gAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACAt                                                                                                                               >  1:1507865/1‑85 (MQ=255)
                                                                                        aaGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGAt                                                             <  2:743147/150‑1 (MQ=255)
                                                                                            aCGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATgcgc                                                         <  1:571583/150‑1 (MQ=255)
                                                                                             cGTGCTGGCGCTTAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGGGACAGGACAGGAGTGCAGTGATCGACGACCACGCCCCCGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCt                                                        >  1:416383/1‑150 (MQ=255)
                                                                                             cGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCt                                                        <  1:347639/150‑1 (MQ=255)
                                                                                                 cTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGt                                                    <  1:294310/150‑1 (MQ=255)
                                                                                                  tGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTg                                                   <  2:1055315/150‑1 (MQ=255)
                                                                                                        tGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACCCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGacc                                             <  1:755302/150‑1 (MQ=255)
                                                                                                         gAAGCCGGTTGAAGTCAAGGGAACCGTGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGacca                                            <  2:1117736/150‑1 (MQ=255)
                                                                                                         gAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGc                                                                  <  2:1308309/128‑1 (MQ=255)
                                                                                                         gAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGc                                                                  >  1:1308309/1‑128 (MQ=255)
                                                                                                              aggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGTTGCTGACCaccaac                                       <  1:874967/150‑1 (MQ=255)
                                                                                                              aggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaac                                       <  2:186490/150‑1 (MQ=255)
                                                                                                               ggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaacc                                      >  1:483401/1‑150 (MQ=255)
                                                                                                                gcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaacca                                     <  2:956264/150‑1 (MQ=255)
                                                                                                                    aaGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATaaa                                 <  1:816951/150‑1 (MQ=255)
                                                                                                                    aaGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATaaa                                 <  2:1371833/150‑1 (MQ=255)
                                                                                                                      gCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGa                               <  2:1543226/150‑1 (MQ=255)
                                                                                                                         aaGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCa                                                                                                    >  1:929355/1‑78 (MQ=255)
                                                                                                                         aaGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCa                                                                                                    <  2:929355/78‑1 (MQ=255)
                                                                                                                         aaGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACAt                            <  2:1472080/150‑1 (MQ=255)
                                                                                                                         aaGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCCGACCACCAACCATAAAGACAt                            >  2:525231/1‑150 (MQ=255)
                                                                                                                          aGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCAt                                    >  1:713118/1‑141 (MQ=255)
                                                                                                                          aGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCAt                                    <  2:713118/141‑1 (MQ=255)
                                                                                                                                   ggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATg                  >  2:155326/1‑150 (MQ=255)
                                                                                                                                   ggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCAAGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATg                  >  1:1017794/1‑150 (MQ=255)
                                                                                                                                    gTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGt                 <  1:1446397/150‑1 (MQ=255)
                                                                                                                                         tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCaa                                                                      >  1:1033579/1‑92 (MQ=255)
                                                                                                                                         tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCaa                                                                      <  2:1033579/92‑1 (MQ=255)
                                                                                                                                         tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCa                       <  2:860723/139‑1 (MQ=255)
                                                                                                                                         tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCa                       >  1:860723/1‑139 (MQ=255)
                                                                                                                                          tcctgCAATCCATTGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGt                 >  2:1200304/5‑144 (MQ=255)
                                                                                                                                          tcctgCAATCCATTGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGt                 <  1:1200304/140‑1 (MQ=255)
                                                                                                                                          cgcgGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCtgt           <  2:1017794/150‑1 (MQ=255)
                                                                                                                                            cgGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTgg         <  2:483401/150‑1 (MQ=255)
                                                                                                                                            cgGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTgg         <  1:907461/150‑1 (MQ=255)
                                                                                                                                            cgGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCAAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTgg         <  2:511838/150‑1 (MQ=255)
                                                                                                                                              gCAATCCATTGATTGCCGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGa                                               >  2:895232/1‑110 (MQ=255)
                                                                                                                                              gCAATCCATTGATTGCCGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGa                                               <  1:895232/110‑1 (MQ=255)
                                                                                                                                              gCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGtt       >  1:1119994/1‑150 (MQ=255)
                                                                                                                                                 atcaatGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAg    >  2:889447/1‑150 (MQ=255)
                                                                                                                                                   gaatAGAGTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGtt  <  2:1411714/149‑1 (MQ=255)
                                                                                                                                                   caatGGATAGCAGGAGACAGGACAGGAGTGCCGTGATCGACGACCACCTCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCGGGGTGCTGAACACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGtt  <  2:416383/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGTT  >  NC_002947/109811‑110107

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