Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 632,051 C→G 6.3% A464P (GCG→CCG)  aldB‑I ← aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947632,0510CG 6.3% 198.4 / 4.2 63A464P (GCG→CCG) aldB‑Ialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base G (2/2);  total (33/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 = 2.08e-01

AAGAAGGCGACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAGCGCTTCGGCTTCGTC  >  NC_002947/631904‑632198
                                                                                                                                                   |                                                                                                                                                   
aagaagGCGACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGcc                                                                                                                                                    >  8:8375/1‑149 (MQ=255)
 agaagGCGACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCg                                                                                                                                                   <  6:89353/149‑1 (MQ=255)
       cGACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGGCGCATGCCGCGTTCGGCGGGTACa                                                                                                                                             <  7:24151/149‑1 (MQ=17)
        gACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGGCGCATGCCGCGTTCGGCGGGTACaa                                                                                                                                            <  3:255165/149‑2 (MQ=17)
        gACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAg                                                                                                                                            >  4:66642/1‑149 (MQ=255)
          cgcgGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACGCGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCACCCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAggt                                                                                                                                          <  6:266517/149‑1 (MQ=255)
                 tAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGCCTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAg                                                                                                                                     >  2:197584/1‑147 (MQ=255)
                  agcagacgCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCa                                                                                                                                   <  4:25078/145‑1 (MQ=255)
                       aagCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGTATGCGCCGGGTACAGGTGGTAGCAGTTGGt                                                                                                                             >  4:142935/1‑149 (MQ=255)
                                      aTGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTCGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGt                                                                                                                             <  5:123743/134‑1 (MQ=255)
                                      aTGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTAGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGt                                                                                                                             >  6:123743/1‑134 (MQ=255)
                                        gTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCt                                                                                                            >  6:213279/1‑149 (MQ=255)
                                         tCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCtt                                                                                                           <  3:66642/149‑1 (MQ=255)
                                           gTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATTCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTg                                                                                                          <  4:45369/148‑1 (MQ=255)
                                           gTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTg                                                                                                          <  2:43894/148‑1 (MQ=255)
                                           gTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGa                                                                                                         >  4:260003/1‑149 (MQ=255)
                                                tCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGccc                                                                                                    >  8:130962/1‑149 (MQ=255)
                                                  aCCAGCAGGTTCTTGGGCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGccc                                                                                                    <  3:211163/147‑1 (MQ=255)
                                                   ccAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCgg                                                                                                 >  4:73599/1‑149 (MQ=255)
                                                    cagcagGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTAc                                                                                                                                              <  7:160258/103‑1 (MQ=255)
                                                    cagcagGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTAc                                                                                                                                              >  8:160258/1‑103 (MQ=255)
                                                    cagcagGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGc                                                                                                >  6:97750/1‑149 (MQ=255)
                                                      gcagGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGcc                                                                                               >  3:49518/1‑148 (MQ=255)
                                                              ttGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCAt                                                                                            >  7:151494/1‑143 (MQ=255)
                                                              ttGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCAt                                                                                            <  8:151494/143‑1 (MQ=255)
                                                              ttGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCAGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGAAGATGGTCCAAACACGACCAGACTTGATGCACAGACACATGCAGt                                                                                       >  4:157476/1‑148 (MQ=17)
                                                                ggTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGt                                                                                       <  7:8375/146‑1 (MQ=255)
                                                                 gTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTAGATGCCCCGGCCCATGCGGTAgg                                                                                    >  5:18909/1‑148 (MQ=255)
                                                                    tgctgATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGTCCCATGCg                                                                                         >  6:190175/1‑140 (MQ=255)
                                                                    tgctgATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCg                                                                                         <  5:190175/140‑1 (MQ=255)
                                                                      ctgATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGcccc                                                                                                   <  7:11555/128‑1 (MQ=255)
                                                                      ctgATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGcccc                                                                                                   >  8:11555/1‑128 (MQ=255)
                                                                      ctgATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCTCCGGCCCATGCGGTAGGCACGGt                                                                              >  6:270108/1‑149 (MQ=255)
                                                                            gTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATg                                                                         <  4:49518/148‑1 (MQ=255)
                                                                             tGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGAt                                                                          <  1:121490/146‑1 (MQ=255)
                                                                             tGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGAt                                                                          >  2:121490/1‑146 (MQ=255)
                                                                                 ccAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCcac                                                                                                                         >  7:6107/1‑95 (MQ=255)
                                                                                 ccAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCcac                                                                                                                         <  8:6107/95‑1 (MQ=255)
                                                                                 ccAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGAGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCgg                                                                   >  6:208542/1‑149 (MQ=255)
                                                                                  cAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCggg                                                                  <  5:213279/149‑1 (MQ=255)
                                                                                  cAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCggg                                                                  <  1:197584/149‑1 (MQ=255)
                                                                                     catcatCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTc                                                                >  1:1897/1‑148 (MQ=255)
                                                                                     catcatCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTc                                                                >  3:51494/1‑148 (MQ=255)
                                                                                     catcatCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTc                                                                <  4:51494/148‑1 (MQ=255)
                                                                                       tcatcTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCcac                                                             <  7:130962/149‑1 (MQ=255)
                                                                                          tcTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGc                                                                                          >  8:52855/1‑117 (MQ=255)
                                                                                          tcTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGc                                                                                          <  7:52855/117‑1 (MQ=255)
                                                                                               tgCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCgg                                                      <  3:73599/148‑1 (MQ=255)
                                                                                                       aCGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATCCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGCCGCCCAGGc                                             >  2:256986/1‑149 (MQ=255)
                                                                                                         gCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGcc                                            >  6:47749/1‑148 (MQ=255)
                                                                                                               cccccGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCg                                     <  5:47749/149‑1 (MQ=255)
                                                                                                                  ccGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGgtg                                  <  2:1897/149‑1 (MQ=255)
                                                                                                                    gACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTc                                >  6:242282/1‑149 (MQ=255)
                                                                                                                        acctGTACCCGGCGCATGCCGCGTTCGGCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTg                            >  5:114652/4‑149 (MQ=17)
                                                                                                                         cttGTACCCGCCTAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATTTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGtt                             <  3:157476/147‑1 (MQ=255)
                                                                                                                          ctGTACCCGGCGCATGCCGCGTTCGGCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTgg                           >  2:116932/2‑148 (MQ=17)
                                                                                                                                ccgccgAACGCGGCGTGCGCCGGGTACAGTTGGTAGCAGTTGGTCAACACAAGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCACGCCGAACTCGgt                                   >  7:227360/1‑134 (MQ=18)
                                                                                                                                ccgccgAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGc                    >  5:250066/1‑149 (MQ=255)
                                                                                                                                ccgacgTACGCGGCTTGCGCGGGGTAGAGGTGGTAGGAGTTGGTCCACACACGGCCAGCCTTGATGCCGCGGAGCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGgt                                   <  8:227360/134‑1 (MQ=17)
                                                                                                                                 cgccgAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGCCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGc                    >  4:225031/1‑148 (MQ=255)
                                                                                                                                  gccgAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGc                    >  6:191030/1‑147 (MQ=255)
                                                                                                                                   ccgAACGCGGTATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAg                 <  3:142935/149‑1 (MQ=255)
                                                                                                                                           ggCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAGCgcttcgg         >  5:55821/1‑149 (MQ=255)
                                                                                                                                             cATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAGCgcttcggct       <  5:97750/149‑1 (MQ=255)
                                                                                                                                                  gcCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAGCGCTTCGGCTTCGTc  >  2:218916/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGAAGGCGACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAGCGCTTCGGCTTCGTC  >  NC_002947/631904‑632198

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