Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,021,444 G→A 6.5% R661Q (CGG→CAG)  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,021,4440GA 6.5% 176.5 / 2.1 61R661Q (CGG→CAG) PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

TTGGCTGACTGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCC  >  NC_002947/3021299‑3021589
                                                                                                                                                 |                                                                                                                                                 
ttGGCTGACTGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGca                                                                                                                                                >  2:214303/1‑149 (MQ=255)
         tGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCAGCATGTGCCggg                                                                                                                                       >  4:211274/1‑149 (MQ=255)
         tGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCAGCATGTGCCggg                                                                                                                                       >  2:389996/1‑149 (MQ=255)
         tGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTCGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTggg                                                                                                                                       >  5:156010/1‑149 (MQ=255)
           ggCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGc                                                                                                                                      <  2:328382/148‑1 (MQ=255)
           ggCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCg                                                                                                                                     >  7:659/1‑149 (MQ=255)
            gcagcaGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCAGCATGTGCCGGGcct                                                                                                                                    <  3:211274/149‑3 (MQ=255)
            gcagcaGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCAGCATGTGCCGGGcct                                                                                                                                    <  1:389996/149‑3 (MQ=255)
                    gcgTGGTGAACATGGCTAGGCTCGCAGGTATTACCAGCGTGTACTTGTCGTGCAGCCTGGGCGCCATGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGACCGGCACATGCTGGGCGTGAAGGTGc                                                                                                                            <  4:323248/149‑1 (MQ=255)
                    gcgTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGc                                                                                                                            <  4:240544/149‑1 (MQ=255)
                      gtggtgAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCca                                                                                                                          <  3:131368/149‑1 (MQ=255)
                         gtgAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGc                                                                                                                            <  5:434318/144‑1 (MQ=255)
                         gtgAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGc                                                                                                                            >  6:434318/1‑144 (MQ=255)
                            aaCATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAAttt                                                                                                                    >  6:298030/1‑149 (MQ=255)
                              cATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGc                                                                                                                  <  6:281654/149‑1 (MQ=255)
                               aTGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTTGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGcc                                                                                                                 <  5:86675/149‑1 (MQ=255)
                                tGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGc                                                                                                                                           >  1:87260/1‑122 (MQ=255)
                                tGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGc                                                                                                                                           <  2:87260/122‑1 (MQ=255)
                                    aaGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCtga                                                                                                            >  1:292951/1‑149 (MQ=255)
                                                tATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGc                                                                                                >  4:375579/1‑149 (MQ=255)
                                                tATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGc                                                                                                >  1:168762/1‑149 (MQ=255)
                                                 aTTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGcc                                                                                               >  8:317528/1‑149 (MQ=255)
                                                    aCCAGCGTGTACTGGTCGTGCAGCCTGGGCGCGCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAg                                                                                            >  1:67655/1‑149 (MQ=255)
                                                        gCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAgcggc                                                                                        <  2:155848/149‑1 (MQ=255)
                                                             tACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTg                                                                                   >  6:403265/1‑149 (MQ=255)
                                                               cTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTggc                                                                                 <  2:292951/149‑1 (MQ=255)
                                                                   tCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGgcggcg                                                                             >  8:376461/1‑149 (MQ=255)
                                                                     gTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGgt                                                                           <  3:375579/149‑1 (MQ=255)
                                                                     gTGCAGCCTGGGCGCCCTGATGCACGCCTAGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGACACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGgt                                                                           >  6:440604/1‑149 (MQ=255)
                                                                        cAGCCTGGGCGCCCTGATGCACTCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGtgtg                                                                        <  6:15320/149‑1 (MQ=255)
                                                                        cAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGtgtg                                                                        >  5:112627/1‑149 (MQ=255)
                                                                              gggCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACcgc                                                                  >  4:149281/1‑149 (MQ=255)
                                                                                gcgcCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCgcgc                                                                >  7:258927/1‑149 (MQ=255)
                                                                                   cccTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTgg                                                             <  8:106217/149‑1 (MQ=255)
                                                                                   cccTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTgg                                                             <  7:376461/149‑1 (MQ=255)
                                                                                         tGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAac                                                       <  8:659/149‑1 (MQ=255)
                                                                                           cACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAaca                                                      <  5:298030/148‑1 (MQ=255)
                                                                                            aCGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAacac                                                     >  4:6211/1‑148 (MQ=255)
                                                                                            aCGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAacac                                                     >  2:92147/1‑148 (MQ=255)
                                                                                              gCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGAt                                                  <  1:214303/149‑1 (MQ=255)
                                                                                              gCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGAt                                                  >  5:401935/1‑149 (MQ=255)
                                                                                                       tGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGcc                                         <  6:401935/149‑1 (MQ=255)
                                                                                                          tGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAgg                                      >  2:290714/1‑149 (MQ=255)
                                                                                                            gAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGgct                                    >  5:398765/1‑149 (MQ=255)
                                                                                                             aaaGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGgctg                                   <  6:398765/149‑1 (MQ=255)
                                                                                                                gcgcCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCtt                                <  7:317528/149‑1 (MQ=255)
                                                                                                                   ccGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGAc                             <  8:258927/149‑1 (MQ=255)
                                                                                                                    cGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGAcc                            <  2:67655/149‑1 (MQ=255)
                                                                                                                        gCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGAt                        >  6:91761/1‑149 (MQ=255)
                                                                                                                            aaGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGcc                    <  2:106247/149‑1 (MQ=255)
                                                                                                                                 aaCCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCgca               <  3:149281/149‑1 (MQ=255)
                                                                                                                                   ccTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCgcag              >  5:50773/1‑148 (MQ=255)
                                                                                                                                   ccTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCgcag              >  8:3111/1‑148 (MQ=255)
                                                                                                                                    cTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCgcag              >  6:192078/1‑147 (MQ=255)
                                                                                                                                      gcagcaGGCCCGGCACATGCTGGGCGTGAAGGTCCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTCGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGa          >  6:331920/1‑149 (MQ=255)
                                                                                                                                       cagcagGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGAc         <  5:440604/149‑1 (MQ=255)
                                                                                                                                         gcagGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGAcgc       <  2:228398/149‑1 (MQ=255)
                                                                                                                                             gCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGAcc   >  8:77601/1‑149 (MQ=255)
                                                                                                                                               ccGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGa     >  2:267590/1‑145 (MQ=255)
                                                                                                                                               ccGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGa     <  1:267590/145‑1 (MQ=255)
                                                                                                                                               ccGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGAccc  >  8:40512/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTGGCTGACTGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCC  >  NC_002947/3021299‑3021589

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