Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT29.6% 103.9 / 14.7 54S64F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (3/35);  new base T (16/0);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 4.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTG  >  NC_000913/4604114‑4604546
                                                                                                                                                                                                                                            |                                                                                                                                                                                                     
gTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                      <  1:252596/238‑1 (MQ=255)
                     cTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                                                                                <  1:429388/239‑1 (MQ=255)
                                     tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACTAGTCAGGAGTCATTTCAGGCTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCTTGCGAAGGGGGGGCAGAGAggggct                                                                                                                                                                         >  2:91192/1‑227 (MQ=255)
                                             aaTGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                                                                          <  1:759981/205‑1 (MQ=255)
                                                       aGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGAGAGCGGCGAGAGGGGTTATTTTGCCCGTCTGGACTGGGGATTaa                                                                                                                                             >  2:595306/1‑239 (MQ=255)
                                                                                      tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAc                                                                                                                                 <  1:678106/220‑1 (MQ=255)
                                                                                       aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGAAGTGAGTGCCGAGAGGGGTTATCTTGCCTGTCTGAACTGGCGGTTAAGTTTTCCCGGACACCGCGTTTTGTttttttgg                                                                                                             >  2:161346/1‑235 (MQ=255)
                                                                                          aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGGGGGTAGCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                                                                           <  1:91192/238‑1 (MQ=255)
                                                                                          aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                                                                              <  1:404271/235‑1 (MQ=255)
                                                                                           aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATGCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATCCCTATACCAGTCAGGAGTGTTTTCAGGAGGCCATGTCGGGGGTCTAGTTTGCGGAGATGATTTTTTCTTTTACTGACATGAGAAGTGAGCGCAGAGAGGGAGTATCTGGCCTGACTATACTCGCGATTAAATTCCCGCGTACCCGGCGTGTAGTTAGCGCggagga                                                                                                         <  2:399066/239‑1 (MQ=255)
                                                                                                cGTCATGAGTAGTATCGGGATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGCATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTCTTTTTCTTTTTTTGCCCTTCTAAAGGAGCGCATAGGGGGGGTATCTTGGCTTGCCGAACTGGTGGATAAGGTTCCCCGTGACCGGCggcta                                                                                                                       >  2:678106/1‑216 (MQ=255)
                                                                                                     tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                                                                               <  1:407194/239‑1 (MQ=255)
                                                                                                                ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAt                                                                                    <  1:161346/239‑1 (MQ=255)
                                                                                                                   aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCCTGAGAAGGTAGCGCTGGGCGGGGTTGTTTTGCCTGACAGAACTGGCGATTAAGTTTCGCCGGACCCGCCGTTTAGTTTTTGCGGAGAGTGATATTGAAGCTCTGGTGATTgg                                                                                 >  2:334002/1‑239 (MQ=255)
                                                                                                                      gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGAGGGCCGGCGGGGGTTTTCTTTTCTTACTTAAATGGCGATTAAAGTTTCGCGTACCCCGGGGTTTGTTATTGCGGGTGTTGATATTTAAGATCTGGTGGTTTGttt                                                                              >  2:365828/1‑238 (MQ=255)
                                                                                                                         agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                                                                                       <  1:709039/132‑1 (MQ=255)
                                                                                                                                   gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGTGAgg                                                                                                                                                                             >  2:709039/1‑132 (MQ=255)
                                                                                                                                     aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                                                                   <  1:2694/235‑1 (MQ=255)
                                                                                                                                     aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                                                                   <  1:595306/235‑1 (MQ=255)
                                                                                                                                             ggttAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACg                                                       <  2:840273/236‑1 (MQ=255)
                                                                                                                                             gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGGGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGAAGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc                                                        <  1:121083/238‑1 (MQ=255)
                                                                                                                                             gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTGCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACg                                                       <  1:30281/239‑1 (MQ=255)
                                                                                                                                                    ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                                                                                   <  1:236691/109‑1 (MQ=255)
                                                                                                                                                    ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                                                                                            <  1:193243/116‑1 (MQ=255)
                                                                                                                                                    ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt                                                                                                                          <  1:172572/165‑1 (MQ=255)
                                                                                                                                                    ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                                                                        <  1:164666/215‑1 (MQ=255)
                                                                                                                                                      ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGtt                                                                                                                         <  1:185850/164‑1 (MQ=255)
                                                                                                                                                      ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                                                                                        <  1:405074/181‑1 (MQ=255)
                                                                                                                                                      ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAATTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc                                                        <  1:401026/229‑1 (MQ=255)
                                                                                                                                                         taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTa                                           <  1:111853/239‑1 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTAGTTATTTCGGATGATTATATTGAAGCTCGGCTGGTTTGTTTATTGTTGCCCTTACCCCTGGGCGGGGTATTTAGGAAAACg                                      >  1:13944/1‑239 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGGAGCGAGGGCAGGGAGGGGTTTTTTTGCCTTACTGAACTTGCGGTTTAGTTTCCGCGTACCCGGCGTTTTGTTATTGTGGGTGATGGTATTTTAGTTCGGCCGATTGGTTCATTGTCTCCATCACCg                                                              >  2:164666/1‑215 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt                                                            <  1:371419/217‑1 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATGAGAAGTGGGGGGGGAGGGGGGGTATCTTGGCTTACTGAAATGGGGGTTAAGTTTTCCCGGGCCCGGCGGTTTGTTTTTg                                                                                                                >  2:172572/1‑165 (MQ=255)
                                                                                                                                                                cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGAGGGGGTATCTTGTCTGTCTGAAATGGCGGTTAAAGTTTCCCGGGACCCGGGGTTTGTTATTGTGGGTGATGATATTTAAg                                                                                              >  2:405074/1‑181 (MQ=255)
                                                                                                                                                                cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                                                                          >  1:753621/1‑90 (MQ=255)
                                                                                                                                                                cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                                                                          <  2:753621/90‑1 (MQ=255)
                                                                                                                                                                   aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                      <  1:551907/75‑1 (MQ=255)
                                                                                                                                                                   aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                      <  1:691050/75‑1 (MQ=255)
                                                                                                                                                                   aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                                                                                  <  1:241761/111‑1 (MQ=255)
                                                                                                                                                                   aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa                                                  <  1:96998/222‑1 (MQ=255)
                                                                                                                                                                    ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                                                                           <  1:200641/164‑1 (MQ=255)
                                                                                                                                                                    ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                                                                                  <  1:687726/173‑1 (MQ=255)
                                                                                                                                                                     cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                                                                                            <  1:233756/162‑1 (MQ=255)
                                                                                                                                                                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATTTTTCCTGTCTGAAc                                                                                                                                                        >  2:241761/1‑111 (MQ=255)
                                                                                                                                                                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTGCCTGAGaa                                                                                                                                                                                            >  2:691050/1‑75 (MQ=255)
                                                                                                                                                                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATGAGaa                                                                                                                                                                                            >  2:551907/1‑75 (MQ=255)
                                                                                                                                                                              tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGGGCGCCGGGAGGGATTATCTTGTCTGACTGAACTGGCGGTTAAGCTTCCGCGGACCCGGCGGTTTGTTATTGCGGATGATGATATTTAAGGTCGGCCGGTTGGTTTATTGTTGCCCTCACCCCTGGGGGGGGTGTTTAGTAAAGCGGCAAccgcgggggttt                      >  2:515565/1‑228 (MQ=255)
                                                                                                                                                                                  gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                                                                          >  1:222162/1‑72 (MQ=255)
                                                                                                                                                                                  gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                                                                          <  2:222162/72‑1 (MQ=255)
                                                                                                                                                                                  gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgc                                                                                                                                     <  1:740459/124‑1 (MQ=255)
                                                                                                                                                                                   aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcg                                                                                                                                      <  1:405079/122‑1 (MQ=255)
                                                                                                                                                                                       aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAgg             <  1:165856/239‑1 (MQ=255)
                                                                                                                                                                                            aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGTGAAGTGAGCGCAGAGAGGGGTTATATTTCCTGTCTGAACTTGCGGTTTAGTTTTCGCGGACCCGGCg                                                                                                                           >  2:740459/1‑124 (MQ=255)
                                                                                                                                                                                             ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTAAGTTTTCGCGTAACCCGc                                                                                                                            >  2:405079/1‑122 (MQ=255)
                                                                                                                                                                                                   cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTg  <  1:365828/238‑1 (MQ=255)
                                                                                                                                                                                                   cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTg  <  1:515565/238‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                     
GTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTG  >  NC_000913/4604114‑4604546

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