Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT31.2% 89.2 / 13.5 48S64F (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.

GTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCA‑ACGCTGGAGATTTTTCATCAGGAACTCTTTTTG  >  NC_000913/4604114‑4604546
                                                                                                                                                                                                                                            |                                                                                                                                                                                                       
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)
                                                                             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                                                                                                                                                                                                        <  2:595117/116‑1 (MQ=255)
                                                                                                                          gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        >  1:595117/1‑116 (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:781460/1‑80 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        <  2:781460/80‑1 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                        <  1:491781/80‑1 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                                                                              <  1:211502/122‑1 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                             >  1:311556/1‑154 (MQ=255)
                                                                                                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCg                                                                                                                             <  2:311556/154‑1 (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                                                                                                                                        <  2:164563/88‑1 (MQ=255)
                                                                                                                                                                                                                     tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCcg                                                                                                                                        >  1:164563/1‑88 (MQ=255)
                                                                                                                                                                                                                                      ctttttctGCCATGAGAAGTGAGAGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACg                                                        >  2:713373/1‑150 (MQ=255)
                                                                                                                                                                                                                                      ctttttctGCCATGAGAAGTGAGAGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACg                                                        <  1:713373/150‑1 (MQ=255)
                                                                                                                                                                                                                                       tttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACg                                                        >  1:800029/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACg                                                        <  2:800029/149‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                       
GTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCG‑CCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCA‑ACGCTGGAGATTTTTCATCAGGAACTCTTTTTG  >  NC_000913/4604114‑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: ‑