Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,684,314 T→C 7.0% F284S (TTT→TCT)  PP_1478 → xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,684,3140TC 7.0% 257.2 / 6.5 85F284S (TTT→TCT) PP_1478xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base T (42/37);  new base C (2/4);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CTGGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCG  >  NC_002947/1684168‑1684462
                                                                                                                                                  |                                                                                                                                                    
ctgGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACTCGCTTTGc                                                                                                                                                   <  1:1810045/150‑1 (MQ=255)
 tgGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCa                                                                                                                                                  <  1:1240200/150‑1 (MQ=255)
    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                                                                                                                                    >  1:1398714/1‑103 (MQ=255)
                                                              tGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                    <  2:1398714/103‑1 (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                                                                                                  >  2:2086501/1‑128 (MQ=255)
                                                                       aCCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATg                                                                                                  <  1:2086501/128‑1 (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                                                                      >  2:2094078/1‑140 (MQ=255)
                                                                                       ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGcaca                                                                      <  1:2094078/140‑1 (MQ=255)
                                                                                       ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAc                                                            >  1:85369/1‑150 (MQ=255)
                                                                                        tccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCAc                                                                                             <  2:94202/116‑1 (MQ=255)
                                                                                        tccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCAc                                                                                             >  1:94202/1‑116 (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                                                                                                           >  2:2234465/1‑81 (MQ=255)
                                                                                                             gCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCAcc                                                                                                           <  1:2234465/81‑1 (MQ=255)
                                                                                                             gCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTc                                      >  1:1616414/1‑150 (MQ=255)
                                                                                                               cGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCt                                    >  2:1331668/1‑150 (MQ=255)
                                                                                                                gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTg                                   >  2:1685425/1‑150 (MQ=255)
                                                                                                                gAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTg                                   <  1:800423/150‑1 (MQ=255)
                                                                                                                   aTTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCt                                <  1:1570217/150‑1 (MQ=255)
                                                                                                                        tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                           <  2:929552/150‑1 (MQ=255)
                                                                                                                        tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                           <  2:4660/150‑1 (MQ=255)
                                                                                                                         aCAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                         >  1:532348/1‑38 (MQ=18)
                                                                                                                         aCAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtc                                                                                                                                         <  2:532348/39‑2 (MQ=18)
                                                                                                                             gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCt                                    >  1:611596/1‑136 (MQ=255)
                                                                                                                             gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCt                                    <  2:611596/136‑1 (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                   >  2:628901/1‑138 (MQ=255)
                                                                                                                                            cgcgCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                   <  1:628901/138‑1 (MQ=255)
                                                                                                                                              cgcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTgg     <  1:2325312/150‑1 (MQ=255)
                                                                                                                                               gcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGc    >  1:121832/1‑150 (MQ=255)
                                                                                                                                               gcTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGc    >  2:2384559/1‑150 (MQ=255)
                                                                                                                                                 tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGc                                 <  1:2238456/119‑1 (MQ=255)
                                                                                                                                                 tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGc                                 >  2:2238456/1‑119 (MQ=255)
                                                                                                                                                 tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                   >  2:1479494/1‑133 (MQ=255)
                                                                                                                                                 tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                   <  1:1479494/133‑1 (MQ=255)
                                                                                                                                                 tttGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTCATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCg  >  2:722417/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGGGTCGAGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCG  >  NC_002947/1684168‑1684462

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: ‑