Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,022,464 T→C 6.5% V1001A (GTG→GCG)  PP_2638 → cellulose synthase operon protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,022,4640TC 6.5% 175.3 / 2.2 62V1001A (GTG→GCG) PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base C (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

TCACTGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGTTCGGTGGGCGTGCCGGTCAG  >  NC_002947/3022317‑3022605
                                                                                                                                                   |                                                                                                                                              
tCACTGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCgtg                                                                                                                                               <  1:171190/149‑1 (MQ=255)
tCACTGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCgtg                                                                                                                                               <  6:280849/149‑1 (MQ=255)
   cTGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCgtggtg                                                                                                                                            <  6:198243/149‑1 (MQ=255)
   cTGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCgtggtg                                                                                                                                            >  4:109479/1‑149 (MQ=255)
    tGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCgtggtg                                                                                                                                            >  3:22701/1‑148 (MQ=255)
    tGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCgtggt                                                                                                                                             >  4:77457/1‑147 (MQ=255)
    tGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCgtggt                                                                                                                                             <  3:77457/147‑1 (MQ=255)
      ccAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCTTGGTGGaa                                                                                                                                         <  8:49828/149‑1 (MQ=255)
       cAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGaa                                                                                                                                         >  1:146913/1‑148 (MQ=255)
         aTGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAg                                                                                                                                      <  5:60020/149‑1 (MQ=255)
           ggCGCCACGCTCGGCCTGAGCTGGGACCAACGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCCTTTTCCACCACGCGGTGGTAAAAGc                                                                                                                                    >  5:78570/1‑149 (MQ=255)
            gcgcCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGcc                                                                                                                                   <  8:164527/149‑1 (MQ=255)
             cgcCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCCTTTTCCACCACGCGGTGGTAAAAGCcg                                                                                                                                  <  6:78570/149‑1 (MQ=255)
              gcCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCcgc                                                                                                                                 >  7:220904/1‑148 (MQ=255)
               ccACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCcgac                                                                                                                                >  4:47733/1‑149 (MQ=255)
                                  ggACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGc                                                                                                              <  6:361813/148‑1 (MQ=255)
                                   gACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGt                                                                                                            <  8:27378/149‑1 (MQ=255)
                                     ccAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGTGCGCCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGa                                                                                                          <  5:201213/149‑1 (MQ=255)
                                      cAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGAc                                                                                                         >  4:197926/1‑149 (MQ=255)
                                           tGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTTCCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGa                                                                                                    <  1:96562/149‑1 (MQ=255)
                                             gCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGAt                                                                                                   <  3:109479/148‑1 (MQ=255)
                                               gACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTTTACGGCGAAAACGTGGCCGACAACCAGCGGACCCTGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCa                                                                                                <  8:114441/149‑1 (MQ=255)
                                               gACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGCGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCCTTTTCCACCACGCGGTGGtaaaag                                                                                                                                     <  5:217957/112‑3 (MQ=38)
                                               gACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGAGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCCTTTTCCACCACGCGGTGGtaaaag                                                                                                                                     >  6:217957/1‑110 (MQ=38)
                                                aCAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCAc                                                                                               <  4:37512/149‑1 (MQ=255)
                                                 ctagggggTGTGGGCCAGCCTGGGGCAGCTTTGGGTGTAGGGGGATAAGCTGGCGGTGAACGTGCGGAGCCGGGCCTTGGGCGGCTTTTACCTGGGCGTGGTGGATAAGGCCGAGGAGGGTGTAGGTGTTGGCGTGACGCTGATGCa                                                                                                <  8:109226/143‑1 (MQ=255)
                                                  aaCGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGa                                                                                                    >  7:216909/1‑142 (MQ=255)
                                                  aaCGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGa                                                                                                    <  8:216909/142‑1 (MQ=255)
                                                  aaCGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTg                                                                                             >  2:316352/1‑149 (MQ=255)
                                                    cGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACCACCAGCGGACCCGGGCCATGGCCGGCTATTACCACCGCGTGCTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGgc                                                                                           >  5:208649/1‑149 (MQ=255)
                                                         ggtgGGCCAGCCTGGGGCACCATAGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCAt                                                                                      <  3:360173/148‑1 (MQ=255)
                                                           tgGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATg                                                                                     <  3:197926/148‑1 (MQ=255)
                                                              gCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACaa                                                                                 <  5:207980/149‑1 (MQ=255)
                                                              gCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGAACCGGGCCATGGCCGGCTATTACCACAGCGAGGTGGAAAAGGCGGACGAGCGTGTACGAGTTGGCGTGACGCTGATGCCCTGGCGGCATgcc                                                                                   >  2:351597/1‑145 (MQ=255)
                                                              gCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACCTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACGACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGAc                                                                                   >  4:199063/1‑147 (MQ=255)
                                                                cAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAg                                                                                >  1:302405/1‑148 (MQ=255)
                                                                             cATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGc                                                                                                                            <  6:362823/91‑1 (MQ=255)
                                                                             cATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGc                                                                                                                            >  5:362823/1‑91 (MQ=255)
                                                                             cATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACgt                                                                                                                    >  2:194049/1‑99 (MQ=255)
                                                                             cATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACgt                                                                                                                    <  1:194049/99‑1 (MQ=255)
                                                                                  gCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGgcgc                                                                                         >  5:345837/1‑121 (MQ=255)
                                                                                  gCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGgcgc                                                                                         <  6:345837/121‑1 (MQ=255)
                                                                                           cGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAgg                                                   >  2:97852/1‑149 (MQ=255)
                                                                                            gAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAggg                                                  >  3:109500/1‑149 (MQ=255)
                                                                                                aCGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGt                                              <  2:302405/149‑1 (MQ=255)
                                                                                                aCGGGGCCGACACCCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGCACCTGGGCGGCTACAGCCTGGTGCAGGGCGGt                                              <  2:121496/149‑1 (MQ=255)
                                                                                                  gTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTta                                            <  3:199063/149‑1 (MQ=255)
                                                                                                  gTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTta                                            <  4:222961/149‑1 (MQ=255)
                                                                                                           aaCCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGccc                                   >  6:247642/1‑149 (MQ=255)
                                                                                                           aaCCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGccc                                   >  6:248402/1‑149 (MQ=255)
                                                                                                             ccAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGc                                 <  3:47733/149‑1 (MQ=255)
                                                                                                                  ggACCCGGGCCATGGCCGGCTTTTACCAGCGGGTGGTGGATAAGTCCGACGTGCGTGTACGTGTTGGCGTGACGCTGATGCTCTGGCGCCATGACAAGGACCTGGGCGGCTACAGGCTGGGGCAGGGCGGTTACTACAGCCCGCAgcgc                            <  6:208649/149‑1 (MQ=255)
                                                                                                                      ccGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAgcg                             <  8:195506/144‑1 (MQ=255)
                                                                                                                      ccGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAgcg                             >  7:195506/1‑144 (MQ=255)
                                                                                                                      ccGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTa                          >  7:355258/1‑147 (MQ=255)
                                                                                                                      ccGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACa                        <  7:307031/149‑1 (MQ=255)
                                                                                                                       cGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTta                                            >  4:184707/1‑128 (MQ=255)
                                                                                                                       cGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTta                                            <  3:184707/128‑1 (MQ=255)
                                                                                                                         ggCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTTTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGtt                     <  4:44677/149‑1 (MQ=255)
                                                                                                                                gCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGTTCGGTgg               <  4:109500/148‑1 (MQ=255)
                                                                                                                                 ccGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACa                        >  2:3217/1‑138 (MQ=255)
                                                                                                                                 ccGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACa                        <  1:3217/138‑1 (MQ=255)
                                                                                                                                  cGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAg                               >  5:305216/1‑130 (MQ=255)
                                                                                                                                  cGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAg                               <  6:305216/130‑1 (MQ=255)
                                                                                                                                  cGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGTTCGGTGGGCg            >  4:38654/1‑149 (MQ=255)
                                                                                                                                         tACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGTTCGGTGGGCGTGCCg       <  1:97852/148‑1 (MQ=255)
                                                                                                                                          accaccGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCACCGCTACAGTTCGGTGGGCGTGCCGGt     >  4:203970/1‑148 (MQ=255)
                                                                                                                                            caccGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGTTCGGTGGGCGTGCCGGTCAg  >  5:71626/1‑149 (MQ=255)
                                                                                                                                            caccGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGTTCGGAGGGCGTGCCGGTCAg  >  7:232723/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCACTGCCAATGGCGCCACGCTCGGCCTGAGCTGGGACCAAGGTGGCGACAACGGGGTGTGGGCCAGCCTGGGGCACCATTGGCTGTACGGCGAAAACGTGGCCGACAACCAGCGGACCCGGGCCATGGCCGGCTTTTACCACCGCGTGGTGGAAAAGGCCGACGAGCGTGTACGTGTTGGCGTGACGCTGATGCACTGGCGCCATGACAAGGACCTGGGCGGCTACAGCCTGGGGCAGGGCGGTTACTACAGCCCGCAGCGCTACAGTTCGGTGGGCGTGCCGGTCAG  >  NC_002947/3022317‑3022605

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: