Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,368,670 G→T 6.3% 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 6.3% 201.8 / 2.0 63R14S (CGC→AGC) greAtranscription elongation factor GreA
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base T (2/2);  total (30/33)
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.71e-01

GCGCTTCGACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGT  >  NC_002947/5368523‑5368810
                                                                                                                                                   |                                                                                                                                            
gcgcTTCGACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGa                                                                                                                                             <  4:175739/149‑1 (MQ=255)
  gcTTCGACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAgc                                                                                                                                           <  2:52729/149‑1 (MQ=255)
               ccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCCGGCGCGGGCGCTCGGGCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGt                                                                                                                               >  6:195391/1‑148 (MQ=255)
               ccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGt                                                                                                                               >  2:22145/1‑148 (MQ=255)
               ccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCGTCCAGGGCGCTCGCGCCCTGGaagaa                                                                                                                               <  5:112884/148‑5 (MQ=255)
               ccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCGTCCAGGGCGCTCGCGCCCTGGaag                                                                                                                                 >  6:112884/1‑144 (MQ=255)
               ccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTAACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCGTCCAGGGCGCTCGCGCCCTGGaagaag                                                                                                                              <  5:66622/149‑6 (MQ=255)
                   ctcctcACGGGCGGCATGGTATTCGGCTTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCTAGGGCGCGAGCGCCCTGGACGGTCATCg                                                                                                                          <  3:163135/149‑1 (MQ=255)
                     cctcACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCggg                                                                                                                        >  4:70473/1‑149 (MQ=255)
                      ctcACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGt                                                                                                                       <  5:248610/149‑1 (MQ=255)
                           gggCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTg                                                                                                                  <  5:66018/149‑1 (MQ=255)
                             gcggcATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGAGCTCTGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGAGCGAGCGCCCTGGACGGTCATAGGGTATTTgg                                                                                                                 >  3:215861/1‑148 (MQ=255)
                                  aTGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGc                                                                                                           >  1:136789/1‑149 (MQ=255)
                                    ggTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCt                                                                                                          >  5:87358/1‑148 (MQ=255)
                                    ggTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCt                                                                                                          >  5:225165/1‑148 (MQ=255)
                                      tATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCAt                                                                                                       <  3:40697/149‑1 (MQ=255)
                                         tCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCt                                                                                                    >  6:251682/1‑149 (MQ=255)
                                           ggCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                  <  8:160574/149‑1 (MQ=255)
                                             cGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCa                                                                                                <  2:136789/149‑1 (MQ=255)
                                              gTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAg                                                                                               >  3:137916/1‑149 (MQ=255)
                                                   ccTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCt                                                                                           <  2:44747/148‑1 (MQ=255)
                                                     ttGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCt                                                                                                                                      <  3:156025/103‑1 (MQ=255)
                                                     ttGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCt                                                                                                                                      >  4:156025/1‑103 (MQ=255)
                                                     ttGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGAGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCTTCGGGTATTTGGTAATGCTCATGCTTTCAGTACtgcg                                                                                        >  2:217711/1‑148 (MQ=255)
                                                       gAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatg                                                                                      >  8:143224/1‑149 (MQ=255)
                                                        aGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatgc                                                                                     <  2:80230/149‑1 (MQ=255)
                                                        aGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatg                                                                                      <  7:198197/148‑1 (MQ=255)
                                                          gTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatg                                                                                      >  7:198215/1‑146 (MQ=255)
                                                             aCCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATc                                                                                >  2:28060/1‑149 (MQ=255)
                                                                 aGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTg                                                                             >  8:170188/1‑148 (MQ=255)
                                                                    tCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAg                                                                         >  2:194114/1‑149 (MQ=255)
                                                                      aCGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAg                                                                         <  8:113558/147‑1 (MQ=255)
                                                                      aCGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGACATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGc                                                                        >  3:142867/1‑148 (MQ=255)
                                                                           cTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGaa                                                                  <  6:225165/149‑1 (MQ=255)
                                                                            ttCACCAATGGCCTGACTCAGGCGCGGGGGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATcc                                                                               >  1:151785/1‑135 (MQ=255)
                                                                            ttCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCa                                                                                                >  5:131882/1‑118 (MQ=255)
                                                                            ttCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCa                                                                                                <  6:131882/118‑1 (MQ=255)
                                                                            ttCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATcc                                                                               <  2:151785/135‑1 (MQ=255)
                                                                            ttCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAAc                                                                 >  6:7172/1‑149 (MQ=255)
                                                                                      gCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTcc                                                       <  8:187266/149‑1 (MQ=255)
                                                                                      gCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTc                                                        <  1:194114/148‑1 (MQ=255)
                                                                                        cTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCg                                                      >  3:25408/1‑148 (MQ=255)
                                                                                           aCTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGAc                                                   <  5:251682/148‑1 (MQ=255)
                                                                                                  cgcgGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAACCGACGAACGGTCTTTTCCGGACCGAATTTc                                           <  1:22145/149‑1 (MQ=255)
                                                                                                    gtggCGCTCGGTCTTGCTCTGGAACAGGTTCTCTTCTTCCCGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAg                                         <  5:228135/147‑1 (MQ=255)
                                                                                                    cgGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAg                                         <  4:142867/149‑1 (MQ=255)
                                                                                                           tCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                  <  4:102351/85‑1 (MQ=255)
                                                                                                           tCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                  >  3:102351/1‑85 (MQ=255)
                                                                                                                    cTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCt                         >  8:76618/1‑149 (MQ=255)
                                                                                                                       aGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCtt                        <  4:137916/147‑1 (MQ=255)
                                                                                                                               gtcctctACTTCCAGGGAGCGAGCGCCCTGGACCGTCAACGGCTATTTGGTAATGCTCATGCTTTCAGTCCTGCATCCACAACCTGCAAGCGACGACCCGTCTGTTCCGGACACAATTTCAGCGCTTCGCAGATCGCATGACCAGCCg               >  6:158777/4‑148 (MQ=255)
                                                                                                                                  ctctTCTTCGAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCaa            <  5:7172/148‑1 (MQ=255)
                                                                                                                                  ctctTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCaa            <  7:76618/148‑1 (MQ=255)
                                                                                                                                   tcttctTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCAc                     <  7:126173/138‑1 (MQ=255)
                                                                                                                                   tcttctTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCAc                     >  8:126173/1‑138 (MQ=255)
                                                                                                                                    cttcttCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATgg         >  3:2343/1‑149 (MQ=255)
                                                                                                                                       cttcCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTa       <  7:170188/148‑1 (MQ=255)
                                                                                                                                       cgtcCAGGGCGCTCGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAg      >  2:183323/3‑149 (MQ=255)
                                                                                                                                          ccAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAgtg    <  6:225632/148‑1 (MQ=255)
                                                                                                                                          ccAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAgtg    >  5:225632/1‑148 (MQ=255)
                                                                                                                                           cAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAgtggt  <  4:25408/149‑1 (MQ=255)
                                                                                                                                            aGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTTGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATg          <  4:75316/140‑1 (MQ=255)
                                                                                                                                            aGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATg          >  3:75316/1‑140 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCGCTTCGACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGT  >  NC_002947/5368523‑5368810

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: