Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,021,449 A→G 6.6% M663V (ATG→GTG)  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,4490AG 6.6% 167.5 / 2.7 60M663V (ATG→GTG) PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (2/2);  total (31/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 = 9.63e-01

TGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCC  >  NC_002947/3021308‑3021589
                                                                                                                                             |                                                                                                                                            
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)
                                                                                                                                             |                                                                                                                                            
TGGGCAGCAGCGCGTGGTGAACATGGCAAGGCGCGCAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCC  >  NC_002947/3021308‑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: ‑