Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,684,277 G→A 9.2% 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,2770GA 9.2% 180.5 / 7.0 66E286K (GAG→AAG) PP_1478xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/4);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

AGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGC  >  NC_002947/1684134‑1684421
                                                                                                                                               |                                                                                                                                                
aGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGAGCTTTGCAGAgcgc                                                                                                                                             >  8:294007/1‑149 (MQ=255)
    cGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTgc                                                                                                                                         <  2:297147/149‑1 (MQ=255)
    cGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtct                                                                                                                                         <  1:227735/149‑3 (MQ=255)
        tGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtctgata                                                                                                                                     <  2:271716/149‑7 (MQ=255)
          gAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCaa                                                                                                                                     <  3:140730/147‑1 (MQ=255)
          gAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTCTGCAAAGCGCGtctgata                                                                                                                                     <  5:27797/147‑7 (MQ=255)
           aaCCCGGATGAAACCGTCGAACTGGATCGTCGCCTGCGGGTGATGGGTGTCGACCTGATCGATGTGTCATCCGGTGGTACCTCGGTCAAAGCCGAAATACCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGATCGCGTGAGCAAgg                                                                                                                                   >  6:21286/1‑148 (MQ=255)
               cGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCa                                                                                                                              <  5:210434/149‑1 (MQ=255)
                ggATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTTTCCTCCGGGGGTACCTCGGTCAACGCCGAAATTCCTACAGTCCCTGTCTATCAGACGCTCTATGCAGAGCGCGTGCGCAAGGAAGCag                                                                                                                             <  5:21286/149‑1 (MQ=255)
                   tGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAgaga                                                                                                                          <  5:185624/149‑1 (MQ=255)
                             gAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCgg                                                                                                                <  4:285369/149‑1 (MQ=255)
                                 tGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCAcc                                                                                                            >  6:88197/1‑149 (MQ=255)
                                   gCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGt                                                                                                          >  2:254492/1‑149 (MQ=255)
                                       gtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCgg                                                                                                       >  4:99911/1‑148 (MQ=255)
                                       gtcgCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCg                                                                                                        <  5:196952/147‑1 (MQ=255)
                                               cGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCa                                                                                               <  8:316562/148‑1 (MQ=255)
                                                 ggTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGTGATCGCCACCGGCACCGTCGGCATGATCACTg                                                                                            <  1:42205/149‑1 (MQ=255)
                                                         gtgtCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGcc                                                                                    >  2:239766/1‑149 (MQ=255)
                                                         gtgtCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGcc                                                                                    >  4:326970/1‑149 (MQ=255)
                                                           gtCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGcc                                                                                    >  2:235041/1‑147 (MQ=255)
                                                            tCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAg                                                                                 >  5:107980/1‑149 (MQ=255)
                                                                    ctCGATGTGTCCTCCGGTGGTACCTCGGTCTACGCCGATATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGGCACCGGCACCGTCGGCATGATCACTGAAGCTGCCCAGGCCGAGc                                                                          <  6:208081/147‑1 (MQ=255)
                                                                    aTCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGca                                                                         >  5:231895/1‑149 (MQ=255)
                                                                    aTCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGc                                                                          <  3:326970/148‑1 (MQ=255)
                                                                     tCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGca                                                                         >  3:58243/1‑148 (MQ=255)
                                                                       gATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACAt                                                                      >  3:108301/1‑149 (MQ=255)
                                                                        aTGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAAGCCGAGCACAt                                                                      >  2:318157/1‑148 (MQ=255)
                                                                           tgtCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACt                                                                                             <  4:221113/122‑1 (MQ=255)
                                                                           tgtCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCt                                                                   <  1:254492/148‑1 (MQ=255)
                                                                             tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                 <  3:99911/148‑1 (MQ=255)
                                                                             tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                 <  5:88197/148‑1 (MQ=255)
                                                                             tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGc                                                                 <  1:239766/148‑1 (MQ=255)
                                                                             tcctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCg                                                                >  2:28077/1‑149 (MQ=255)
                                                                               ctccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTa                                                              >  1:169272/1‑149 (MQ=255)
                                                                                 ccGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGACATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTAcc                                                            >  6:201712/1‑149 (MQ=255)
                                                                                      ggtACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGCAGCAGAGATCGCCACCGGCACCTTCGGCATGATCACTGAACCTTCCCAGGCATCGCAAATTCTGAGTACCggg                                                         >  4:202812/1‑147 (MQ=255)
                                                                                      ggtACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTAGCAGACGCGCTTTGCAGCGCGCGTGCGCAAGGACGCAGAGATCGCCACCGGCACCGTCGGCATGAACCCTGCACCTGCCCCGGCCGCGCCCCTTCTGCGTACCggg                                                         >  4:207079/1‑147 (MQ=255)
                                                                                         aCCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGc                                                    >  7:308637/1‑149 (MQ=255)
                                                                                          ccTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAgg                                                     >  7:206576/1‑147 (MQ=255)
                                                                                          ccTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAgg                                                     <  8:206576/147‑1 (MQ=255)
                                                                                          ccTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGAGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAggac                                                   >  5:214283/1‑147 (MQ=255)
                                                                                              ggTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGAc                                                <  1:28077/148‑1 (MQ=255)
                                                                                              ggTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                               >  5:166708/1‑149 (MQ=255)
                                                                                              ggTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                               <  1:246073/149‑1 (MQ=255)
                                                                                              ggTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAAGCCGAGCACATTCTGCGTACCGGGCAGGCCGa                                                 <  1:318157/147‑1 (MQ=255)
                                                                                              ggTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGACATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACg                                               <  5:201712/149‑1 (MQ=255)
                                                                                               gTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGt                                              <  6:231895/149‑1 (MQ=255)
                                                                                               gTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGt                                              >  1:320993/1‑149 (MQ=255)
                                                                                                tCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGt                                              >  4:13169/1‑148 (MQ=255)
                                                                                                   aCGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATc                                          <  6:166708/149‑1 (MQ=255)
                                                                                                                tACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCgg                              >  1:159990/1‑148 (MQ=255)
                                                                                                                 aCAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTcc                                      >  5:240128/1‑139 (MQ=255)
                                                                                                                 aCAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTcc                                      <  6:240128/139‑1 (MQ=255)
                                                                                                                  cAGGCCCTGGCTCTCCGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGCGCTCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCCGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                            <  3:207079/148‑1 (MQ=255)
                                                                                                                  cAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGa                            <  4:108301/148‑1 (MQ=255)
                                                                                                                  cAGGCCCTGGCTATCAGACGCGCTGTGCATAGCGCGTGCGCAATGAAGCAGAGATCGCCACAGGCACCGTCGGCATTATCACTGAACCTGCCCAGGCCGATCACATTCTGCGTACCGGGCAGTCCGACGTGATCTTCCTGGCTCGGGa                            <  3:202812/148‑1 (MQ=255)
                                                                                                                     gCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAgctg                        <  8:308637/149‑1 (MQ=255)
                                                                                                                          ggCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCAGCCACGGCAGCAGCTgcg                   >  4:1845/1‑149 (MQ=255)
                                                                                                                           gCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCtgctgc                    <  6:107980/147‑1 (MQ=255)
                                                                                                                               tCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATc              <  2:320993/149‑1 (MQ=255)
                                                                                                                                 gcaCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATc                                          >  3:35342/3‑119 (MQ=255)
                                                                                                                                 gcaCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATc                                          <  4:35342/117‑1 (MQ=255)
                                                                                                                                   aCGCGCTCTGCAAAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTa          >  5:137686/1‑149 (MQ=255)
                                                                                                                                     gcgcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACt        >  3:119558/1‑149 (MQ=255)
                                                                                                                                       gcTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTgg      <  3:13169/149‑1 (MQ=255)
                                                                                                                                            gCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGc  <  2:169272/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AGGACGGCTGGAACCCGGATGAAACCGTCGAACTGGCTCGTCGCCTGCGGGTGCTGGGTGTCGACCTGATCGATGTGTCCTCCGGTGGTACCTCGGTCAACGCCGAAATTCCTACAGGCCCTGGCTATCAGACGCGCTTTGCAGAGCGCGTGCGCAAGGAAGCAGAGATCGCCACCGGCACCGTCGGCATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGC  >  NC_002947/1684134‑1684421

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