Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,684,319 G→A 5.7% E286K (GAG→AAG)  PP_1478 → xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,684,3190GA 5.7% 206.3 / 2.0 70E286K (GAG→AAG) PP_1478xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAT  >  NC_002947/1684171‑1684468
                                                                                                                                                    |                                                                                                                                                     
ggTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCaga                                                                                                                                                      >  2:369642/1‑150 (MQ=255)
ggTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCaga                                                                                                                                                      <  2:876197/150‑1 (MQ=255)
      ggACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGt                                                                                                                                                <  2:1153062/150‑1 (MQ=255)
       gACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTg                                                                                                                                               >  1:1187680/1‑150 (MQ=255)
          ggCTGGATCCCGGATGAAACCGTCGAACGGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCACGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtctga                                                                                                                                            <  2:754294/150‑5 (MQ=255)
             tGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCCACCTGATCGATGTGTCCTCCTGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATAAGAAACGCTTTGCAGAGCACGTGCGCaac                                                                                                                                         >  2:374418/1‑149 (MQ=255)
              ggAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtct                                                                                                                                              <  2:1844862/144‑3 (MQ=255)
              ggAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtct                                                                                                                                              >  1:1844862/1‑142 (MQ=255)
                   ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCa                                                                                                                                   <  1:483372/150‑1 (MQ=255)
                           aaCCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGc                                                                                                                           >  1:474583/1‑150 (MQ=255)
                              cGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCAc                                                                                                                        <  1:369642/150‑1 (MQ=255)
                               gTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCaga                                                                                                                                                      >  2:685752/1‑119 (MQ=255)
                               gTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCaga                                                                                                                                                      <  1:685752/119‑1 (MQ=255)
                               gTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCAcc                                                                                                                       >  1:1252242/1‑150 (MQ=255)
                                        gCTCGTCGCCCGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTc                                                                                                              <  1:1753500/150‑1 (MQ=255)
                                          tcgtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCgg                                                                                                            >  1:880482/1‑150 (MQ=255)
                                            gtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCa                                                                                                          >  1:1730095/1‑150 (MQ=255)
                                            gtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCa                                                                                                          >  1:232852/1‑150 (MQ=255)
                                                 cTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGa                                                                                                       <  2:465318/148‑1 (MQ=255)
                                                 cTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGa                                                                                                       >  1:465318/1‑148 (MQ=255)
                                                  tGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCa                                                                                                    <  1:741698/150‑1 (MQ=255)
                                                               tgtCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGcccc                                                                                       >  1:583735/1‑149 (MQ=255)
                                                               tgtCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCa                                                                                       >  1:1016011/1‑150 (MQ=255)
                                                                  cGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGc                                                                                    <  2:1252242/150‑1 (MQ=255)
                                                                   gACCTGATCGATGTGTCCTCCGGTGGTACCCCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGcc                                                                                   >  1:1508403/1‑150 (MQ=255)
                                                                    aCCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGaa                                                                                               <  2:292936/137‑1 (MQ=255)
                                                                    aCCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGaa                                                                                               >  1:292936/1‑137 (MQ=255)
                                                                     ccTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGa                                                                                 <  1:374418/150‑1 (MQ=255)
                                                                       tGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGCAATTCATCCAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGc                                                                               <  2:1187680/150‑1 (MQ=255)
                                                                        gATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAGATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGca                                                                              >  1:1340498/1‑150 (MQ=255)
                                                                         aTCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGtc                                                                                                                                               >  2:595292/1‑83 (MQ=255)
                                                                         aTCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGtc                                                                                                                                               <  1:595292/84‑2 (MQ=255)
                                                                         aTCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGcac                                                                             >  1:775282/1‑150 (MQ=255)
                                                                              tgtgTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAATCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCt                                                                        <  1:693119/150‑1 (MQ=255)
                                                                                tgtCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                      <  2:880482/150‑1 (MQ=255)
                                                                                    ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                                  >  1:1261578/1‑150 (MQ=255)
                                                                                                 tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGaa                                                                                               >  2:1484174/1‑108 (MQ=255)
                                                                                                 tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGaa                                                                                               <  1:1484174/108‑1 (MQ=255)
                                                                                                 tCGGTCAACGCCGAAATGCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCAATGAACCTGTCCAGAACGATCACATTCTGCGTACCGGGCAGGCCgcc                                                     >  2:990220/1‑148 (MQ=255)
                                                                                                     tCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGa                                                 >  1:478187/1‑150 (MQ=255)
                                                                                                     tCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGa                                                 <  2:775282/150‑1 (MQ=255)
                                                                                                     tCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGa                                                 >  1:1608544/1‑150 (MQ=255)
                                                                                                          gCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTc                                            >  1:1440226/1‑150 (MQ=255)
                                                                                                            cGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGCGATCTTCCt                                          <  2:474583/150‑1 (MQ=255)
                                                                                                             gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATc                                                                                                                             <  1:363392/66‑1 (MQ=255)
                                                                                                             gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATc                                                                                                                             >  2:363392/1‑66 (MQ=255)
                                                                                                                aTTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCa                                                            <  2:458394/128‑1 (MQ=255)
                                                                                                                aTTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCa                                                            >  1:458394/1‑128 (MQ=255)
                                                                                                                 ttCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTc                                     >  2:967312/1‑150 (MQ=255)
                                                                                                                   ccTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCgg                                   <  1:967312/150‑1 (MQ=255)
                                                                                                                   ccTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCgg                                   <  2:868322/150‑1 (MQ=255)
                                                                                                                     tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                                 >  1:1897557/1‑150 (MQ=255)
                                                                                                                     tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                                 <  2:1730095/150‑1 (MQ=255)
                                                                                                                      aCAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAg                                >  2:1018463/1‑150 (MQ=255)
                                                                                                                          gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCg                                                                                                             <  1:162175/69‑1 (MQ=255)
                                                                                                                          gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCg                                                                                                             >  2:162175/1‑69 (MQ=255)
                                                                                                                          gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                      <  2:172056/108‑1 (MQ=255)
                                                                                                                          gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                      >  1:172056/1‑108 (MQ=255)
                                                                                                                           cccTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAgctgct                           >  1:923097/1‑150 (MQ=255)
                                                                                                                             cTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                         >  2:1831193/1‑150 (MQ=255)
                                                                                                                             cTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                         <  1:1018463/150‑1 (MQ=255)
                                                                                                                               ggCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGCACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCATCCTGGCTCGGGAGCTGCTgcgc                       >  1:1761090/1‑150 (MQ=255)
                                                                                                                               ggCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTgcgc                       <  2:1311464/150‑1 (MQ=255)
                                                                                                                                gCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGcgcg                      >  2:1393291/1‑150 (MQ=255)
                                                                                                                                 cTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGa                     <  1:1940664/150‑1 (MQ=255)
                                                                                                                                  ttgcgcaCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTTGCTCGGTAGCTGCTCCGCGAt                    >  2:1859409/7‑150 (MQ=255)
                                                                                                                                     cAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCt                 >  1:430146/1‑150 (MQ=255)
                                                                                                                                      aGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCtt                <  2:1897557/150‑1 (MQ=255)
                                                                                                                                          gcgcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCt                 >  1:1178535/1‑145 (MQ=255)
                                                                                                                                          gcgcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCt                 <  2:1178535/145‑1 (MQ=255)
                                                                                                                                           cgcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTgg           <  2:1440226/150‑1 (MQ=255)
                                                                                                                                               ttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGc                                                         <  2:450565/100‑1 (MQ=255)
                                                                                                                                               ttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGc                                                         >  1:450565/1‑100 (MQ=255)
                                                                                                                                                 gCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTg     >  1:890939/1‑150 (MQ=255)
                                                                                                                                                    gagCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTa               >  2:198095/1‑137 (MQ=255)
                                                                                                                                                    gagCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTa               <  1:198095/137‑1 (MQ=255)
                                                                                                                                                    gagCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAt  <  2:430146/150‑1 (MQ=255)
                                                                                                                                                    gagCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAt  >  1:153338/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAT  >  NC_002947/1684171‑1684468

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