Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,655,631 T→G 9.2% intergenic (‑86/+66) 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,6310TG 9.2% 180.4 / 8.4 64intergenic (‑86/+66)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base G (3/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAA  >  NC_002947/5655483‑5655772
                                                                                                                                                         |                                                                                                                                             
cAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGttt                                                                                                                                               <  7:108256/149‑1 (MQ=255)
 aGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTc                                                                                                                                              <  7:255847/149‑1 (MQ=255)
     cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCAGCTCCCACCGGGGTTCCTGc                                                                                                                                          >  6:240330/1‑149 (MQ=255)
        aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCTGGGCTCACTAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGAGGTGTGT‑‑‑‑CATTAGCGGGCAAGCCAGCTCCCACATGGTTTCATTCAAc                                                                                                                                       >  7:124167/1‑149 (MQ=255)
          cGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACct                                                                                                                                     <  5:66307/149‑1 (MQ=255)
            gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctct                                                                                                                                   >  3:253668/1‑149 (MQ=255)
              agATACGGACATATACGTTGAGTGGAACTTTTTTGGGACGGTTGGGGATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAATCCCGCTCCCTCAGTGTTTCATGCAACATCtgt                                                                                                                                 <  4:34225/148‑1 (MQ=255)
                    ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTGATGCAACCTCTGTGGGAGc                                                                                                                           >  8:148835/1‑149 (MQ=255)
                     gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCg                                                                                                                          <  1:74538/149‑1 (MQ=255)
                      ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                             <  4:34647/149‑1 (MQ=21)
                        aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                       <  2:5325/133‑1 (MQ=255)
                        aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                       >  1:5325/1‑133 (MQ=255)
                           tcagtgtaGTGCAACTTTGTAGGGACGGTGGGGCTTTCGTGGTTCACGAGACAGTTCAATTGTCGTCAGCATGAGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCAGAGGTTTTCATGCAACCTCTGTGGGAGCGGGCtt                                                                                                                     <  8:124167/141‑1 (MQ=38)
                           aCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCtt                                                                                                                     <  7:277326/148‑1 (MQ=255)
                            ccGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                 <  8:109758/135‑1 (MQ=255)
                            ccGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                 >  7:109758/1‑135 (MQ=255)
                                    tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATg                                                                                                                                           >  6:216737/1‑118 (MQ=255)
                                    tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATg                                                                                                                                           <  5:216737/118‑1 (MQ=255)
                                    tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCTGGCTTGCCCGCGAAt                                                                                                           >  5:334431/1‑149 (MQ=255)
                                       aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                          >  7:40200/1‑147 (MQ=255)
                                            tGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCaa                                                                                                                                        >  4:83292/1‑112 (MQ=255)
                                            tGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCaa                                                                                                                                        <  3:83292/112‑1 (MQ=255)
                                            tGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCCGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGGTTCCACAGGGTTTCATGCAACCTCTGTGGGATCGGGCTTGCCCGCGAATGAAGGc                                                                                                     <  5:241277/147‑1 (MQ=255)
                                              tAGGGACGGGGGGGCATTCGGGGCGCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGac                                                                                                  <  7:204274/148‑1 (MQ=255)
                                                    cGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGCTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgccc                                                                                           <  5:240330/149‑1 (MQ=255)
                                                           gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCCCGCCCCCACGGGTTGATGGAACCACTGTGGGAGCCGGGGTGGCCGCGAATTAAgg                                                                                                      >  1:339691/1‑131 (MQ=255)
                                                                 ggggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCCCACCCACAGGGTTTAATGCAACCTCTGTGGGAGCGGGCGTGCCCGCGAATGAAGGCGACACCGCCCGCAGATGGTTaa                                                                               >  8:41029/1‑148 (MQ=255)
                                                                 ggggCTAACGAGACATTTCAATTGTCGTCTGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGGGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGGCACCGCCCGCAGATCTTTaaa                                                                              <  8:40200/149‑1 (MQ=38)
                                                                  gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAAt                                                                             <  7:41029/149‑1 (MQ=255)
                                                                        aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCACAGTGGGAGCGGGCTTTCCCGCGAAAGAATGCGCCACCGCCCGCAGATCGTTAAATAAGTCt                                                                       >  8:160897/1‑149 (MQ=38)
                                                                        aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAACCCCCCACCCACAGGGGTTCATGGAACCTCTGTGGGACCCGGGATTCCCGCCAATGTAGGCGACACCGCCCGCAGATCTTTAAATCattgt                                                                       >  2:214148/1‑145 (MQ=255)
                                                                          gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                     <  1:214148/149‑1 (MQ=255)
                                                                           agaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCCCAGGGTTTCATACCCCCGCGGTGGGAGCGGGCTTGCCCGCGGATGACGGCGCCGCCGCCCGCGGAGCGGTAAAGGAGTCtggg                                                                    >  4:102550/1‑146 (MQ=255)
                                                                                 ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCagcagc                                                              >  5:333718/1‑149 (MQ=255)
                                                                                 ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCagcagc                                                              >  6:331782/1‑149 (MQ=255)
                                                                                 ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAAACCCCACCCACAGGGGTTTAAGGAAACGCTGTGGGAGCGGGGGTGGCCCCGGAATAAGGCGAACCCGCCCGCAGATCTTTTAATAATTCTTGcagcag                                                               >  6:24945/1‑148 (MQ=255)
                                                                                         cgttctgCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCTAGCCCGCTCCCACAGGGTCTCATGGAAGCTCTGTGTGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCgt                                                       <  7:217104/142‑2 (MQ=38)
                                                                                           gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                       <  2:172512/82‑1 (MQ=255)
                                                                                           gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                     <  7:160897/148‑1 (MQ=255)
                                                                                           gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAg                                                    >  6:342809/1‑149 (MQ=255)
                                                                                           gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACACGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGGCCGCGAATGAAGGCGACACCGCCCGCCGATCGTTAAATAAGTCTTCCAGCAGCTCGTCGGCAg                                                    >  6:1302/1‑149 (MQ=38)
                                                                                              aGCCTGCGGAATTTGCGGTGGGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTa                                                  <  3:6127/148‑1 (MQ=255)
                                                                                                  tGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGCTCTTTAAATCAGTCTTCCAGCCGCTCTTCGGCAGTACCGAg                                             >  5:79002/1‑149 (MQ=255)
                                                                                                    cGGAATTTGCGGGGTGG‑‑‑‑CATTCGCGGTCAAGCCCGCTCCCACAGGGGTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGat                                             <  6:334431/147‑2 (MQ=255)
                                                                                                         tttGCGGGGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtc                                      <  5:24945/149‑1 (MQ=255)
                                                                                                               gtgtGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                <  3:169023/149‑1 (MQ=255)
                                                                                                                  tgG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGt                              <  5:331782/148‑1 (MQ=255)
                                                                                                                          aTTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                >  6:346575/1‑142 (MQ=255)
                                                                                                                          aTTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                <  5:346575/142‑1 (MQ=255)
                                                                                                                              gcgGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                     <  7:148835/149‑1 (MQ=255)
                                                                                                                              gcgGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                     <  4:253668/149‑1 (MQ=255)
                                                                                                                                      gCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCg              >  3:104466/1‑148 (MQ=255)
                                                                                                                                       cccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGa             >  4:52096/1‑148 (MQ=255)
                                                                                                                                        ccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa            >  5:295871/1‑148 (MQ=255)
                                                                                                                                           cTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt        >  3:219671/1‑149 (MQ=255)
                                                                                                                                            tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAAT‑‑‑‑GCCACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg    >  2:226121/1‑148 (MQ=21)
                                                                                                                                            tCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTATTCGAACAATg       >  5:317204/1‑149 (MQ=255)
                                                                                                                                             cccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGc      >  3:138222/1‑149 (MQ=38)
                                                                                                                                              ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGCAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                <  2:122932/138‑1 (MQ=37)
                                                                                                                                              ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc     <  5:342809/149‑1 (MQ=255)
                                                                                                                                                acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTcc                                   <  2:294963/117‑1 (MQ=255)
                                                                                                                                                acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTcc                                   >  1:294963/1‑117 (MQ=255)
                                                                                                                                                 caGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAACAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa  >  4:333608/1‑149 (MQ=255)
                                                                                                                                                  aGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                >  7:31472/1‑118 (MQ=255)
                                                                                                                                                  aGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                <  8:31472/118‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                             
CAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAA  >  NC_002947/5655483‑5655772

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: