Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 632,043 G→C 6.4% A466A (GCC→GCG)  aldB‑I ← aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947632,0430GC 6.4% 196.8 / 3.9 62A466A (GCC→GCG) aldB‑Ialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base C (2/2);  total (33/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 = 2.78e-01

AAGAAGGCGACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAGCGCTTCGG  >  NC_002947/631904‑632191
                                                                                                                                           |                                                                                                                                                    
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                                                                                                                                       >  8:160258/1‑103 (MQ=255)
                                                    cagcagGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTAc                                                                                                                                       <  7:160258/103‑1 (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                                                                   >  2:121490/1‑146 (MQ=255)
                                                                             tGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGAt                                                                   <  1:121490/146‑1 (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)
                                                                                                                                           |                                                                                                                                                    
AAGAAGGCGACGCGGTCTAGAAGAAGCCCAGCGGATTGATGTCGTAGCTCACCAGCAGGTTCTTGGTCTGCTGATAGTGATCCAGCATCATCTTGTGCGTCTCACGCCCAACCCCCGACTTCTTGTACCCGCCGAACGCGGCATGCGCCGGGTACAGGTGGTAGCAGTTGGTCCACACACGGCCAGCCTTGATGCCCCGGCCCATGCGGTAGGCACGGTTGATGTCGCGGGTCCACACGCCGGCGCCCAGGCCGAACTCGGTGTCGTTGGCAATTGCCAGCGCTTCGG  >  NC_002947/631904‑632191

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