Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT31.2% 90.0 / 13.7 48S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (4/29);  new base T (15/0);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 3.55e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCA‑ACGCTGGAGATTTTTCATCAGGAACTCTTTTTG  >  NC_000913/4604107‑4604546
                                                                                                                                                                                                                                             |                                                                                                                                                                                                             
gcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttc                                                                                                                                                                                                               <  1:283050/238‑1 (MQ=255)
       gTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        <  2:790056/238‑1 (MQ=255)
                                          gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCCCGTTTGCGGCGGGCGATTTTTCTTTTTCTTCCTTTAGAAAGGAACGCCGAGAAGGAGTAt                                                                                                                                                                         >  2:726802/1‑234 (MQ=255)
                                                              aGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATTAGAAGGTGGCTGCGAGAGGGGTTTTCTTTCCTTGTGGAACT‑GGCgggtta                                                                                                                                               >  2:859297/1‑234 (MQ=255)
                                                                          aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                                                                                      <  1:305208/173‑1 (MQ=255)
                                                                                aaTGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTAccc                                                                                                                                 <  1:814124/235‑1 (MQ=255)
                                                                                aaTGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                             <  1:138327/239‑1 (MQ=255)
                                                                                    tCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                                                                                                          <  1:432969/159‑1 (MQ=255)
                                                                                    tCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGGCTTTTTTTTTTTTTTTCCCTGAGaaat                                                                                                                                                                                            >  2:305208/1‑171 (MQ=255)
                                                                                          tCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                                                                                   <  1:254499/239‑1 (MQ=255)
                                                                                           cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                                                                                  <  1:378259/239‑1 (MQ=255)
                                                                                              aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCATtag                                                                                                                                                                                                >  2:432969/1‑159 (MQ=255)
                                                                                                 aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                                                                                                <  1:622779/235‑1 (MQ=255)
                                                                                                  aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                                                                                           <  2:302684/239‑1 (MQ=255)
                                                                                                                   taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGGAGAGAGGGGTTATCTTGTCTGGCTGAACT‑TGGGGTTAAGTTTTCTCGTACCCGGCGGTTAGTTGTTTTGGATGATGATAATTAAAGTCGGc                                                                                          >  2:837834/1‑239 (MQ=255)
                                                                                                                       ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAt                                                                                      <  1:859297/239‑1 (MQ=255)
                                                                                                                                 gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        >  1:595117/1‑116 (MQ=255)
                                                                                                                                 gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        <  2:595117/116‑1 (MQ=255)
                                                                                                                                           aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                                                                                      <  1:859639/108‑1 (MQ=255)
                                                                                                                                               tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                                                                                              <  1:528192/128‑1 (MQ=255)
                                                                                                                                                     ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCACGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAgg                                                                                                                                                                                            >  2:859639/1‑107 (MQ=255)
                                                                                                                                                         aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGAAGGGGGGGCCGGGAGGGGTTAtttttg                                                                                                                                                                    >  2:528192/1‑126 (MQ=255)
                                                                                                                                                           ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                                                                                   <  1:265274/111‑1 (MQ=255)
                                                                                                                                                               ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCAt                                                                    <  1:302618/216‑1 (MQ=255)
                                                                                                                                                                taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                                                                                      <  2:823188/87‑1 (MQ=255)
                                                                                                                                                                taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                                                                                      >  1:823188/1‑87 (MQ=255)
                                                                                                                                                                    cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                                                                                              <  1:293454/107‑1 (MQ=255)
                                                                                                                                                                    cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGt                                                                                                                                             <  1:572455/139‑1 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttttttttct                                                                                                                                                                                                        >  2:491781/1‑80 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        <  1:491781/80‑1 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        >  1:781460/1‑80 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        <  2:781460/80‑1 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                                                                              <  1:211502/122‑1 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                             <  2:311556/154‑1 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                             >  1:311556/1‑154 (MQ=255)
                                                                                                                                                                     tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCGGAGAGGGGtttt                                                                                                                                                                         >  2:265274/1‑109 (MQ=255)
                                                                                                                                                                      ccATACCTATGCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTAAAGCGt                                      <  1:837834/239‑1 (MQ=255)
                                                                                                                                                                           ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGgt                                                      <  1:434896/218‑1 (MQ=255)
                                                                                                                                                                            cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                                                                                 <  1:324527/175‑1 (MQ=255)
                                                                                                                                                                              ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCCGGGAGGGGGTATTTTGCCTGACTGAAATGGGGGATTAAGTTTcccggaa                                                                                                                                    >  2:572455/1‑134 (MQ=255)
                                                                                                                                                                               taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTaa                                           <  1:48006/225‑1 (MQ=255)
                                                                                                                                                                                    gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGAc                                                         <  1:441005/206‑1 (MQ=255)
                                                                                                                                                                                    gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGt                                             <  1:116924/218‑1 (MQ=255)
                                                                                                                                                                                      cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGAGGGGGTATCTTGCCTGACTGAAAT‑GGCGGTTAAGTTTTCCCGGGACCGGCGGTTTTTTTTTTCGGAGGAGGATATTTAAAATCGGcggg                                                                                       >  2:324527/1‑172 (MQ=255)
                                                                                                                                                                                         gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTTAACT‑GGCGATTAAGTTTTCGCGTACCCCGCGTTTAGTTTTTTCGGGTGTTGATATTTAAGACCGGCTGGTTGGGTTATTTTGT‑CCCGATCCCCTGGGCGGGGTGTTAAGTGAAGCGTGA‑ACg                                >  2:48006/1‑225 (MQ=255)
                                                                                                                                                                                             cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAGGGGTTTTCTTGCCTGTCTGGAAT‑GGCGATTAAGTTTTCGCGTACCCGGCGGTTTGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGGTCCTTGTGG‑CCGTCACCACCGGACGGCGGATTTAGTTAAACGGCACACGCGGGGGAtatttatcca               >  2:225830/1‑234 (MQ=255)
                                                                                                                                                                                              aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTGGAAGTGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACT‑GGCGGTTAAGTTTCCGCGTACCCGGCGGTTTGTTTTTGCGGATGGTGATATTTAAGCTCGGCTGGTTGTGTCATTGTTG‑CCATCACCGCTGGACGGTGGTTTaa                                               >  2:441005/1‑206 (MQ=255)
                                                                                                                                                                                              aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACT‑GGCGGTTTAGGTTTCGCGTACCCCGGGGTTTGTTATTGCGGGGGGTGATATTTAAACTCGGCTGGTTGGTTTATTTTTGCCCTTCACC‑CCGGGCGGGGTATTTAGTTAAACGGCa‑a                                  >  2:116924/1‑218 (MQ=255)
                                                                                                                                                                                                     gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                                                                             <  1:893182/59‑1 (MQ=255)
                                                                                                                                                                                                       tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGAc                                                         <  1:148659/187‑1 (MQ=255)
                                                                                                                                                                                                       tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATGTTTTCTGTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCA‑ACGCTGGAGATTTTTCATCAGGAACTCttt     <  1:872977/238‑1 (MQ=255)
                                                                                                                                                                                                          cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCA‑ACGCTGGAGATTTTTCATCAGGAACTCTTTTTg  <  1:520373/238‑1 (MQ=255)
                                                                                                                                                                                                               cgcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCaga                                                                                                                                                                                   >  2:893182/1‑59 (MQ=255)
                                                                                                                                                                                                                 cgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCCGAGGGGGGTTTTCTTGTCTGACTGAACT‑GGCGGTTAAGGTTTCGCGTACCCGGCGGTTTGTTTTTTCGGGTGGTGATATTGAAGCTCGGCTGGTTTGTTCATTTTGT‑CCATCACCCCTGGGCGgtgggttta                                               >  2:148659/1‑181 (MQ=255)
                                                                                                                                                                                                                  ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTAcc                                                                                                                                  <  1:664198/104‑1 (MQ=255)
                                                                                                                                                                                                                            tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCcg                                                                                                                                        >  1:164563/1‑88 (MQ=255)
                                                                                                                                                                                                                            tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCcg                                                                                                                                        <  2:164563/88‑1 (MQ=255)
                                                                                                                                                                                                                                             ctttttctGCCATGAGAAGTGAGAGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACg                                                        >  2:713373/1‑150 (MQ=255)
                                                                                                                                                                                                                                             ctttttctGCCATGAGAAGTGAGAGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACg                                                        <  1:713373/150‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                             
GCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCA‑ACGCTGGAGATTTTTCATCAGGAACTCTTTTTG  >  NC_000913/4604107‑4604546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑