Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT31.4% 59.3 / 10.3 36S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (3/21);  new base T (10/1);  total (14/22)
Fisher's exact test for biased strand distribution p-value = 1.53e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAG  >  NC_000913/4604095‑4604591
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
gggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAAGGCGTCATGAGTAGTATCGATATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCATGATGCCATGTCTCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                                                                                                                                                                                         <  1:261800/250‑1 (MQ=255)
     aaaTAGCGCGCAGGGTGCAAGAGCACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTACGAAGAGAAATGTCGCACTCATAGAAAACTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGCTCTCGTTTGCGGCGGTCTCtttttttttttt                                                                                                                                                                                                                                                    >  2:154555/1‑250 (MQ=255)
                         gagaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGTAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                                                                                                                                    <  1:143860/247‑1 (MQ=255)
                          agaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                                                                                                                               <  2:353786/251‑1 (MQ=255)
                                             gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                                                                                                                            <  1:209329/251‑1 (MQ=255)
                                                 aTGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                                                                                                                         <  1:343247/250‑1 (MQ=255)
                                                           cGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAt                                                                                                                                                                                              <  1:467831/251‑1 (MQ=255)
                                                              aGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGAATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTAATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGTCATGGGAAGTGAGCGTAGAGAGGGCATATCTTGTATGGCTGAAATGGCgggtg                                                                                                                                                                                            >  2:83388/1‑246 (MQ=255)
                                                                                aTTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                                                                         <  2:63375/251‑1 (MQ=255)
                                                                                  tAAGAAGAGAAATGTCGCACGCATAGAAAAATGCGTCATGAGTAGTAGCGGTATTGAGAGTCTATTCAGAAAGTTTGCGGGTAACCCTGATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTACCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt                                                                                                                                                                        <  2:195871/250‑1 (MQ=255)
                                                                                   aagaagAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGGCATGAGGAGTGAGGGCAGGGGGGGGTTATTTTTCCTGACTTAACTGGCGGTTAAGTTTTCTCGTACCCGGCggtt                                                                                                                                                                      >  2:92847/1‑248 (MQ=255)
                                                                                   aagaagAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt                                                                                                                                                                      <  2:100339/251‑1 (MQ=255)
                                                                                     gaagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt                                                                                                                                                                      >  1:363156/1‑249 (MQ=255)
                                                                                      aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGATCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                                                                                                                   <  2:90084/251‑1 (MQ=255)
                                                                                                 cGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                                                                                                                        <  1:65970/251‑1 (MQ=255)
                                                                                                    aCTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGTGAAGGGAGCGCAGGGGGGGGTTATCTTGCCTGTCTTGACTGGCGGTTACATTTCCCCGTAccccgggttttgttgttttgggtgatg                                                                                                                                                     >  2:156887/1‑227 (MQ=255)
                                                                                                       cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGAGCAGAGAGGGATTATCTTGCCTGACTGAAATGGCGGTTAAGGTTCCGCGGACCCGGCGTTTTGTTTTTGCGgatgatg                                                                                                                                                     >  1:262940/1‑248 (MQ=255)
                                                                                                          aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCAATTTTTTCGTTTTTTTCCATCAGGAGAGTGCGCGGCGAGGGCTTATCTTGCCTGGCTTGACTGGCGATTAAGTTTCCCCGTTCACGCCGATTAGTTATTTCgggtggtggtgatg                                                                                                                                               >  1:192030/1‑239 (MQ=255)
                                                                                                                 tGCGTCATGAGTAGTATCGGTATTGGGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGGGAGCGCCGAGAGGGGTTATCTTTCCTTACTGAACTGGCGGTTAAATTTCCGCGTACCCGGCGGTTTGTTATTGCGGGTGATGATATTTAAGCTCg                                                                                                                                        >  1:170912/1‑251 (MQ=255)
                                                                                                                    gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTGCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                                                                                                                                     <  1:363796/251‑1 (MQ=255)
                                                                                                                      cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTTCCCTGAGGAATGAGGGCAGAGAGGGGTTATTTTTCCTGCCTGAAATGGCGATTCAGATTTCGCGTGCCCCGCGGTTTGGTATTGCTGATGATGATATTTAAgaccagctg                                                                                                                                    >  1:18017/1‑242 (MQ=255)
                                                                                                                               taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGGGGAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGGGATTAAGGTTTCGCGTAACCCGCGTTTTATTATTGCCGATGCTGGTATTTAAACTCGGCCGATTgggttat                                                                                                                          >  1:205373/1‑246 (MQ=255)
                                                                                                                                                    cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTAGTCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAGCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCAGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc                                                                                                      <  2:18017/250‑1 (MQ=255)
                                                                                                                                                        aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGTTTTTATTTTTCTGACATGAGAAGTGAGCGCAGAGAGGGATTATCTTGTCTGACTGGACTGGCGATTAAGTTTTCCCGTACCCCGCGTTTTGTTATTTCGGGAGTTGATATTTAAGTTCGGCTTATTGGTTCATTTTCGcccttccccctggccggggg                                                                                                 >  1:428254/1‑233 (MQ=255)
                                                                                                                                                        aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtgt                                                                                                 <  1:291281/251‑1 (MQ=255)
                                                                                                                                                                           ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACg                                                                              <  1:45653/251‑1 (MQ=255)
                                                                                                                                                                                 tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCCTGAGAAGTGGGCGCCGAGGGGGATTATTTTTCCTGTCTGAACTGGCGGATAAGTTTTCGCGTACCCGGCGGTTTGTTTTTTCGGGGGGTGATATTTAAGCCCGGCTGGTTGGGTCATTTTTGCCATCCCCGCTGGGAGGTGTATTAAGTAAAGCATCCACGccgg                                                                          >  2:60184/1‑246 (MQ=255)
                                                                                                                                                                                  ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTGTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAGTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGGTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGga                                                                        <  1:154555/249‑1 (MQ=255)
                                                                                                                                                                                       ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGAttttt                                                                 <  1:194625/251‑1 (MQ=255)
                                                                                                                                                                                                  cAGGAGTCATTTCAGGATGCGCTGTCGCGGATCTCGTTTGCGGAGGTCATTTTTTTTTTTTCGGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACGGAAATGGATATTAAGTTGCCGCGTCCCCGGCGTTTAGTTATTGCGGAGGATGATATGGAAGGTCGGCTGATTGGTTCATTGTCGCTATCGCCTCTGGACGGTGTATTAAGTAGAGCGTCAACGC‑TTGAGATTTTTCATCAGGCAct                                                      <  2:363755/251‑1 (MQ=255)
                                                                                                                                                                                                     gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTTCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCtt                                                   <  2:134766/251‑1 (MQ=255)
                                                                                                                                                                                                       gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCtttt                                                 <  1:320901/251‑1 (MQ=255)
                                                                                                                                                                                                         cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTGCATGAGGAGGGAGCGCAGAGGGGGGTTTTCTTTTCTGACTGAACTGGCGGTTTAGTTTTCGCGGACCCGGCGGTTTGTTATTGCGGGTGATGGTATTTAAACCCGGCGGATTTGGTCATTGTTGCCCTTACCCCTGGGCGGGGGATTAAGAAAAACGACAACGCGTGG‑GAtttttttttcggaacactttttt                                               >  2:6996/1‑233 (MQ=255)
                                                                                                                                                                                                         cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCttttt                                                <  1:73995/250‑1 (MQ=255)
                                                                                                                                                                                                              cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCTCCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                                          <  2:325421/251‑1 (MQ=255)
                                                                                                                                                                                                               aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCCTGAGAAGTGAGCGGAGGGGGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGGTTTCGCGGGCCCGGCGGTTAGTTATTTCGGGTGATGAAATTGAAGCTCGGCTGGTTGGTTCTTTGTTGCCATCCCCGCTGGGCGGGGGATTAAGTAAAAGGTCCACGC‑CGGGGGTTTTTCATCAGGGACTCttttttttattt                                         >  2:279974/1‑244 (MQ=255)
                                                                                                                                                                                                                                agctcGTTTGCGGCGGTCATGTTTTCTGTTTCCGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTGGCATGACTGAACTTGCGATTAAGTTACCGCGTTCCCGGCGTTTAGTTATTGCGGATGATGGTAGTGAAGCTCGGCTGATTGGTTCATTGTCGCCATTACCGCTGGACGGTGTATTAAGTAAAGCGACAACGC‑TGGAGAGTTTTCATCAGGAAGTCTTTTTGTCAGTAAATGGTGTACGTCAggt                        <  2:24433/249‑2 (MQ=255)
                                                                                                                                                                                                                                           ggcggGCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgcccgaccgac             >  2:252844/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      tttctttttctGCCATGAGAAATGAGCGCAGCGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCATGCGACCGACCGACTGAACAATccg  >  2:54436/1‑249 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
TGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGC‑TGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAG  >  NC_000913/4604095‑4604591

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