Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,647,083 C→A 6.2% V92L (GTG→TTG)  gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,647,0830CA 6.2% 209.7 / 3.8 65V92L (GTG→TTG) gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base A (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

GCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAG  >  NC_002947/1646937‑1647228
                                                                                                                                                   |                                                                                                                                                      
gctggcCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCAc                                                                                                                                                        <  4:345581/147‑1 (MQ=255)
  tggcCAGCGCGCAAGCGCCCGAAGCGACGACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGt                                                                                                                                                    >  1:206800/1‑149 (MQ=255)
  tggcCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGt                                                                                                                                                    >  1:201144/1‑149 (MQ=255)
        gcgcgcAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGccac                                                                                                                                          >  8:138719/1‑149 (MQ=255)
            gcAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGAGGGAGCACCAGAAAGCCAGCCGTGGCACCAGG‑‑‑‑TGGCAACAGt                                                                                                                                     >  1:307830/1‑149 (MQ=255)
             cAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCAACTGG‑‑‑‑TGccacggctg                                                                                                                                     <  4:256614/149‑6 (MQ=255)
                            aCCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCGCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGc                                                                                                                               >  1:266820/1‑140 (MQ=255)
                            aCCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGc                                                                                                                               <  2:266820/139‑1 (MQ=255)
                            aCCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGaa                                                                                                                     <  7:138719/149‑1 (MQ=255)
                                 tGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCg                                                                                                                       <  5:357752/142‑1 (MQ=255)
                                 tGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCg                                                                                                                       >  6:357752/1‑142 (MQ=255)
                                 tGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAggg                                                                                                                >  8:177828/1‑149 (MQ=255)
                                   ccACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAggg                                                                                                                <  8:330620/147‑1 (MQ=255)
                                   ccACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAggg                                                                                                                <  1:141224/147‑1 (MQ=255)
                                                tgcagtgTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGc                                                                                                                                             >  5:121270/1‑106 (MQ=255)
                                                tgcagtgTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGc                                                                                                                                             <  6:121270/106‑1 (MQ=255)
                                                       ttGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCcagcag                                                                                          <  5:44103/149‑1 (MQ=255)
                                                             ttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAg                                                                                    >  7:136498/1‑149 (MQ=255)
                                                              tGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAgg                                                                                   >  5:79199/1‑149 (MQ=255)
                                                                  ggctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCa                                                                                >  1:94585/1‑148 (MQ=255)
                                                                       ggCCGATCAGCGGGCGACGGGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAgg                                                                          <  1:354840/149‑1 (MQ=255)
                                                                           cGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGccc                                                                       <  8:332210/149‑1 (MQ=255)
                                                                            gATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTAAGGGGGCACCATTTACAGGCCAAATTAGAAAAGGGCCCCACCGGTGGTGCCCACCAGCCCCAGGCCaaaaag                                                                           >  1:114933/1‑144 (MQ=255)
                                                                                aGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAggg                                                                   <  6:79199/148‑1 (MQ=255)
                                                                                   gggCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGAt                                                               >  6:280386/1‑149 (MQ=255)
                                                                                    ggCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATc                                                              >  3:105243/1‑149 (MQ=255)
                                                                                    ggCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATc                                                              >  3:301555/1‑149 (MQ=255)
                                                                                             cccACGGCAGCAGCAACACCACGCCGATGGCGAACCATATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTCCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                     >  6:84987/1‑149 (MQ=255)
                                                                                             cccACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                     <  2:94585/149‑1 (MQ=255)
                                                                                             cccACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                     <  2:3762/149‑1 (MQ=255)
                                                                                              ccACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTggg                                                    <  5:280386/149‑1 (MQ=255)
                                                                                               cACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACTCGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGc                                                   <  2:307830/149‑1 (MQ=255)
                                                                                               cACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCTGGGCAGCAcc                                                                                      <  2:347305/114‑1 (MQ=255)
                                                                                               cACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCTGGGCAGCAcc                                                                                      >  1:347305/1‑114 (MQ=255)
                                                                                                   gcagcagcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagc                                               >  3:15055/1‑149 (MQ=255)
                                                                                                      gcagcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagc                                            <  7:177828/149‑1 (MQ=255)
                                                                                                        agcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagcag                                          <  8:136498/149‑1 (MQ=255)
                                                                                                           aaCACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGAc                                       >  7:29395/1‑149 (MQ=255)
                                                                                                           aaCACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGAc                                       >  4:78235/1‑149 (MQ=255)
                                                                                                             caccacGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGc                                     >  5:225582/1‑149 (MQ=255)
                                                                                                                cacGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGCTCTGGCGTTGGGCCAGCAGCAGGACGccc                                   >  4:262371/1‑148 (MQ=255)
                                                                                                                 acGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCTTAACGCCCg                                  >  1:128202/1‑148 (MQ=255)
                                                                                                                     cGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                             >  5:266801/1‑149 (MQ=255)
                                                                                                                       aTGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                             >  4:1850/1‑147 (MQ=255)
                                                                                                                       aTGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                             <  3:1850/147‑1 (MQ=255)
                                                                                                                            gAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGc                      <  4:105243/149‑1 (MQ=255)
                                                                                                                            gAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGc                      >  7:278182/1‑149 (MQ=255)
                                                                                                                               ccAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCaa                   >  6:121227/1‑149 (MQ=255)
                                                                                                                                cAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCaa                   >  1:234399/1‑148 (MQ=255)
                                                                                                                                  gATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagc                                            <  8:102041/121‑1 (MQ=255)
                                                                                                                                  gATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagc                                            >  7:102041/1‑121 (MQ=255)
                                                                                                                                      gccagccGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGc            <  2:234399/149‑1 (MQ=255)
                                                                                                                                      gccacccGTTGCACCAGA‑‑‑‑TGACACCAGCCCACGCCTACTTCGAACCGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCTATCTGGCGTTGGTCCAGCAGCAGGACGCCCGAAAGCGCAAAGCCAATACCGGc            <  8:282580/149‑1 (MQ=255)
                                                                                                                                       ccagccGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCGAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCg           >  7:229939/1‑149 (MQ=255)
                                                                                                                                       ccagccGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCGGGACGCCCGAAAGGGCAAACCCAATACCGGc            >  4:225296/1‑148 (MQ=255)
                                                                                                                                         agccGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCCGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGAt         <  2:119699/149‑1 (MQ=255)
                                                                                                                                           tggTGGCAACTGG‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagca                                              <  3:107269/108‑1 (MQ=255)
                                                                                                                                           tggTGGCAACTGG‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagca                                              >  4:107269/3‑110 (MQ=255)
                                                                                                                                           tggTGGCAACTGG‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGAt         >  3:167128/3‑147 (MQ=255)
                                                                                                                                             gTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                     >  2:41663/1‑101 (MQ=255)
                                                                                                                                             gTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                     <  1:41663/101‑1 (MQ=255)
                                                                                                                                             gTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAgg     <  5:84987/149‑1 (MQ=255)
                                                                                                                                               ggCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgta   <  2:206800/149‑1 (MQ=255)
                                                                                                                                                gCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgtag  >  6:28514/1‑149 (MQ=255)
                                                                                                                                                gCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgta   <  2:201144/148‑1 (MQ=255)
                                                                                                                                                gCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgta   <  5:121227/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGT‑‑‑‑TGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAG  >  NC_002947/1646937‑1647228

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: