Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,368,670 G→T 7.8% R14S (CGC→AGC)  greA ← transcription elongation factor GreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,368,6700GT 7.8% 178.4 / 2.1 60R14S (CGC→AGC) greAtranscription elongation factor GreA
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base T (2/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGTGCAGTAAA  >  NC_002947/5368530‑5368818
                                                                                                                                            |                                                                                                                                                    
gACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCTCGCGCCCt                                                                                                                                              <  5:248263/149‑1 (MQ=255)
gACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCt                                                                                                                                              <  5:54937/149‑1 (MQ=255)
gACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCt                                                                                                                                              <  6:338811/149‑1 (MQ=255)
   cATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGa                                                                                                                                           <  2:287777/149‑1 (MQ=255)
       cccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCTCGCGCCCTGGaaga                                                                                                                                        <  5:338391/148‑4 (MQ=255)
       cccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCTCGCGCCCTGGaag                                                                                                                                         >  5:58801/1‑145 (MQ=255)
         cTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCTCGCGCCCTGGaagaaga                                                                                                                                     <  4:239030/149‑7 (MQ=255)
               ctcACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCggg                                                                                                                                >  8:316582/1‑148 (MQ=255)
               ctcACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCgg                                                                                                                                 <  1:37520/147‑1 (MQ=255)
                tcACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGt                                                                                                                               >  7:23866/1‑148 (MQ=255)
                tcACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTa                                                                                                                              >  5:65223/1‑149 (MQ=255)
                  aCGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGt                                                                                                                               <  8:23866/146‑1 (MQ=255)
                     ggcggcATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTgg                                                                                                                         <  8:70053/149‑1 (MQ=255)
                          cATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGTCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACg                                                                                                                                         <  6:152264/128‑1 (MQ=255)
                          cATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACg                                                                                                                                         >  5:152264/1‑128 (MQ=255)
                           aTGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGc                                                                                                                   >  8:1949/1‑149 (MQ=255)
                                 ttCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGc                                                                                                             >  4:117202/1‑149 (MQ=255)
                                    ggCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                          <  7:316582/149‑1 (MQ=255)
                                      cGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCCCCAAGGGCCTGCCTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACCGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCa                                                                                                        <  6:41391/149‑1 (MQ=255)
                                      cGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCa                                                                                                        <  3:199796/149‑1 (MQ=255)
                                              ttGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGATTTCAGTCAtgaa                                                                                                >  2:323938/1‑147 (MQ=255)
                                                gAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgca                                                                                                <  4:6111/147‑1 (MQ=255)
                                                   gTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatg                                                                                              >  6:14853/1‑146 (MQ=255)
                                                      aCCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATc                                                                                        >  3:125730/1‑149 (MQ=255)
                                                              cACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGc                                                                                <  5:14853/149‑1 (MQ=255)
                                                                 gcgcTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTcc                                                                                                    >  2:86219/1‑126 (MQ=255)
                                                                 gcgcTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTcc                                                                                                    <  1:86219/126‑1 (MQ=255)
                                                                     ttCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACg                                                                                                                                         >  8:290845/1‑85 (MQ=255)
                                                                     ttCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACg                                                                                                                                         <  7:290845/85‑1 (MQ=255)
                                                                     ttCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGAGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGAAAGCGACGAAc                                                                         >  3:324864/1‑149 (MQ=255)
                                                                          cAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCGTCGGGTATTTGGTAATTCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCt                                                                    >  7:310413/1‑149 (MQ=255)
                                                                               gCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTc                                                                <  1:1363/148‑1 (MQ=255)
                                                                                      tCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACg                                                                        >  5:96055/1‑133 (MQ=255)
                                                                                      tCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACg                                                                        <  6:96055/133‑1 (MQ=255)
                                                                                        aGGCGCGGTCGCTCGGTCTTGCTCAGGAATAGGTGGTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCGTCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGaa                                                       <  8:310413/148‑1 (MQ=255)
                                                                                                  gcTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCtt                                            <  4:158517/149‑1 (MQ=255)
                                                                                                    tCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCg                                          >  8:23465/1‑149 (MQ=255)
                                                                                                    tCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCg                                          >  2:260656/1‑149 (MQ=255)
                                                                                                            gCTCAGGAACTGGTCCTCTTCTTGCAGGGCCCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGAt                                     <  4:125730/146‑1 (MQ=255)
                                                                                                            gCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGAt                                     <  1:125428/146‑1 (MQ=255)
                                                                                                            gCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATgg                                   <  3:117202/148‑1 (MQ=255)
                                                                                                             cTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGc                                  <  6:65223/148‑1 (MQ=255)
                                                                                                                aGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGCCGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCa                              >  8:214253/1‑149 (MQ=255)
                                                                                                                aGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCa                              >  3:219902/1‑149 (MQ=255)
                                                                                                                 ggAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCAc                             >  7:24159/1‑149 (MQ=255)
                                                                                                                   aaCAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGGATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCGAGCGACGAACGGGCTGTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCGTCAc                             >  2:241499/1‑147 (MQ=255)
                                                                                                                        gTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCg                       >  6:313205/1‑148 (MQ=255)
                                                                                                                            tcttctTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATg                  >  3:226199/1‑149 (MQ=255)
                                                                                                                             cttcttCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATgg                 >  2:140008/1‑149 (MQ=255)
                                                                                                                             cttcttCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATgg                 >  2:20552/1‑149 (MQ=255)
                                                                                                                             cttcttCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATgg                 >  6:326279/1‑149 (MQ=255)
                                                                                                                              ttgttcCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTTTTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATg                  <  3:230495/144‑1 (MQ=255)
                                                                                                                              ttcttcCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATg                  <  8:24159/147‑1 (MQ=255)
                                                                                                                                cttcCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTa               <  1:20552/148‑1 (MQ=255)
                                                                                                                                 ttcCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAgt             >  6:289012/1‑149 (MQ=255)
                                                                                                                                  tcCAGGGCGCTCGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAgtg            >  8:21648/1‑149 (MQ=255)
                                                                                                                                      gggCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGTGc        <  5:326279/149‑1 (MQ=255)
                                                                                                                                      gggCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGCTCTTTTCCGCACCGAATTACAGCGCTTCGCAGATGGCTTCACCAGACGCAATGGTAgtcgtg         >  8:83159/1‑148 (MQ=255)
                                                                                                                                        gcgcGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAgtggtg         >  4:320/1‑146 (MQ=255)
                                                                                                                                        gcgcGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGTGCa       <  3:320/148‑1 (MQ=255)
                                                                                                                                          gcgAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGTGCAGt     >  6:14725/1‑148 (MQ=255)
                                                                                                                                            gAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAAGTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGGGGTGCAGTaaa  >  4:262591/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGTGCAGTAAA  >  NC_002947/5368530‑5368818

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: