Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,655,636 C→A 10.2% intergenic (‑91/+61) epd ← / ← tktA D‑erythrose 4‑phosphate dehydrogenase/transketolase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,655,6360CA10.2% 143.5 / 7.9 58intergenic (‑91/+61)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base A (3/3);  total (29/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.26e-01

CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCC  >  NC_002947/5655488‑5655781
                                                                                                                                                     |                                                                                                                                                   
cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGc                                                                                                                                                     <  2:310012/149‑1 (MQ=255)
  caGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCaa                                                                                                                                                   <  6:177321/149‑1 (MQ=255)
   aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                                  >  8:254163/1‑149 (MQ=255)
         cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                            <  5:298219/149‑1 (MQ=255)
         cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGCGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                            <  8:144322/149‑1 (MQ=255)
              cGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGa                                                                                                                                        <  3:138237/148‑1 (MQ=255)
                  cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCggg                                                                                                                                   <  6:128199/149‑1 (MQ=255)
                             aGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGTAATTTGCGGGGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGGGGGAGCGGGCTTGCCcgcg                                                                                                                         <  8:21141/148‑1 (MQ=255)
                              gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCgg                                                                                                                                    >  4:315385/1‑136 (MQ=37)
                              gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCgg                                                                                                                                    <  3:315385/136‑1 (MQ=37)
                                 cAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGAGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                     >  2:300917/1‑148 (MQ=255)
                                    cTTTGTAGGGACGGGGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGGGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGTGCGGGCTTGCCCGCGAATGAAgg                                                                                                                <  3:214992/149‑1 (MQ=255)
                                     tttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGc                                                                                                               >  3:124452/1‑149 (MQ=255)
                                      ttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGc                                                                                                               <  4:124452/148‑1 (MQ=255)
                                         tAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGCCAAGCCCGCTCCCACAGGGCTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGaca                                                                                                           <  3:286059/149‑1 (MQ=255)
                                            ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                                      <  8:282261/136‑1 (MQ=255)
                                            ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                                      >  7:282261/1‑136 (MQ=255)
                                                     ggCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCACCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCCACACCGCCCGCaca                                                                                                >  1:199074/1‑146 (MQ=255)
                                                      gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATc                                                                                              <  2:238016/149‑1 (MQ=255)
                                                         ttCGGGGCTCACGAGACATTTCACTTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCttt                                                                                           <  1:41028/149‑1 (MQ=255)
                                                         ttCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCtt                                                                                            >  5:279229/1‑148 (MQ=255)
                                                             gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCGGTGGGAGCGGGCGTGCCCGCGAATGAAGGCGCCCCCGCCCGCAGCTCTTTaaa                                                                                        >  7:265299/1‑148 (MQ=38)
                                                              ggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAAt                                                                                       >  1:326517/1‑148 (MQ=255)
                                                              ggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGAGCTTTTAAt                                                                                       >  7:304581/1‑148 (MQ=255)
                                                               gCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCa                                                                                     <  5:351928/149‑1 (MQ=255)
                                                                  cACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTc                                                                                  <  3:160827/149‑1 (MQ=255)
                                                                     gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACCGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGGCCCCGCCCGCAGCTCGTTAAAtagatcg                                                                                 >  4:301847/1‑141 (MQ=37)
                                                                     gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTaaa                                                                                        <  3:301847/140‑1 (MQ=255)
                                                                        aCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCATGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                             <  1:300917/148‑1 (MQ=255)
                                                                        aCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcag                                                                            <  6:343073/149‑1 (MQ=255)
                                                                                    tcgtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCgg                                                                <  1:293880/148‑1 (MQ=255)
                                                                                       tcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAATCCCCCTACCACAAGGTTTAATTCAACCACTGTAGGAGCGGGCTTTGCCCCGAATGAAGGCGACACCGCCCCCCGATCTTTAAATAATTATTACACAAGCTCATAGGCAg                                                             >  8:306973/1‑148 (MQ=255)
                                                                                            cTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACCCCGCCCGCAGCTCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTAcc                                                         >  5:167987/1‑147 (MQ=255)
                                                                                               cGGAATTTGCAGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                      <  8:251770/147‑1 (MQ=255)
                                                                                                     ttGCGGTGTGGCTTTCGCGTGAAAGCCCGATCCCACAGGTTTTCATGAAACCTCTTTGTGAGCGGGCTTGCCCGCGAATGAAGGCTACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtc                                               <  7:306973/148‑1 (MQ=38)
                                                                                                       gCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtctc                                             >  4:265181/1‑148 (MQ=255)
                                                                                                       gCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCGCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATAAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGttggc                                             >  8:18773/1‑145 (MQ=255)
                                                                                                          gtgtGGCATTCGCGGGAAAGCCCGCTCCCACATGGTTTCATGCAACCTCTGTGGTAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATGAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                         <  2:199074/149‑1 (MQ=255)
                                                                                                           tgtgGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCtt                                                                                            >  1:300669/1‑98 (MQ=255)
                                                                                                           tgtgGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCtt                                                                                            <  2:300669/98‑1 (MQ=255)
                                                                                                                     gcgGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                              >  5:312172/1‑149 (MQ=255)
                                                                                                                     gcgGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                              >  7:249436/1‑149 (MQ=255)
                                                                                                                     gcgGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                               >  4:232080/1‑148 (MQ=255)
                                                                                                                       gggCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAAct                            >  8:73123/1‑149 (MQ=255)
                                                                                                                         ccatgCCCCCTCCCACAGGGTTTATTGCACCCTCTGTGGTCGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAActct                          <  8:304581/145‑1 (MQ=255)
                                                                                                                          cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                         <  6:279229/149‑1 (MQ=255)
                                                                                                                          cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                         <  2:326517/149‑1 (MQ=255)
                                                                                                                           aaGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                        >  1:150586/1‑149 (MQ=255)
                                                                                                                           aaGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                        <  3:232080/149‑1 (MQ=255)
                                                                                                                                cGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                     >  6:288989/1‑148 (MQ=17)
                                                                                                                                  cTCCCACAGGGTTTCATGCACCCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCTGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt                 >  8:4782/1‑149 (MQ=255)
                                                                                                                                      cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGcc                                 <  5:136318/123‑1 (MQ=37)
                                                                                                                                      cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGcc                                 >  6:136318/7‑129 (MQ=37)
                                                                                                                                      cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCAAGCAGCTCATCTAGCAGTAAAGAGGATGTTATCAAGCGTGAAGACAAACTCTTCGAACaaagaa              >  1:303882/1‑145 (MQ=255)
                                                                                                                                       acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg             >  3:338832/1‑148 (MQ=255)
                                                                                                                                       acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg             >  8:74202/1‑148 (MQ=255)
                                                                                                                                                tCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc     <  7:254163/148‑1 (MQ=255)
                                                                                                                                                 cATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgcgc   >  3:341706/1‑149 (MQ=255)
                                                                                                                                                  aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGcc  >  8:111658/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCC  >  NC_002947/5655488‑5655781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: