Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,976,352 G→T 6.6% A187D (GCC→GAC)  PP_3506 ← magnesium chelatase subunit ChII

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,976,3520GT 6.6% 190.7 / 2.3 60A187D (GCC→GAC) PP_3506magnesium chelatase subunit ChII
Reads supporting (aligned to +/- strand):  ref base G (23/33);  new base T (2/2);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTGCCGATCAGTACGAAGCGGGCGCT  >  NC_002947/3976203‑3976497
                                                                                                                                                     |                                                                                                                                                 
cggcgTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAgg                                                                                                                                                   <  1:633321/150‑1 (MQ=255)
cggcgTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAgg                                                                                                                                                   <  2:104963/150‑1 (MQ=255)
 ggcgTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGc                                                                                                                                                  <  1:120716/150‑1 (MQ=255)
  gcgTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGcc                                                                                                                                                 >  2:488075/1‑150 (MQ=255)
   cgTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCa                                                                                                                                                <  2:42504/150‑1 (MQ=255)
           cgcTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGAcggcga                                                                                                                                        <  1:1306976/150‑2 (MQ=255)
              tCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGa                                                                                                                                     >  1:523931/1‑150 (MQ=255)
                  tAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGa                                                                                                                                 >  2:9264/1‑150 (MQ=255)
                             cttgcCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCgg                                                                                                                      <  2:559524/147‑1 (MQ=255)
                             aGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCgg                                                                                                                      <  2:1126777/150‑1 (MQ=255)
                                   tGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCg                                                                                                                <  2:1348902/150‑1 (MQ=255)
                                   tGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCg                                                                                                                <  2:523931/150‑1 (MQ=255)
                                    ggTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCgg                                                                                                               >  2:1342596/1‑150 (MQ=255)
                                    ggTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCgg                                                                                                               <  1:504118/150‑1 (MQ=255)
                                        aTCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGTCCCAGTGTGCCTAGCACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTc                                                                                                           >  1:678103/1‑150 (MQ=255)
                                        aTCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTc                                                                                                           >  1:1382950/1‑150 (MQ=255)
                                         tCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTcc                                                                                                          >  2:641449/1‑150 (MQ=255)
                                              ggCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAAt                                                                                                                                   >  2:1152138/1‑120 (MQ=255)
                                              ggCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAAt                                                                                                                                   <  1:1152138/120‑1 (MQ=255)
                                              ggCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCaa                                                                                                     <  1:1132298/150‑1 (MQ=255)
                                              ggCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCaa                                                                                                     >  2:1352865/1‑150 (MQ=255)
                                                cAATGTCGGCCAAGGCCTGGCGCGCAGCGGTGCAGCGCTCGCGCAGTTGCGCCTGGGACTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGc                                                                                                   >  1:1328583/1‑150 (MQ=255)
                                                        gCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCg                                                                                                         >  2:142944/1‑136 (MQ=255)
                                                        gCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCg                                                                                                         <  1:142944/136‑1 (MQ=255)
                                                        gCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAgc                                                                                           <  1:834354/150‑1 (MQ=255)
                                                          cAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAgcgc                                                                                         <  2:892772/150‑1 (MQ=255)
                                                           aaGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCATAGGTCCAGGGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGcc                                                                                        >  1:1170243/1‑150 (MQ=255)
                                                           aaGGCCTGGCGCGCCGCGTGGCAGCGCTAGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGTCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGcc                                                                                        <  1:317219/150‑1 (MQ=255)
                                                               ccTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGt                                                                                    >  1:348290/1‑150 (MQ=255)
                                                                  ggCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCa                                                                                 <  1:9264/150‑1 (MQ=255)
                                                                   gcgcgcCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCaa                                                                                <  1:488075/150‑1 (MQ=255)
                                                                   gcgcgcCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCaa                                                                                <  1:1342596/150‑1 (MQ=255)
                                                                       gccgcGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCaa                                                                                                     >  2:1260751/1‑125 (MQ=255)
                                                                       gccgcGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCaa                                                                                                     <  1:1260751/125‑1 (MQ=255)
                                                                       gccgcGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGAGCCGAACGCCTGAGATTGTTTTTCCACGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCg                                                                            <  1:1236986/150‑1 (MQ=255)
                                                                             tGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACgg                                                                      <  1:641449/150‑1 (MQ=255)
                                                                                 aGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGa                                                                  <  1:407202/150‑1 (MQ=255)
                                                                                 aGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTAAAGAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGa                                                                  <  1:726778/150‑1 (MQ=255)
                                                                                                    ctgggcctgggcCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTc                                               >  1:1227011/1‑150 (MQ=255)
                                                                                                       ggcctgggccCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGcc                                            <  2:1382950/150‑1 (MQ=255)
                                                                                                               cccAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCaa                                                              <  1:486676/124‑1 (MQ=255)
                                                                                                               cccAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCaa                                                              >  2:486676/1‑124 (MQ=255)
                                                                                                                ccATTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCgg                                   <  2:1005337/150‑1 (MQ=255)
                                                                                                                ccAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGcc                                            <  1:1278395/141‑1 (MQ=255)
                                                                                                                ccAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGcc                                            >  2:1278395/1‑141 (MQ=255)
                                                                                                                 cAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCggg                                  >  2:1376973/1‑150 (MQ=255)
                                                                                                                       gcgcAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATg                            <  1:656562/150‑1 (MQ=255)
                                                                                                                       gcgcAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATg                            <  2:529552/150‑1 (MQ=255)
                                                                                                                         gcAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTc                                               >  1:836558/1‑129 (MQ=255)
                                                                                                                         gcAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTc                                               <  2:836558/129‑1 (MQ=255)
                                                                                                                          cAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTg                         <  2:1227011/150‑1 (MQ=255)
                                                                                                                          cAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTg                         >  1:787443/1‑150 (MQ=255)
                                                                                                                           aGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTGc                        <  2:678103/150‑1 (MQ=255)
                                                                                                                                       ggTCATTGTCGAAGTCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGa                                                                  <  1:904384/96‑1 (MQ=255)
                                                                                                                                       ggTCATTGTCGAAGTCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGa                                                                  <  1:194635/96‑1 (MQ=255)
                                                                                                                                       ggTCATTGTCGAAGTCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGa                                                                  >  2:904384/1‑96 (MQ=255)
                                                                                                                                       ggTCATTGTCGAAGTCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGa                                                                  >  2:194635/1‑96 (MQ=255)
                                                                                                                                        gTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTc                                                                    >  1:118833/1‑93 (MQ=255)
                                                                                                                                        gTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTc                                                                    <  2:118833/93‑1 (MQ=255)
                                                                                                                                            ttGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTGCCGATCAGTACGAAGCgg       <  1:1376973/150‑1 (MQ=255)
                                                                                                                                               tCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTGCCGATCAGTACGAAGCGGgcg    >  2:230228/1‑150 (MQ=255)
                                                                                                                                                 gAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTGCCGATCAGTACGAAGCGGGCGCt  >  2:498897/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGCCTGGGCCCAGTGTGCGCAGAACGCCTGAGGGTCATTGTCGAAGGCCAGACGGCGGCGAATGATCTGCTGCCGGGCCTCGGGCTCCGGCAAGCCCTCCAGCGCCACGTTCAAGCCGAAACGGTCGAGCAACTGCGGGCGCAGTTCGCCCTCTTCCGGGTTCATGGTGCCGATCAGTACGAAGCGGGCGCT  >  NC_002947/3976203‑3976497

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: