Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 473,116 T→C 6.0% E623E (GAA→GAG)  dnaG ← DNA primase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947473,1160TC 6.0% 193.2 / 2.5 67E623E (GAA→GAG) dnaGDNA primase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (2/2);  total (32/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.11e-01

AAAACAAGCCGAGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGTGCCATGCCAACGTGCGATC  >  NC_002947/472968‑473263
                                                                                                                                                    |                                                                                                                                                   
aaaaCAAGCCGAGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAgttgt                                                                                                                                                     >  1:12784/1‑149 (MQ=255)
aaaaCAAGCCGACGATTATACACGAGCATGGGCCTCACGCGCCAGAGGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCTCTTTGACGTGATTTCCTGAGCAgttgt                                                                                                                                                     <  4:93242/149‑1 (MQ=255)
   aCAAGCCGAGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTcc                                                                                                                                                  <  7:57389/149‑1 (MQ=255)
      aGCCGAGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAgg                                                                                                                                               <  7:326103/149‑1 (MQ=255)
        ccGAGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCt                                                                                                                                             >  1:327764/1‑149 (MQ=255)
          gAGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTgc                                                                                                                                           >  2:334761/1‑149 (MQ=255)
           aGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTgcg                                                                                                                                          >  7:127778/1‑149 (MQ=255)
                  taCCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGt                                                                                                                                   <  7:311365/149‑1 (MQ=255)
                        aGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAAGATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTg                                                                                                                                            >  1:252196/1‑134 (MQ=255)
                        aGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTg                                                                                                                                            <  2:252196/134‑1 (MQ=255)
                           aTGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTggcgggcg                                                                                                                          <  1:7704/149‑1 (MQ=255)
                               gCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACa                                                                                                                      <  3:104852/149‑1 (MQ=255)
                                 cTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAg                                                                                                                    >  8:25257/1‑149 (MQ=255)
                                 cTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAg                                                                                                                    <  8:190215/149‑1 (MQ=255)
                                    aCGCGCGAGATGAGGTCGGCGTTTGTGCCGGATGATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGAGGTGATTTCCTGAGCTGTTGTTCCAGGGTGCGCTCGCGTTGGCGGGCGGACAAGc                                                                                                                   <  1:11492/147‑1 (MQ=255)
                                    aCGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTa                                                                                                                 <  7:25257/149‑1 (MQ=255)
                                    aCGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCACTCGCGTTGGCGGGCGGACAAGc                                                                                                                   <  1:51326/147‑1 (MQ=255)
                                        gcCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTta                                                                                                             >  5:303787/1‑149 (MQ=255)
                                          cAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTTGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTtata                                                                                                           <  7:120273/149‑1 (MQ=255)
                                          cAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTtata                                                                                                           <  8:127778/149‑1 (MQ=255)
                                            gATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCt                                                                                                                  >  4:178257/1‑140 (MQ=255)
                                            gATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCt                                                                                                                  <  3:178257/140‑1 (MQ=255)
                                            gATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAgt                                                                                                         <  6:228571/149‑1 (MQ=255)
                                             aTGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAg                                                                                                                >  4:333612/1‑141 (MQ=255)
                                             aTGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAg                                                                                                                <  3:333612/141‑1 (MQ=255)
                                                aGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTc                                                                                                      <  8:105018/148‑1 (MQ=255)
                                                aGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTc                                                                                                      <  2:327764/148‑1 (MQ=255)
                                                       cGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACt                                                                                              <  8:62620/149‑1 (MQ=255)
                                                            gtgCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTta                                                                                                             >  2:315269/1‑129 (MQ=255)
                                                            gtgCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTta                                                                                                             <  1:315269/129‑1 (MQ=255)
                                                            gtgCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACtgttg                                                                                          <  7:295227/148‑1 (MQ=255)
                                                                 ggAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAg                                                                                    >  8:1309/1‑149 (MQ=255)
                                                                  gcaCATTTCGGCGCAGCACGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAgg                                                                                   <  6:303787/147‑1 (MQ=255)
                                                                        ttCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGAAAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCg                                                                             >  6:353645/1‑149 (MQ=255)
                                                                         tCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCg                                                                             <  2:12784/148‑1 (MQ=255)
                                                                                agcagGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGa                                                                      >  1:100880/1‑148 (MQ=255)
                                                                                               tGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTc                                                      >  6:162772/1‑149 (MQ=255)
                                                                                                gggTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGc                                                                                                                   <  6:245956/87‑1 (MQ=255)
                                                                                                gggTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGc                                                                                                                   >  5:245956/1‑87 (MQ=255)
                                                                                                gggTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCAtt                                                             >  1:321234/1‑141 (MQ=255)
                                                                                                gggTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCAtt                                                             <  2:321234/141‑1 (MQ=255)
                                                                                                   ttttttCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGc                                                   >  3:118798/1‑148 (MQ=255)
                                                                                                       ttCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAggg                                              >  4:214906/1‑149 (MQ=255)
                                                                                                       ttCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAggg                                              >  3:347342/1‑149 (MQ=255)
                                                                                                           tCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGcgcg                                          >  8:68433/1‑149 (MQ=255)
                                                                                                            cGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGcgcg                                          >  5:34459/1‑148 (MQ=255)
                                                                                                             gCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGcgcg                                          >  5:30832/1‑147 (MQ=255)
                                                                                                             gCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGt                                         <  4:347342/148‑1 (MQ=255)
                                                                                                               tGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAg                                       >  7:176090/1‑148 (MQ=255)
                                                                                                                    aaTTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAggcg                                 >  3:312760/1‑149 (MQ=255)
                                                                                                                    aaTTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAggcg                                 >  3:89724/1‑149 (MQ=255)
                                                                                                                     aTTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAggcgg                                <  4:312760/149‑1 (MQ=255)
                                                                                                                     aTTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAggcgg                                <  4:118798/149‑1 (MQ=255)
                                                                                                                      ttCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAggcggc                               >  3:70153/1‑149 (MQ=255)
                                                                                                                        cGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGcc                              <  2:100880/148‑1 (MQ=255)
                                                                                                                         gCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCt                            <  6:34459/149‑1 (MQ=255)
                                                                                                                           tttGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAggcgg                                >  6:284604/1‑143 (MQ=255)
                                                                                                                           tttGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAggcgg                                <  5:284604/143‑1 (MQ=255)
                                                                                                                            ttGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGtt                         >  4:46003/1‑149 (MQ=255)
                                                                                                                              gACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCg                       >  7:264294/1‑149 (MQ=255)
                                                                                                                                cGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGt                     <  3:214906/149‑1 (MQ=255)
                                                                                                                                      ttCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGtgccatg               >  5:56825/1‑149 (MQ=255)
                                                                                                                                         cTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGtgccatgcca            <  5:162772/149‑1 (MQ=255)
                                                                                                                                         cTGAGCAGTTGCTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAGCTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGtgccatgcc             >  1:24907/1‑148 (MQ=255)
                                                                                                                                         cTGAGCAGTTGCTCCAGGCTGCGCTCGCGTNGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAGCTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGtgccatgcc             >  1:30382/1‑148 (MQ=255)
                                                                                                                                          tGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCg                       >  4:60440/1‑137 (MQ=255)
                                                                                                                                          tGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCg                       <  3:60440/137‑1 (MQ=255)
                                                                                                                                            aGCAGTTGCTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAGCTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGTGCCATGCCAAc          <  2:24907/148‑1 (MQ=255)
                                                                                                                                            aGCAGTTGCTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAGCTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGTGCCATGCCAAc          <  2:30382/148‑1 (MQ=255)
                                                                                                                                                 ttgttCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGTGCCATGCCAACGTGCGa    <  5:353645/149‑1 (MQ=255)
                                                                                                                                                    ttCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGTGCCATGCCAACGTGCGATc  >  1:76205/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAAACAAGCCGAGGATTATACCCGAGCATGGGCCTCACGCGCCAGATGAGGTCGGCGTTTGTGCGGGAACATTTCGGCTCAGCAGGGCGAGGAGCTGGGTTTTTTCCTCGCTGGTCAATTCGCTTTGACGTGATTTCCTGAGCAGTTGTTCCAGGCTGCGCTCGCGTTGGCGGGCGGACAAGCTAGTTATAGTGTCGAAAAACTGTTGTTCAAGGTTGTCGGCCACGATAAGCCATTCCTTTTCCGCCAGGGCGCGTAGCAGGCGGCCCTGTTCGGTGCCATGCCAACGTGCGATC  >  NC_002947/472968‑473263

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