Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,642,761	C→T	A124A (GCG→GCA)	rsmJ ←	16S rRNA m(2)G1516 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,642,761	0	C	T	91.7%	121.9 / NA	49	A124A (GCG→GCA)	rsmJ	16S rRNA m(2)G1516 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (23/21);  minor base G (0/4);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATT  >  NC_000913/3642522‑3643005                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      aGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCAc                                                                                                                                                                                                                                            <  1:178960/251‑1 (MQ=255) aGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGGGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCACCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGGCGTCGCGTAGTGCGGCAACCAc                                                                                                                                                                                                                                            >  2:234801/1‑251 (MQ=255)       ccATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCGTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGAtt                                                                                                                                                                                                                                      <  1:372282/251‑1 (MQ=255)              cATCAAGATCCGGTCCCACCAGCGACTGAAAGACTCGCTTCTCTTTTTTCACTAGCTCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCGGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTcc                                                                                                                                                                                                                               <  1:251797/251‑1 (MQ=255)                   aGATCCGGTCCCACCAGCGACTGAAAGACCCGCATCTCTTTTTTCACTAGCCCGCTTTTCTGCTTATGCGGGACCATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGAGGCAACCACTGCATTACGCTCCAGCac                                                                                                                                                                                                                          >  2:212962/1‑250 (MQ=255)                     aTCCGGTCCCACCAGCGACTTAAAGCCACTCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCCAAATAAACCCCCTGCGGGCGCGGGGTAATACCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATGTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATc                                                                                                                                                                                                                         <  2:146303/250‑1 (MQ=255)                         ggTCCCCCCAGCGCCGGAAAGACACGCAGCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACAGCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGGCTAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCAc                                                                                                                                                                                                                    <  2:175536/251‑1 (MQ=255)                            cccACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCGCCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCAcgcg                                                                                                                                                                                                                 <  2:102564/251‑1 (MQ=255)                             ccaccaGCTACCAAAAGACACGCTTCTCTTTTTGCACTAGCGCTCTTTTCGGCTTATGCGGGAACATCTGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGCCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTTCACCCAGCCGCCGATTTCCGCACCCGCATAACCACGCGCCAGGCCGTCGTCGAGGATTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCgg                                                                                                                                                                                                                <  2:200158/251‑1 (MQ=255)                                ccaGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGc                                                                                                                                                                                                               >  1:296248/1‑249 (MQ=255)                                    cGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGTGTCGAGATAAACCACCTGCGGGCGCGGGGTAATAACAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATCACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCCCTGGATTACGCTCCAGCATCCGCACGCGACAGCCgc                                                                                                                                                                                                         >  2:395865/1‑250 (MQ=255)                                          aaaGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCCCCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAg                                                                                                                                                                                                   <  1:233446/251‑1 (MQ=255)                                                 cGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCCCGCGGCAGCCGACTGAAGCCAGTAc                                                                                                                                                                                            >  2:124720/1‑251 (MQ=255)                                                          tttttCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCAt                                                                                                                                                                                    <  2:213117/250‑1 (MQ=255)                                                                aCTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCACCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCCACCACTGGATTACGCACCAGCATCAGCACGCGGCAGCCGACTGAAGCCAGTCCAAAGGCATCGCGccc                                                                                                                                                                             >  2:18123/1‑251 (MQ=255)                                                                 cTAGCGAGCTTTTATGCTAATGCGGGAACATAGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGACGTCAGGCTGGAGGCGTGACTTAACGGCAA‑CCGCTCCAGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCAAGGCCGTCGTCGAGTCGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGCCTGAAGCCAGTACAAAGGCATCGCGccc                                                                                                                                                                             >  2:278825/1‑250 (MQ=255)                                                                       cgcTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCGGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACTAAGGCATCGCGCCCAAGT‑CCt                                                                                                                                                                     <  1:83016/251‑1 (MQ=255)                                                                           tttCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAACCAGTACAAAGGCATCGCGCCCAAGT‑CCTGc                                                                                                                                                                   >  2:119772/1‑249 (MQ=255)                                                                             tCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGt                                                                                                                                                                <  1:85178/250‑1 (MQ=255)                                                                             tCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTg                                                                                                                                                               >  2:128894/1‑251 (MQ=255)                                                                                gCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGc                                                                                                                                                             >  1:132688/1‑250 (MQ=255)                                                                                  ttATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATc                                                                                                                                                          >  2:100696/1‑251 (MQ=255)                                                                                              aTCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTACCTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCAAGGACGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGA‑CCAGCAGTGGCATCCACCACATcagc                                                                                                                                              >  1:165547/1‑248 (MQ=255)                                                                                                cGGGACGAAATAAACCCCCTGGGTGCGTGGGTTAATATCCGCCAGCGCGTTAAGGCTGGAGGCGTGAATTAAC‑GCAAACCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATACCCCCGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGACGCCAGTACAAAGGCATCGCGCCCAAGT‑TCTGAAGTGGCATCCACCTCATCCGGCa                                                                                                                                            <  2:239284/251‑1 (MQ=255)                                                                                                   gccggtataatcccccTGGGGGGGCGGGGTAATATCAGTCAGGGCCGCCAGGCTGGAGGCTTGAATACCCTGCAA‑CCGCCCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCGCTGGATTACGCTCCAGCCTCCGCACGCGGCAGCCGGCGGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTTCAGTGGCATCCACCACATCCGGCAaat                                                                                                                                         <  2:246988/237‑1 (MQ=255)                                                                                                   gTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAaat                                                                                                                                         >  1:139373/1‑251 (MQ=255)                                                                                                          tAAACCACCTGCGGGCGCGGGGTAATATCAGTCAACGCCGTCAAGCTGGAGGCGTGAATTAACTGAAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCACCCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATCATCGcc                                                                                                                                  >  2:383263/1‑251 (MQ=255)                                                                                                            aaCCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGA‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCtt                                                                                                                                >  1:309271/1‑251 (MQ=255)                                                                                                                accTGCGGGCGCGGGGGAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAAt                                                                                                                            >  1:84429/1‑251 (MQ=255)                                                                                                                 ccTGCGGGGGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGGGAATTAACTGCAA‑CCGCACCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAATCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATg                                                                                                                           <  1:278825/251‑1 (MQ=255)                                                                                                                       ggCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGGCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCAcc                                                                                                                     <  1:284166/251‑1 (MQ=255)                                                                                                                           cgGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCtt                                                                                                                 >  1:281378/1‑251 (MQ=255)                                                                                                                              ggTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCGGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCg                                                                                                              <  1:378066/251‑1 (MQ=255)                                                                                                                                                                 gAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCAGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCAGCAAATAATCGCCTTTAATGACCACCGCTTTCGCCACCGCCTCACCACGACCCCCGCCGACTTTGCgc                                                                           >  2:328538/1‑250 (MQ=255)                                                                                                                                                                             a‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATc                                                               >  1:386908/1‑251 (MQ=255)                                                                                                                                                                              cccGCTCCTGCAACCAGCCGCCGATTTCCGTATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCTCCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATc                                                               <  1:18123/250‑1 (MQ=255)                                                                                                                                                                                 gCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCCGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGCATTTCCGTCGGTGCGCCCTCGCt                                                            >  2:15346/1‑251 (MQ=255)                                                                                                                                                                                    ccTGCAACCAGCCCCCGTTTTCCGCATCCGCATACCCGCGCGCCAGGCCTTAGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCCCGCGGCAGCCGACTGAAGCCAGTACAAAGGCATGGCGCCCAAGT‑CCTGCAGTGGCTTCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCTCCATCGCTc                                                           <  1:347328/249‑1 (MQ=255)                                                                                                                                                                                      ggcaccaaGCCGCCGATTTCCGCACCCGCCTAACCACGCGCCAGGCCGGCGTCGAGTAGGGCGGAAACCACTGGATTCCGCTCCAGCATCCGCACGGGGCAGCCGACTGAAGCCAGTACAAAGGCAGCGCGCCAAAGT‑CCTGCAGTGGCATCCACCACATCCGGAAAATAATCGCCTTTAATCCCCACCGCTTTGGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGcc                                                       <  2:106949/244‑1 (MQ=255)                                                                                                                                                                                              gccgccGATTTCCGCATCCGCATAACCACGCGCCAGTCCGTCGTCGAGTAGGCCGGCAACCACGGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACGTAAGCCAGTACAAAGGCATCGCGCCCAAGC‑CCTGCATTGGTATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTGTCGCCACCGCCTCACCCCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAAt                                               <  1:171063/251‑1 (MQ=255)                                                                                                                                                                                                         ccccacccgcatAACCACGCGCCAGGCCGGCGTCGAGTAGGGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCCCCACATCCGGCAAATAATCGCCTTTAAGGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATg                                    <  2:287518/245‑1 (MQ=255)                                                                                                                                                                                                          cgcatccgcatAACCACGCGCCAGGCCGTCGTCGATTAGTGCGGCAACCGCTGGATTACGCTCCAGCATCCGCACGCGGCGGCCGACTGAAGCCAGTACAAAGGCTTCGCGCCCAAGT‑CCGGCAGTGGCACCCCCCACATCCGGCAAATATTCGCCTTTAATGCCCCCCGTTTTCGCCCCCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATgc                                   <  1:394084/251‑1 (MQ=255)                                                                                                                                                                                                                cgcatACCCACGCGCCGGGCCGTCGTCGAGTAGTGCGGCACCCCCTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCaa                             <  1:234801/251‑1 (MQ=255)                                                                                                                                                                                                                     aaCCACGCGCCAGGCCGTCGTCGAGTACGGCGGCAACCACAGTATGACGCCCCAGCTTCCGCGCGCGGCAGCCGCCTGAAGCCAGTCCAAAGGCACCGCGCC‑TAGTGCGTGCAGTGGTATCCACCACAGCCGTCAAAAAATCGCCGTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATGTGCGTCGGTGCGTCATCGCTCCGTCAACAAAATCAACAAAGATGCCGCCAAGtttt                        <  2:74517/251‑1 (MQ=255)                                                                                                                                                                                                                             gcCAGGCCGTCGTCGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCtcatc                >  1:278487/1‑251 (MQ=255)                                                                                                                                                                                                                               cAGGCCGCCGCCGAGTTGCGCGACACCCACTGGATTCCGCCCCCGCCCCCGCAGGGGGCGGCCGCCTGAAGCCAGGACAAAGGCACCCCGCCCAAGT‑CCTGCAGTGGCACCCACCCCATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACTAAATCAACAAAGATGCCGCCAAGGGTTGGCTCATCAc              <  2:146715/251‑1 (MQ=255)                                                                                                                                                                                                                                     gtcgccGAGTAGTGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCACACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATACCGCCAAGATTTGGCtcatcaagcctgc        >  1:339566/1‑244 (MQ=255)                                                                                                                                                                                                                                       cgccgAGTAGTGCGGCAACCACTGGATTCCGCTCCAGCATCCGCCCGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTgcgc      <  1:119772/248‑1 (MQ=255)                                                                                                                                                                                                                                           gAGTAGTGCGCCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGCATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGCCTCATCACGCTTGCGCAAtt  >  2:189160/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      AGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACCACCTGCGGGCGCGGGGTAATATCAGTCAGCGCCGTCAGGCTGGAGGCGTGAATTAACTGCAA‑CCGCTCCTGCAACCAGCCGCCGATTTCCGCATCCGCATAACCACGCGCCAGGCCGTCGTCGAGTAGCGCGGCAACCACTGGATTACGCTCCAGCATCCGCACGCGGCAGCCGACTGAAGCCAGTACAAAGGCATCGCGCCCAAGT‑CCTGCAGTGGCATCCACCACATCCGGCAAATAATCGCCTTTAATGCCCACCGCTTTCGCCACCGCCTCACCGCGACCACCGCCGAATTTGCGTCGGTGCGCCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATT  >  NC_000913/3642522‑3643005

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