Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,261,969 C→G 7.9% E265Q (GAG→CAG) ‡ PP_3734 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,261,9690CG 7.9% 176.6 / 2.5 58E265Q (GAG→CAG) ‡PP_3734ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base G (2/3);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCAGTTCAAGATAGCGGGC  >  NC_002947/4261825‑4262116
                                                                                                                                                 |                                                                                                                                                   
aCGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTTCTCCAgc                                                                                                                                                 <  2:279932/149‑1 (MQ=255)
aCGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAgc                                                                                                                                                 >  5:161940/1‑149 (MQ=255)
     gCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAgcgc                                                                                                                                               >  1:21368/1‑146 (MQ=255)
      ccAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGGAGCAGCGccgtg                                                                                                                                           <  3:6012/149‑4 (MQ=255)
       cAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGAc                                                                                                                                          <  2:21368/149‑1 (MQ=255)
        aaTGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCCCGCGCGCAGGTTACGGATCAGCCCCCGGCGCTCCTCCAGCGCCAGAc                                                                                                                                          >  1:241333/1‑149 (MQ=255)
        aaTGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGAcc                                                                                                                                         >  7:83194/1‑149 (MQ=255)
        aaTGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGAcc                                                                                                                                         >  7:302381/1‑149 (MQ=255)
             gTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTg                                                                                                                                    >  6:148677/1‑149 (MQ=255)
               acccACTGGCCACCCCGGCCAGCCCACCCAGGGCCGCGAACAAGCCCAGCTCCAGCACCCGGGCGCACAGCCGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTgcg                                                                                                                                  <  8:313395/147‑1 (MQ=255)
                aCCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTgcg                                                                                                                                  <  3:170534/148‑1 (MQ=255)
                  cACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTgcggcg                                                                                                                               >  7:179967/1‑149 (MQ=255)
                       gCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTgg                                                                                                                            <  5:148677/147‑1 (MQ=255)
                         cACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCaa                                                                                                                        <  1:200123/149‑1 (MQ=255)
                          aCCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAAt                                                                                                                       <  3:109530/149‑1 (MQ=255)
                                     cccacccaGTGCCGCGAACAAGCCCAGCTCCGGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCcg                                                                                                             >  5:166438/1‑148 (MQ=255)
                                     cccacccaGTGCCGCGAACAAGCCCAGCTCCGGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCcg                                                                                                             <  6:166438/148‑1 (MQ=255)
                                            aGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGaa                                                                                                     <  6:161940/149‑1 (MQ=255)
                                             gTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAAc                                                                                                    <  8:179967/149‑1 (MQ=255)
                                                ccGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTACAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGc                                                                                                  >  1:241289/1‑148 (MQ=255)
                                                 cgcgAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCca                                                                                                <  8:302381/149‑1 (MQ=255)
                                                  gcgAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCcag                                                                                               >  3:1603/1‑149 (MQ=255)
                                                         aaGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGcc                                                                                         >  7:254091/1‑149 (MQ=255)
                                                             ccAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAgc                                                                                     >  2:19473/1‑149 (MQ=255)
                                                                gCTCCAGCACCAGGGCGCACAGCTGGGTCTTCAGGATGACGCCGCACGCGCGCAGGGTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGcgcg                                                                                  <  6:59236/149‑1 (MQ=255)
                                                                  tCCAGCACCAGTGCGCACAGCAGGGTCTTCAGGCTGACTACGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGGTCGAGCGCCAGACCGATTGCGGCGTTGGCAATGAACAATCCGACGACGAACGCCAGCAGGCCGAGCTCGGt                                                                                <  2:241333/149‑1 (MQ=255)
                                                                   ccAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGt                                                                                >  1:201170/1‑148 (MQ=255)
                                                                    cAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCcagcag                                                                                            <  7:224048/135‑1 (MQ=255)
                                                                    cAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCcagcag                                                                                            >  8:224048/1‑135 (MQ=255)
                                                                    cAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGt                                                                                >  4:280354/1‑147 (MQ=255)
                                                                        aCCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTc                                                                          >  6:23209/1‑149 (MQ=255)
                                                                          cAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGcgcg                                                                                  >  5:276826/1‑139 (MQ=255)
                                                                          cAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGcgcg                                                                                  <  6:276826/139‑1 (MQ=255)
                                                                          cAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGAGCAGGTTACGGATCAGCCAACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAACCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCacgac                                                                          >  7:209098/1‑143 (MQ=255)
                                                                               cgcACAGCAGGGTCTTCAGGCTGACGCCTCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCTGTCAGGTCGAGGTgg                                                                   <  8:83194/149‑1 (MQ=255)
                                                                                  acaGCAGGGTCTTCAGGCTGACGCCTCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTACAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGCGACCGGACAGGTCGAGTTGGaaa                                                                >  7:114364/1‑149 (MQ=255)
                                                                                   cagcagGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAAc                                                               <  2:241289/149‑1 (MQ=255)
                                                                                     gcagGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTg                                                             >  3:229076/1‑149 (MQ=255)
                                                                                      cagGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGt                                                            >  2:88504/1‑149 (MQ=255)
                                                                                               cAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCa                                                      <  3:280354/146‑1 (MQ=255)
                                                                                                    tGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTg                                               >  8:32093/1‑148 (MQ=255)
                                                                                                      aCGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCaa                                            <  6:90858/149‑1 (MQ=255)
                                                                                                        gccgcACGCTCGCAGGTTACGGATCAGCCCACGGCGCTGCTGCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGc                                          <  8:114364/149‑1 (MQ=255)
                                                                                                            cACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGccacc                                      <  8:254091/149‑1 (MQ=255)
                                                                                                            cACGCGCGCAGGGTACGGATCAGCCCACGGCGCTGCTCGAGCGCCAGAGCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGccacc                                      <  7:257513/149‑1 (MQ=255)
                                                                                                                   gcAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTcc                               >  4:237297/1‑149 (MQ=255)
                                                                                                                       gTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGcc                           >  2:63446/1‑149 (MQ=255)
                                                                                                                          aCGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCcgcg                        >  8:178186/1‑149 (MQ=255)
                                                                                                                            ggATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGc                      >  2:89780/1‑149 (MQ=255)
                                                                                                                               tCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAATCTGTCGGTAAGGCGCTGCAAGCCTCCTGCGTCCTGCCGCGGCTg                    >  4:178423/1‑148 (MQ=255)
                                                                                                                                 aGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCa                  <  5:23209/148‑1 (MQ=255)
                                                                                                                                 aGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCa                  <  7:32093/148‑1 (MQ=255)
                                                                                                                                     tctggCGCTGGAGCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCa                  >  8:261926/4‑144 (MQ=255)
                                                                                                                                     tctggCGCTGGAGCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCa                  <  7:261926/141‑1 (MQ=255)
                                                                                                                                     cACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCAGTTc              >  6:199381/1‑148 (MQ=255)
                                                                                                                                      ctggCGCTGGAGCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCAGTTCaa            >  7:270844/3‑149 (MQ=255)
                                                                                                                                       cGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCAGTTCAAg           <  2:201170/149‑1 (MQ=255)
                                                                                                                                             tgctCCATCGCCAGTCCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCGAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGc                      <  5:284088/132‑1 (MQ=255)
                                                                                                                                                 ccAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCAGTTCAAGATAGCGGGc  >  3:294598/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCGCACAGCAGGGTCTTCAGGCTGACGCCGCACGCGCGCAGGTTACGGATCAGCCCACGGCGCTGCTCCAGCGCCAGACCGATTGCGGCGTGGGCAATGAACAAGCCGACGACGAACGCCAGCAGGCCGAGCGCGGTCAGGTCGAGGTGGAAACTGTCGGTCAGGCGCTGCAAGCCACCATCGTCCTGCCGCGGCTGCAGTTCAAGATAGCGGGC  >  NC_002947/4261825‑4262116

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: