Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,684,319 G→A 6.9% 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 6.9% 258.9 / 6.7 87E286K (GAG→AAG) PP_1478xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (44/37);  new base A (2/4);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAT  >  NC_002947/1684172‑1684468
                                                                                                                                                   |                                                                                                                                                     
gTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCagag                                                                                                                                                     >  1:1756903/1‑150 (MQ=255)
      gACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                               <  2:117784/150‑2 (MQ=255)
      gACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                               <  2:223463/150‑2 (MQ=255)
      gACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGCGTGTCGACCTTATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                               <  1:286514/150‑2 (MQ=255)
         ggCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTTCAGAGCGCGTgcgc                                                                                                                                            >  1:929552/1‑150 (MQ=255)
             ggAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAgg                                                                                                                                        <  1:582496/150‑1 (MQ=255)
                  ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATGGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCa                                                                                                                                   >  2:800423/1‑150 (MQ=255)
                  ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCa                                                                                                                                   <  2:2058541/150‑1 (MQ=255)
                  ccGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCa                                                                                                                                   <  1:2151491/150‑1 (MQ=255)
                        gAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATc                                                                                                                             <  2:753796/150‑1 (MQ=255)
                            ccGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCa                                                                                                                         >  2:1570217/1‑150 (MQ=255)
                             cGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCAc                                                                                                                        >  2:1981859/1‑150 (MQ=255)
                              gTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCAcc                                                                                                                       >  1:984046/1‑150 (MQ=255)
                                         tcgtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCgg                                                                                                            >  2:296504/1‑150 (MQ=255)
                                         tcgtcgCCTGCGGGTGCTGGGTGTCGACCGGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACcggcgg                                                                                                            >  1:1109564/1‑150 (MQ=255)
                                            tcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCAt                                                                                                         <  2:1756903/150‑1 (MQ=255)
                                              gCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGa                                                                                                       >  1:1213525/1‑150 (MQ=255)
                                                      gTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGaa                                                                                               <  1:662453/150‑1 (MQ=255)
                                                          tGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                          <  2:1398714/103‑1 (MQ=255)
                                                          tGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                          >  1:1398714/1‑103 (MQ=255)
                                                                tCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAgg                                                                                     <  1:752304/150‑1 (MQ=255)
                                                                tCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAgg                                                                                     >  2:2277297/1‑150 (MQ=255)
                                                                  gACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGAt                                                                                                      <  1:553956/131‑1 (MQ=255)
                                                                  gACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGAt                                                                                                      >  2:553956/1‑131 (MQ=255)
                                                                  gACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGcc                                                                                   >  2:785372/1‑150 (MQ=255)
                                                                  gACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGcc                                                                                   >  1:366227/1‑150 (MQ=255)
                                                                   aCCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATg                                                                                                        <  1:2086501/128‑1 (MQ=255)
                                                                   aCCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATg                                                                                                        >  2:2086501/1‑128 (MQ=255)
                                                                     cTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAg                                                                                <  1:1242191/150‑1 (MQ=255)
                                                                        aTCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGcac                                                                             >  2:553902/1‑150 (MQ=255)
                                                                        aTCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGcac                                                                             >  1:70210/1‑150 (MQ=255)
                                                                            aTGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTc                                                                         >  1:2304231/1‑150 (MQ=255)
                                                                               tgtCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                      >  2:1947326/1‑150 (MQ=255)
                                                                               tgtCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                      <  2:2304231/150‑1 (MQ=255)
                                                                                 tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                          >  2:2130954/1‑80 (MQ=255)
                                                                                 tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                          <  1:2130954/80‑1 (MQ=255)
                                                                                 tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATg                                                                                                        >  2:789624/1‑114 (MQ=255)
                                                                                 tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATg                                                                                                        <  1:789624/114‑1 (MQ=255)
                                                                                 tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGt                                                                    >  2:2325312/1‑150 (MQ=255)
                                                                                   ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGcaca                                                                            <  1:2094078/140‑1 (MQ=255)
                                                                                   ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGcaca                                                                            >  2:2094078/1‑140 (MQ=255)
                                                                                   ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                                  >  1:85369/1‑150 (MQ=255)
                                                                                    tccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCAc                                                                                                   >  1:94202/1‑116 (MQ=255)
                                                                                    tccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCAc                                                                                                   <  2:94202/116‑1 (MQ=255)
                                                                                     ccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCACGGAAGCAGAGATCGGGACTTGGATCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCg                                                                <  1:296504/150‑1 (MQ=255)
                                                                                          ggtACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAg                                                           <  2:366227/150‑1 (MQ=255)
                                                                                           gtACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGc                                                                                    <  2:493535/124‑1 (MQ=255)
                                                                                           gtACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGc                                                                                    >  1:493535/1‑124 (MQ=255)
                                                                                               cTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGa                                                      >  2:2015733/1‑150 (MQ=255)
                                                                                                tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTgc                                                                                                                                              <  1:912291/61‑1 (MQ=255)
                                                                                                tCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTgc                                                                                                                                              >  2:912291/1‑61 (MQ=255)
                                                                                                tCGGTCAACGCCGAAATTCCTACAGCCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGAc                                                     <  1:1134619/150‑1 (MQ=255)
                                                                                                   gTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAgg                                                          <  1:1492318/142‑1 (MQ=255)
                                                                                                   gTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAgg                                                          >  2:1492318/1‑142 (MQ=255)
                                                                                                    tCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGa                                                 <  1:356555/150‑1 (MQ=255)
                                                                                                         gCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCAcc                                                                                                                 <  1:2234465/81‑1 (MQ=255)
                                                                                                         gCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCAcc                                                                                                                 >  2:2234465/1‑81 (MQ=255)
                                                                                                         gCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTc                                            >  1:1616414/1‑150 (MQ=255)
                                                                                                           cGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCt                                          >  2:1331668/1‑150 (MQ=255)
                                                                                                            gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTg                                         <  1:800423/150‑1 (MQ=255)
                                                                                                            gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTg                                         >  2:1685425/1‑150 (MQ=255)
                                                                                                               aTTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCt                                      <  1:1570217/150‑1 (MQ=255)
                                                                                                                    tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                                 <  2:929552/150‑1 (MQ=255)
                                                                                                                    tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                                 <  2:4660/150‑1 (MQ=255)
                                                                                                                     aCAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                               <  2:532348/39‑2 (MQ=18)
                                                                                                                     aCAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                               >  1:532348/1‑38 (MQ=18)
                                                                                                                         gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCt                                          <  2:611596/136‑1 (MQ=255)
                                                                                                                         gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCt                                          >  1:611596/1‑136 (MQ=255)
                                                                                                                           ccTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACGGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGctgctg                          >  2:2195762/1‑150 (MQ=255)
                                                                                                                           ccTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGctgctg                          >  1:1858467/1‑150 (MQ=255)
                                                                                                                              ggCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACt                                                                                                  <  1:2137198/75‑1 (MQ=255)
                                                                                                                              ggCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACt                                                                                                  >  2:2137198/1‑75 (MQ=255)
                                                                                                                              ggCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTgcgc                       >  2:217699/1‑150 (MQ=255)
                                                                                                                                cTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGa                     <  2:1616414/150‑1 (MQ=255)
                                                                                                                                 tATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGAt                    >  1:1560981/1‑150 (MQ=255)
                                                                                                                                  aTCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATc                   <  1:2015733/150‑1 (MQ=255)
                                                                                                                                   tCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGccc                                                                                        <  2:419614/80‑1 (MQ=255)
                                                                                                                                   tCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGccc                                                                                        >  1:419614/1‑80 (MQ=255)
                                                                                                                                        cgcgCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                         <  1:628901/138‑1 (MQ=255)
                                                                                                                                        cgcgCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                         >  2:628901/1‑138 (MQ=255)
                                                                                                                                          cgcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTgg           <  1:2325312/150‑1 (MQ=255)
                                                                                                                                           gcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGc          >  1:121832/1‑150 (MQ=255)
                                                                                                                                           gcTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGc          >  2:2384559/1‑150 (MQ=255)
                                                                                                                                             tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGc                                       >  2:2238456/1‑119 (MQ=255)
                                                                                                                                             tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGc                                       <  1:2238456/119‑1 (MQ=255)
                                                                                                                                             tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                         >  2:1479494/1‑133 (MQ=255)
                                                                                                                                             tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                         <  1:1479494/133‑1 (MQ=255)
                                                                                                                                             tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTCATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCg        >  2:722417/1‑150 (MQ=255)
                                                                                                                                                   gagCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAgctgct                           >  2:608171/1‑125 (MQ=255)
                                                                                                                                                   gagCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAgctgct                           <  1:608171/125‑1 (MQ=255)
                                                                                                                                                   gagCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAt  >  1:1282215/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCAT  >  NC_002947/1684172‑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: ‑