Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,674,750:1 +G 100% intergenic (+112/+28) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,674,7501.G98.3% 193.0 / ‑3.5 58intergenic (+112/+28)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (29/28);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGC‑G‑CGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947/5674610‑5674890
                                                                                                                                              |                                                                                                                                             
gcCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAgc‑ggccccct                                                                                                                                       >  8:73929/1‑142 (MQ=255)
       gctgctGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCAGCAGAGAAATTCTTCATATCTAAAGGAGCGGCCATGTGTCGTGATAGG‑GCCGCAAAGC‑GTCCCCCTCGGttt                                                                                                                                 >  8:91103/1‑148 (MQ=37)
         tgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAgc‑g                                                                                                                                              <  4:201327/133‑1 (MQ=255)
         tgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAgc‑g                                                                                                                                              >  3:201327/1‑133 (MQ=255)
           ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGc                                                                                                                            <  2:317328/149‑1 (MQ=255)
              aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTgg                                                                                                                         >  6:394116/1‑149 (MQ=255)
                 cGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTgg                                                                                                                         >  5:46457/1‑146 (MQ=255)
                  gCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAg                                                                                                                      >  6:380444/1‑148 (MQ=255)
                      gggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAAt                                                                                                                  >  7:361596/1‑148 (MQ=255)
                      gggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAAt                                                                                                                  >  7:195555/1‑148 (MQ=255)
                       ggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTa                                                                                                                <  1:314100/149‑1 (MQ=255)
                         cTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGCGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTAc                                                                                                               <  1:269875/149‑1 (MQ=37)
                         cTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACt                                                                                                              >  1:361798/1‑149 (MQ=255)
                           gATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTc                                                                                                            >  1:254022/1‑149 (MQ=255)
                            aTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTc                                                                                                            <  1:38432/148‑1 (MQ=255)
                            aTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCCGTTTTCAGCTGGCAGAAATTACTTCt                                                                                                           >  4:13314/1‑149 (MQ=255)
                                aaGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                                                                                                       >  5:185870/1‑149 (MQ=255)
                                  gTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTca                                                                                                     <  2:254022/149‑1 (MQ=255)
                                        aGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                               <  1:52815/149‑1 (MQ=255)
                                         gggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                              <  8:361596/149‑1 (MQ=255)
                                            ggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGac                                                                                           <  7:64276/149‑1 (MQ=255)
                                               tgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt                                                                                        >  5:151614/1‑149 (MQ=255)
                                                       aGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtcca                                                                                <  3:168265/149‑1 (MQ=255)
                                                         cccTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTTAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccagg                                                                              >  2:135106/1‑149 (MQ=255)
                                                                aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                        <  3:13314/148‑1 (MQ=255)
                                                                    cacaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTg                                                                   <  6:46457/149‑1 (MQ=255)
                                                                     acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                  <  7:91103/149‑1 (MQ=255)
                                                                       aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCg                                                                 >  1:68305/1‑148 (MQ=255)
                                                                          gCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                              >  4:321234/1‑148 (MQ=255)
                                                                            cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGa                                                                <  6:330031/144‑1 (MQ=255)
                                                                            cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGa                                                                >  5:330031/1‑144 (MQ=255)
                                                                             gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAg                                                           >  6:213464/1‑148 (MQ=255)
                                                                                  tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCg                                                     >  5:157546/1‑149 (MQ=255)
                                                                                    cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                    <  2:361798/148‑1 (MQ=255)
                                                                                       agaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                             >  5:136288/1‑136 (MQ=255)
                                                                                       agaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                             <  6:136288/136‑1 (MQ=255)
                                                                                       agaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCt                                                >  3:341264/1‑149 (MQ=255)
                                                                                        gaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat                                                                                     >  4:78249/1‑111 (MQ=255)
                                                                                        gaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat                                                                                     <  3:78249/111‑1 (MQ=255)
                                                                                         aaTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCttt                                              <  7:73929/149‑1 (MQ=255)
                                                                                             tCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAt                                           <  5:380444/148‑1 (MQ=255)
                                                                                                 gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGc                                         <  3:321234/146‑1 (MQ=255)
                                                                                                 gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCg                                       >  3:344129/1‑148 (MQ=255)
                                                                                                   tctcTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                     <  3:345557/148‑1 (MQ=255)
                                                                                                    ctctACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGt                                   >  2:77960/1‑149 (MQ=255)
                                                                                                     tctACAGGAGCGGCCTTGTCTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGACTGCGAAGTAGCACCCGGGCCTTTGATTCCGAAGTc                                  >  4:160799/1‑149 (MQ=255)
                                                                                                         cAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgtt                               <  1:77960/148‑1 (MQ=255)
                                                                                                             aGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTTATATCCAGGTCCAGGGTCTGCTAAGTAGCACCCGGGCCTTTGATGCCGAATTCgttcaggt                          <  6:185870/149‑1 (MQ=39)
                                                                                                              gCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggtt                         <  5:213464/149‑1 (MQ=255)
                                                                                                                  ccTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccag                                                                               >  4:18827/1‑91 (MQ=255)
                                                                                                                  ccTTGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccag                                                                               <  3:18827/91‑1 (MQ=255)
                                                                                                                    ttGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgt                   >  7:213826/1‑149 (MQ=255)
                                                                                                                    ttGTGTCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgt                   >  7:223142/1‑149 (MQ=255)
                                                                                                                      gtgtCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtgg                 >  5:327403/1‑149 (MQ=255)
                                                                                                                         tCGTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg               <  6:327403/148‑1 (MQ=255)
                                                                                                                           gTGATAGG‑GCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg            >  3:150137/1‑149 (MQ=255)
                                                                                                                                    gCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa    <  5:394116/149‑1 (MQ=255)
                                                                                                                                    gCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa    <  8:195555/149‑1 (MQ=255)
                                                                                                                                    gCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa    <  2:38798/149‑1 (MQ=255)
                                                                                                                                     ccGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAg   >  3:228900/1‑149 (MQ=255)
                                                                                                                                      cGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGCCACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc  >  7:180125/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGG‑GCCGCAAAGC‑G‑CGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947/5674610‑5674890

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