Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,056,115	0	A	C	30.9%	58.7 / 12.3	44	G75G (GGT→GGG)	glnL	protein histidine kinase NtrB
Reads supporting (aligned to +/- strand):  ref base A (17/12);  new base C (1/12);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 2.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCTTTGCTGAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTCCTGAACAGGTTGC  >  NC_000913/4055865‑4056363                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            aCGCTTTGCTGAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTaaaa                                                                                                                                                                                                                                                             <  2:125340/251‑1 (MQ=255)      ttGCTGAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCTTCACTTCGTT‑ATCGGTAAAA‑CCt                                                                                                                                                                                                                                                         <  1:347212/249‑1 (MQ=255)           gAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTGG‑CC‑CCg                                                                                                                                                                                                                                                <  1:378454/251‑1 (MQ=255)               agcTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTc                                                                                                                                                                                                                                            >  1:119372/1‑251 (MQ=255)                       cgcgccACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACttt                                                                                                                                                                                                                                    >  1:164865/1‑251 (MQ=255)                         cgccACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTACATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCt                                                                                                                                                                                                                                  >  2:374171/1‑251 (MQ=255)                             acgaaacccGCCAGGCGAATTTTTAATCACATGTGCCAGCCCGCGCACTAAATCCCGGGCAGCACCCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGTTGGTTATCCATCGGAGCCATCTCCAGCAGGACCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCTTCACTTCTTT‑ATCGGTAAAC‑CCTTC‑CC‑CCGCCTCCAGACTTACTTGCa                                                                                                                                                                                                                              <  1:124294/244‑1 (MQ=255)                                  aaCCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATTAGc                                                                                                                                                                                                                         >  2:311018/1‑251 (MQ=255)                                      cccAAGCGGTTTTTTACTCTTAGGCGCCAGGCCGCGCACTAAAACTCGGGCAGCATCCGT‑C‑GGGGCGTGCTGTGGCTGTTCCTGACT‑TAGGCGGCGGTGGTTATCCATGGGGGCCACCTCCGGCAGGATCAGGCCGTCCGGCATTTGCTGGGCCGTCACATAATGGATATGCGAGCGCCC‑TTCGATGACC‑ACCTTCTCTT‑GTTGACCGGGAGAG‑CTTCG‑CT‑CCGCCTCCAGACTTTCTTGCATCAGCTCGa                                                                                                                                                                                                                     <  1:195768/250‑1 (MQ=255)                                       ccAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCACCCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑AAAGCGGCGCTGTTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGTATATGCGGGCGCCC‑GTCGATGCCC‑AGCGTCACTTCTTT‑AGCGGTAAAC‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGat                                                                                                                                                                                                                    <  1:211852/251‑1 (MQ=255)                                         aaGCGGATATTTAATCTCATGAGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGATGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGGTATCCATCGGAGCCATCTCCAGCAGGATCAAGCCGTCCGGCAGACGCTGGACCGTCACAGTAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACGTCGTT‑ATCGGGAAAC‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGatat                                                                                                                                                                                                                  >  2:17317/1‑251 (MQ=255)                                          aggggATTTTTAATCTCAGGTGCCAGGCCGCGCACTAAACCACGGGCAGCAACCTG‑C‑GGGGGGGGCTGTAGCTGTTCCTGACT‑TAGGCGGCGCTGGTTACCCACCGGAGCCATCACCAGCAGGATCATGCCGTCCGGCATCCGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCTCTTCGTT‑ATCGGTAAAC‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATAtt                                                                                                                                                                                                                 <  2:164865/248‑1 (MQ=255)                                            ggaatttttactccctgttcccgggCCGCGCACTAAACCACGGGCAGCA‑CCGG‑CGTGGGCGTGCTGGAGCTTTTCCTGTCT‑TAAGCGGCGGTGGTTTCCCATCGGCGCCATCCCCAGCGGGTTCTGGCCGTCCGGCATCCGCTGGGCCGTCACAGAAAGTATATGCGACCCCCC‑GCCGATGACC‑AGCGTCACTTCGTT‑ATCGGAAAAC‑CCTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTa                                                                                                                                                                                                               <  2:390974/228‑1 (MQ=255)                                                         ccATGTGCCAGGCCGCCCAATAAACCACGGGCAGCACCCTG‑C‑GGGGCGTGCTGGAGCTGTTCCTGATT‑TAAGCGGCGCTGGTTATCCACGGGAGCCATCCCCAGCAGGATCATGCCGTCGGGCATCCGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCCCTTCTTT‑ATCGGTAAAC‑CTTTC‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCt                                                                                                                                                                                                  <  1:57708/250‑1 (MQ=255)                                                          cATGGGCCAGGCCCCCCCCTAAATCACGGGCATCAACCGG‑C‑TGGGCGGGCTGTAGCTGTTCCTGTCT‑TAAGCGGCGTTGGTTATCCATCGGAGCCATCTCCAGCAGGATCAGCCCGCCCGGCATCCGCTGGGCCTTCACAGAAAGGATATGCGCGCGCCC‑GTCGAGGACC‑AGCGTCCTTTCTTT‑ATCGGTAAAC‑CCTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                                                                                                                                                                                                    <  1:311018/248‑1 (MQ=255)                                                                  aGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGAAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACATCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTAAACAGTTc                                                                                                                                                                                         >  2:248241/1‑251 (MQ=255)                                                                     ccGCGCACTAAATCCCGGGCAGCACCCTG‑C‑TGGGGGTGCTGTAGCTGTTCCTGACT‑TAGGCGGCGCTGTTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACTCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CTTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCgg                                                                                                                                                                                      <  1:299233/251‑1 (MQ=255)                                                                     ccGCGCACTAAATCACGGGCAGCAAC‑TGTC‑GGGGCGTGCTGTCGCGGTTCCTG‑CTATACGCGGCGCTGGTTATCCGGCGGGGCCATCCCAAGCAGGATCATGCCGCCCGGCATCCGCGGGCCCGTCACAAAAAGGTGATGCGCGCCCCC‑GTCGAGGCCC‑AGCTTCACCTCGTT‑ATCTTTAAAC‑CCTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCgg                                                                                                                                                                                      <  2:419686/251‑1 (MQ=255)                                                                           aCTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGATCCTGACT‑TAAGCGGCGCTGGTTATCAATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGCAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑CTCGGTAACA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACgg                                                                                                                                                                                >  2:25002/1‑251 (MQ=255)                                                                                 ccACGGGCAGCACCCGG‑C‑TGGGCGTGCTTTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCTCCGGGGCCATCCCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAC‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTAcc                                                                                                                                                                          <  1:110327/250‑1 (MQ=255)                                                                                      ggCAGCACCCGG‑C‑GGGGCGTGCTGTAGCGGTTCCTGACT‑TAAGCGGCGCTGTTTATCCACCGGAGCCACCCCCAGCAGGATCAGGCCTCCCGGCATACCCGGGGCCGTCACAGAAAGGAAATGCGAGCGCCC‑GTCGATGACC‑AGCTTCACTTCTTT‑ATCGGAAAAC‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACa                                                                                                                                                                     <  1:263666/251‑1 (MQ=255)                                                                                                     ggggCTTGCTGTAGTTGTTCCTAACT‑TAAGCGGCGCTGGTTTTCCACCGGGGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCTTCATTTGTTT‑ATCGGTAAAA‑CCTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCt                                                                                                                                                        <  1:289593/250‑1 (MQ=255)                                                                                                         cggtctttAGCTGTTCCTCATT‑GAAGCCGCGCTGGTTATCCTTCGGGCCCATGCCCGGCATGACCACCCCCCCCGGCATCCGCGGGGCGGTCACAGAAATGATAGGCGCGCTCCC‑GGCGGTACCC‑AGGGTCCTTTCGTT‑ATCGGGAAAA‑CTTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTTCGGGAGCTTTgg                                                                                                                                                    <  2:184743/244‑1 (MQ=255)                                                                                                             cTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGACTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCg                                                                                                                                                 >  1:432565/1‑250 (MQ=255)                                                                                                                      ttCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgtt                                                                                                                                       >  2:182858/1‑251 (MQ=255)                                                                                                                                 aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAgg                                                                                                                             >  1:20020/1‑251 (MQ=255)                                                                                                                                    cggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGtt                                                                                                                          >  1:160189/1‑251 (MQ=255)                                                                                                                                                            atccccAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTc                                                                                                  >  1:303408/5‑251 (MQ=255)                                                                                                                                                              ctcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGCAAGGATATGCGCGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAACCAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCAt                                                                                                >  2:48349/1‑251 (MQ=255)                                                                                                                                                                      ggatctgcccgcccGGCATACGCGGGGCGGTCACAAAAAGGATA‑GCGGGCCCCCCGGCGA‑GACCCGGCGTCCCTTGGTT‑ATTGG‑AAAA‑CCTTG‑CCGCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAAGTTGCGGGAGCTGTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAAc                                                                                        <  2:195072/239‑1 (MQ=255)                                                                                                                                                                      ggaccATGCCGCCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGGGCGCCC‑GTCGATGCCC‑AGCGTCACTTCGTT‑ATCGGTAAAC‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAAc                                                                                        <  1:213248/247‑1 (MQ=255)                                                                                                                                                                        aTCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACa                                                                                       >  2:256880/1‑250 (MQ=255)                                                                                                                                                                                     gtatcctgtGGGCCGTCGCAGAAAGGATAGGCGACCCCCC‑GTCGATCGCC‑AGCTTCACTTCGTT‑ATCGGTAAAC‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATc                                                                         <  2:77969/243‑1 (MQ=255)                                                                                                                                                                                            cTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCTTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTGGCTCaa                                                                                                                                                                                               <  1:134488/126‑1 (MQ=255)                                                                                                                                                                                            cTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCaa                                                                                                                                                                                               >  2:134488/1‑126 (MQ=255)                                                                                                                                                                                             tGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGCCC‑AGCGACACTTCCTT‑ATCGGTAAAA‑CATTG‑CC‑CCGCCTCCCGACTTACTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGATGTCATCGATTAACAAACTACTGTTAATCAGCGAGtt                                                                 >  1:198106/1‑251 (MQ=255)                                                                                                                                                                                                          aaaGGATAGGCGAGCGCCG‑GTAGATGACC‑AGCGTCCTTTCTTT‑TTCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCa                                                    <  1:439386/251‑1 (MQ=255)                                                                                                                                                                                                                 atGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGGGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCgg                                             >  1:358251/1‑251 (MQ=255)                                                                                                                                                                                                                       gagcgccctgcaatgcccatcttCACTTCTTT‑A‑CGGTAAAAACTTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCTATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGTCCAGCATCGGGCTg                                         <  1:188786/229‑1 (MQ=255)                                                                                                                                                                                                                        gcccccgtgcgagaaccggCGTCCCTTTGTT‑CTCGGAAAAT‑CCTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGc                                        <  2:303408/233‑1 (MQ=255)                                                                                                                                                                                                                         gcccccgtggataacaatcGTCACTTCTTT‑ATCG‑TAAAA‑CCCTGTCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGGAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTGGTCATCGATTAACAGAATACTGTTAATCAGCGAGTTGAGGATCGGCCCAGCATTGGGCTGCg                                       <  2:289833/233‑1 (MQ=255)                                                                                                                                                                                                                                      gcccggcgcctttttgttattgttAAAC‑CCTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGAGCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCa                           <  2:397085/228‑1 (MQ=255)                                                                                                                                                                                                                                       caggCGTCACTTCGTA‑GACGGTAAAA‑CTTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATa                         <  2:20020/248‑1 (MQ=255)                                                                                                                                                                                                                                               aCTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGtc                 >  2:101897/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cGTT‑AGCGGTAAAA‑CCTGG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTCCt             <  1:244988/251‑1 (MQ=255)                                                                                                                                                                                                                                                               aaa‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTACCTCAACAGCTCCTGTAACGGTGTACCAAACAATTTGCGGAAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCCTCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTCCTGAACAGGTTGc  >  1:273313/1‑251 (MQ=255)                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            ACGCTTTGCTGAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTG‑C‑TGGGCGTGCTGTAGCTGTTCCTGACT‑TAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCC‑GTCGATGACC‑AGCGTCACTTCGTT‑ATCGGTAAAA‑CCTTG‑CC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTCCTGAACAGGTTGC  >  NC_000913/4055865‑4056363

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