Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 444,596 T→A 5.1% intergenic (+75/‑99) bioD → / → PP_0367 dethiobiotin synthetase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947444,5960TA 5.1% 541.9 / 12.2 176intergenic (+75/‑99)bioD/PP_0367dethiobiotin synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (85/80);  new base A (4/5);  total (90/86)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 181 total aligned reads displayed.

GGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGACATGGAAATCACC  >  NC_002947/444447‑444745
                                                                                                                                                     |                                                                                                                                                     
gggCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGttt                                                                                                                                                       >  1:1219159/1‑150 (MQ=255)
    gggTGCCCAAGCGCAAGCAGGCCAGGGCCGACATGGGGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTAccc                                                                                                                                                   <  1:614851/150‑1 (MQ=255)
      gTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCcgc                                                                                                                                                 <  2:1308963/150‑1 (MQ=255)
          ccAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaag                                                                                                                                             >  1:2134865/1‑150 (MQ=255)
          ccAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaag                                                                                                                                             >  2:228227/1‑150 (MQ=255)
            aaGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaagaa                                                                                                                                           <  1:32147/150‑1 (MQ=255)
                tCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCa                                                                                                                                       >  2:347439/1‑150 (MQ=255)
                  aaGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGt                                                                                                                                     >  2:1975628/1‑150 (MQ=255)
                    gCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGc                                                                                                                                   >  2:1197269/1‑150 (MQ=255)
                      aGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAg                                                                                                                                 <  1:287956/150‑1 (MQ=255)
                       ggCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAgg                                                                                                                                <  2:1985830/150‑1 (MQ=255)
                        gCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCcgc                                                                                                                                                 <  2:1947039/132‑1 (MQ=255)
                        gCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCcgc                                                                                                                                                 >  1:1947039/1‑132 (MQ=255)
                         ccAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGaa                                                                                                                              >  2:2192728/1‑150 (MQ=255)
                         ccAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACACGCGAAGATGCCAGTGCAGGaa                                                                                                                              <  2:643575/150‑1 (MQ=255)
                               ccGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCa                                                                                                                        >  2:1619869/1‑150 (MQ=255)
                                cGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGc                                                                                                                                   <  2:1206697/138‑1 (MQ=255)
                                  aCATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAgg                                                                                                                     >  2:2004980/1‑150 (MQ=255)
                                  aCATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAgg                                                                                                                     >  1:822801/1‑150 (MQ=255)
                                              aGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTgg                                                                                                         <  1:518874/150‑1 (MQ=255)
                                                 aTCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCAt                                                                                                      <  2:1286377/150‑1 (MQ=255)
                                                  tCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATAcc                                                                                                                          <  1:1410061/129‑1 (MQ=255)
                                                    tGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGa                                                                                                   >  2:1746253/1‑150 (MQ=255)
                                                     gCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaaga                                                                                                                                            <  2:2325021/108‑1 (MQ=255)
                                                     gCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaaga                                                                                                                                            >  1:2325021/1‑108 (MQ=255)
                                                     gCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAg                                                                                                  <  1:1803496/150‑1 (MQ=255)
                                                         cTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCAt                                                                                              >  1:2202778/1‑150 (MQ=255)
                                                         cTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCAt                                                                                              >  2:458899/1‑150 (MQ=255)
                                                           ggATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTa                                                                                            <  1:1975628/150‑1 (MQ=255)
                                                            gATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAg                                                                                           >  2:769087/1‑150 (MQ=255)
                                                             aTCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAgg                                                                                          <  2:455159/150‑1 (MQ=255)
                                                                tgctgGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCga                                                                                                                                               <  2:278851/94‑1 (MQ=255)
                                                                tgctgGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGAt                                                                                       <  1:2004980/150‑1 (MQ=255)
                                                                   tgGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaag                                                                                                                                             >  2:2406316/1‑93 (MQ=255)
                                                                       ttAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATAc                                                                                                                           <  1:1508102/107‑1 (MQ=255)
                                                                       ttAGTCTGGGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACg                                                                                <  2:541936/150‑1 (MQ=255)
                                                                         aGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACggg                                                                              <  2:884316/150‑1 (MQ=255)
                                                                          gTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaaga                                                                                                                                            <  1:1351089/87‑1 (MQ=255)
                                                                          gTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCgaaga                                                                                                                                            >  2:1351089/1‑87 (MQ=255)
                                                                           tCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGa                                                                                                   <  2:2257235/127‑1 (MQ=255)
                                                                           tCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGcc                                                                            >  2:399306/1‑150 (MQ=255)
                                                                           tCTGTTGCAGGCCCGGCCTCTTCCCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGAGAAGAAGCCCCGGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGcc                                                                            >  1:70457/1‑150 (MQ=39)
                                                                             tGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACg                                                                                <  1:581041/144‑1 (MQ=255)
                                                                             tGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACg                                                                                >  2:581041/1‑144 (MQ=255)
                                                                             tGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCtt                                                                          >  2:220177/1‑150 (MQ=255)
                                                                             tGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCtt                                                                          >  2:1459447/1‑150 (MQ=255)
                                                                             tGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCtt                                                                          >  1:1530839/1‑150 (MQ=255)
                                                                              gTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAg                                                                                                                                 <  1:1935056/94‑1 (MQ=255)
                                                                                     gCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGcccc                                                                  >  2:150188/1‑150 (MQ=255)
                                                                                       ccGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGGGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTa                                                                >  1:2407246/1‑150 (MQ=255)
                                                                                        cGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAg                                                               >  1:2178946/1‑150 (MQ=255)
                                                                                        cGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAg                                                               >  2:2302176/1‑150 (MQ=255)
                                                                                        cGGCCTCTTCGCGGGTAAAACCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAg                                                               >  2:2301718/1‑150 (MQ=255)
                                                                                         ggCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGc                                                              >  2:1018050/1‑150 (MQ=255)
                                                                                         ggCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGc                                                              <  1:482274/150‑1 (MQ=255)
                                                                                           ccTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCAt                                                                                              <  1:802919/116‑1 (MQ=255)
                                                                                           ccTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCAt                                                                                              >  2:802919/1‑116 (MQ=255)
                                                                                               ttCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGttt                                                        <  2:1136005/150‑1 (MQ=255)
                                                                                                  gcgGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAggg                                                                                         >  1:328591/1‑114 (MQ=255)
                                                                                                    gggTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAAtt                                                   <  1:399306/150‑1 (MQ=255)
                                                                                                     ggTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTg                                                  >  1:2017123/1‑150 (MQ=255)
                                                                                                      gTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCAAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTgg                                                 <  1:2111362/150‑1 (MQ=255)
                                                                                                      gTAAACCCGCTCCCACGAGGGAACGCATACCGTTGTGGGCGACGGCCAACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTgg                                                 <  1:1467352/150‑1 (MQ=255)
                                                                                                       tAAACCCGCTCCTACGAGGGAACGCATCCCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                                                <  1:536105/150‑1 (MQ=255)
                                                                                                       tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                                                <  1:125434/150‑1 (MQ=255)
                                                                                                       tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                                                <  1:228227/150‑1 (MQ=255)
                                                                                                       tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                                                <  1:458899/150‑1 (MQ=255)
                                                                                                       tAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGACATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTggg                                                <  2:2217014/150‑1 (MQ=255)
                                                                                                           cccGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTg                                            >  2:547193/1‑150 (MQ=255)
                                                                                                            ccGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCGATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAGTTGGGGCTGt                                           >  1:365474/1‑150 (MQ=255)
                                                                                                                tCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGt                                       >  1:894468/1‑150 (MQ=255)
                                                                                                                        gggAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGACCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGc                               <  1:2225869/150‑1 (MQ=255)
                                                                                                                          gAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGt                             >  1:1437715/1‑150 (MQ=255)
                                                                                                                          gAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATACAGCGt                             <  2:489612/150‑1 (MQ=255)
                                                                                                                             ccaaTACCGTCGTGGCCGCGGGTTCACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCaa                          <  2:2029405/147‑1 (MQ=255)
                                                                                                                              gCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAAt                         >  2:599498/1‑150 (MQ=255)
                                                                                                                              gCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAAt                         <  1:769087/150‑1 (MQ=255)
                                                                                                                               cATACCGTTGTGGGCGCGGGTTTACCCGCGAATAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATg                        >  2:1905175/1‑150 (MQ=255)
                                                                                                                                 tACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGTCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGa                      <  1:2302176/150‑1 (MQ=255)
                                                                                                                                 tACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGGCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGa                      <  1:2301718/150‑1 (MQ=255)
                                                                                                                                    cGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACgg                                                                               <  2:1604443/90‑1 (MQ=255)
                                                                                                                                     gTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTg                                                  <  2:465752/118‑1 (MQ=255)
                                                                                                                                      ttGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCt                 >  2:1374537/1‑150 (MQ=255)
                                                                                                                                      ttGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCt                 >  1:361314/1‑150 (MQ=255)
                                                                                                                                      ttGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCt                 >  1:1614224/1‑150 (MQ=255)
                                                                                                                                      ttGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCt                 <  1:1018050/150‑1 (MQ=255)
                                                                                                                                      tcttcttcgcgGGTAAACCCGCGCAGAAGCTAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGCCTTTAAACGGGCCTTTTCACGCCTACGCTGGTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCt                 <  1:420653/143‑1 (MQ=255)
                                                                                                                                        ttcttcgcgGGTAAACCCGCGCAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGa               >  2:420653/6‑150 (MQ=255)
                                                                                                                                        gtgGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGa               <  2:361314/150‑1 (MQ=255)
                                                                                                                                         tcttcgcgGGTAAACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGa                      <  2:1434029/138‑1 (MQ=255)
                                                                                                                                            gcgcGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCttt                                                                         <  2:1502622/88‑1 (MQ=255)
                                                                                                                                            gcgcGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCttt                                                                         >  1:1502472/1‑88 (MQ=255)
                                                                                                                                             cgcgtatttATCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGACATgg          <  1:1707669/144‑1 (MQ=255)
                                                                                                                                             cgcgGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGACATgg          <  1:1197269/150‑1 (MQ=255)
                                                                                                                                             cgcgGGTAAACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGACATgg          <  2:2115282/150‑1 (MQ=255)
                                                                                                                                              gcgGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCaa                                                                                                                       >  1:2355974/1‑40 (MQ=255)
                                                                                                                                              gcgGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCaa                                                                                                                       <  2:2355974/40‑1 (MQ=255)
                                                                                                                                              gcgGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCt                 >  1:1049328/1‑142 (MQ=255)
                                                                                                                                                 ggTTTACCCGCGACGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGACATGGAAAt      >  2:1244196/1‑150 (MQ=255)
                                                                                                                                                     tACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGACATGGAAATCAcc  <  1:547193/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGCGGGTTCCCAAGCTCAAGCAGGCCAGTGCCGACATGGTGGCTGAGCATCTGCAACTGGATCTGCTGGATTAGTCTGTTGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGAACGCATACCGTTGTGGGCGCGGGTTTACCCGCGAAGAAGCCAGTGCAGGAATACCCAAGGCCTATCAACTGGCATTGAGCCATTAGGGATTTAAACGGGCCTTTTCGCCCCTAGCCTGTTTTAATTGGGGCTGTTCGTCATCCAGCGTCAATGGAGTCCTGACATGGAAATCACC  >  NC_002947/444447‑444745

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