Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 852,166 G→C 5.3% A315A (GCG→GCC) ‡ hemA → glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947852,1660GC 5.3% 242.9 / 3.1 76A315A (GCG→GCC) ‡hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base C (2/2);  total (41/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 = 8.88e-01

ACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCCTGC  >  NC_002947/852027‑852314
                                                                                                                                           |                                                                                                                                                     
acaAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCgg                                                                                                                                                      >  5:143221/1‑141 (MQ=255)
acaAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCgg                                                                                                                                                      <  6:143221/141‑1 (MQ=255)
acaAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAggcg                                                                                                                                              >  3:196799/1‑149 (MQ=255)
     cccATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCg                                                                                                                                          >  6:321294/1‑148 (MQ=255)
      ccATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGaa                                                                                                                                        <  7:137276/149‑1 (MQ=255)
      ccATCTTCATGATCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCg                                                                                                                                          >  4:7950/1‑147 (MQ=255)
        aTCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCTGGGCCGCGCCCTggcggctcttga                                                                                                                                      <  3:22913/149‑7 (MQ=255)
           ttCATGGTCGACATCGCTGTACCTCGCGATAGCGCAACCGAAGTCGGCGAGCTGGATGATGTCTCCCTGTACACCGTCGACGACCTGCACGATGGAGGGGCGGAACACCTCAAGCGCCGCCAGGGCGCGGCCCCGGCGGCCGCAGCGCt                                                                                                                                   >  6:245900/1‑149 (MQ=255)
             cATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGc                                                                                                                                    <  7:175968/146‑1 (MQ=255)
              aTGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCCCGATGTAGGGGCGGAAACCCTCACGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTgg                                                                                                                                 >  6:237168/1‑148 (MQ=255)
                    gACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGt                                                                                                                          >  7:241161/1‑149 (MQ=255)
                     aCATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTg                                                                                                                         >  3:214327/1‑149 (MQ=255)
                      cATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTgg                                                                                                                       <  4:214327/149‑1 (MQ=255)
                              tACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGaa                                                                                                               <  8:241161/149‑1 (MQ=255)
                               aCCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCg                                                                                                                 >  2:109470/1‑146 (MQ=255)
                                      gATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTgg                                                                                                                       >  8:58501/1‑133 (MQ=255)
                                      gATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTgg                                                                                                                       <  7:58501/133‑1 (MQ=255)
                                        tatCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCAtgc                                                                                                     >  7:146719/1‑149 (MQ=255)
                                           cGAAACCGAACTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACCACCTCCACGATGAAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCAAGGCGGCCGAACAGCTGCTATCCGTGGGCGCCGAACACTTCAtccacc                                                                                                  >  8:200451/1‑144 (MQ=255)
                                            gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGctgcg                                                                                                 >  2:189252/1‑149 (MQ=255)
                                            gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGctgcg                                                                                                 >  6:43315/1‑149 (MQ=255)
                                            gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGctgcg                                                                                                 >  2:157612/1‑149 (MQ=255)
                                            gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGctgcg                                                                                                 >  4:301412/1‑149 (MQ=255)
                                              aaCCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGctgcgcc                                                                                               >  1:80906/1‑149 (MQ=255)
                                                   aaGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGctgcgcctgcgc                                                                                          <  5:321294/149‑1 (MQ=255)
                                                    aGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGCGGCGCCGAAGACTTCATGctgcgcctgcgc                                                                                          <  1:189252/149‑1 (MQ=255)
                                                        ggCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACt                                                                                     >  6:46663/1‑149 (MQ=255)
                                                         gCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTg                                                                                    <  1:162248/149‑1 (MQ=255)
                                                             gCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTg                                                                                >  7:80228/1‑149 (MQ=255)
                                                                 gatgatGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGt                                                                            <  5:43315/149‑1 (MQ=255)
                                                                 gatgatGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGt                                                                            <  3:7950/149‑1 (MQ=255)
                                                                     atgTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                                        <  5:82988/149‑1 (MQ=255)
                                                                     atgTCTACCTGTAAACCGTAGACGACCTGCACGATGTAGTGTCGGAAAACCTAAAGAGCCGCCATGGCGCGGCCCAGGAGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCATGCGCGAAc                                                                                      >  2:293590/1‑135 (MQ=255)
                                                                     atgTCTACATGTACACCGTCGACGCCCTGCACTATGTAGTGGCGTAAAACCTCAAGAGCCGCCATGTCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAATACTTCATGCTGCGCCTGCGCTAAc                                                                                      <  1:293590/135‑1 (MQ=255)
                                                                        tCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGt                                                                      <  7:85495/148‑1 (MQ=255)
                                                                         cTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTgg                                                                          <  8:60576/143‑1 (MQ=255)
                                                                         cTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTgg                                                                          >  7:60576/1‑143 (MQ=255)
                                                                           aCCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCt                                                                   >  2:159228/1‑148 (MQ=255)
                                                                             cTGTACCCCGTCGACGCCCTGCCCGAGGGAGTGGCGGAAACCCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGGATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCGGCAGTGGACGTGCTCAg                                                                <  5:245900/149‑1 (MQ=255)
                                                                              tGTACTCGGTCGACGACCTGCACGTTGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAgg                                                               <  8:146719/149‑1 (MQ=255)
                                                                               gTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTAt                                                                                                                              <  3:166810/86‑1 (MQ=255)
                                                                               gTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTAt                                                                                                                              >  4:166810/1‑86 (MQ=255)
                                                                               gTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGAGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                              >  7:244402/1‑149 (MQ=255)
                                                                                tACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAgg                                                               <  7:200451/147‑1 (MQ=255)
                                                                                tACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGc                                                             <  8:80228/149‑1 (MQ=255)
                                                                                  cacCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGcc                                                            >  5:289225/1‑148 (MQ=255)
                                                                                   acCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTa                                                          >  2:52340/1‑149 (MQ=255)
                                                                                         gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGc                                                      >  3:65116/1‑147 (MQ=255)
                                                                                         gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGc                                                      >  1:140205/1‑147 (MQ=255)
                                                                                         gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCa                                                    <  3:301412/149‑1 (MQ=255)
                                                                                              ccTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCa                                                 >  7:234915/1‑147 (MQ=255)
                                                                                               cTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGCCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAg                                              >  3:147318/1‑149 (MQ=255)
                                                                                                  cACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGt                                             >  3:286541/1‑147 (MQ=255)
                                                                                                    cGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAgc                                         <  1:52340/149‑1 (MQ=255)
                                                                                                     gATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGTGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAgcg                                        <  8:244402/149‑1 (MQ=255)
                                                                                                     gATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAgcg                                        >  1:95320/1‑149 (MQ=255)
                                                                                                           gTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAg                                                                >  4:240230/1‑119 (MQ=255)
                                                                                                           gTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAg                                                                <  3:240230/119‑1 (MQ=255)
                                                                                                              gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                               >  1:218791/1‑149 (MQ=255)
                                                                                                              gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCACGCTGCGCCTGCGCGAACTGGCTGCAGTGGAGGTGCTCATGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                               >  5:222252/1‑149 (MQ=255)
                                                                                                              gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCACGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCATGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                               <  6:222252/149‑1 (MQ=255)
                                                                                                                   aaaCCTCAAGAGCCGCCAGGGCGCGGCCCATGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACgaag                          <  6:61975/149‑1 (MQ=255)
                                                                                                                      ccTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAAc                       <  4:147318/149‑1 (MQ=255)
                                                                                                                       cTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACt                      >  1:312679/1‑149 (MQ=255)
                                                                                                                       cTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACt                      >  2:316172/1‑149 (MQ=255)
                                                                                                                       cTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCGGGGATACGTGGGCGCCGACGACTTCATGCTGCGCCTGCGCGTACTGGATGCTGTGGCAGTGCGCAGGGCCTAGAGCGCGCGAAGTGAGCTCCTGCGCGACGACGAGCt                      >  6:188327/1‑149 (MQ=255)
                                                                                                                          aaGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTg                     <  1:170660/147‑1 (MQ=255)
                                                                                                                           agagCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCaa                  <  8:234915/149‑1 (MQ=255)
                                                                                                                            gagCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                                        >  7:275548/1‑94 (MQ=255)
                                                                                                                            gagCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                                        <  8:275548/94‑1 (MQ=255)
                                                                                                                             agCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCaaaa                <  2:95320/149‑1 (MQ=255)
                                                                                                                             agCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCaaaa                <  2:312679/149‑1 (MQ=255)
                                                                                                                               ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                              <  1:181103/101‑1 (MQ=255)
                                                                                                                               ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                              >  2:181103/1‑101 (MQ=255)
                                                                                                                               ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAgg              >  7:272782/1‑149 (MQ=255)
                                                                                                                                cgccAGGGCGCGGCGCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGc             <  4:196799/149‑1 (MQ=255)
                                                                                                                                  ccAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCCCCGAACACTTCATGCTGCGCCTGCCCGAACTCGCTGCAGTGGAc                                                                        >  6:222980/1‑88 (MQ=255)
                                                                                                                                  acAGGGCGCGGCCCTGGCGGCCTAATAGCTGGTATCGGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCTGTGGAc                                                                        <  5:222980/87‑1 (MQ=255)
                                                                                                                                           ggCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCctgc  <  4:65116/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
ACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGGGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCCTGC  >  NC_002947/852027‑852314

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