Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,515,879	C→T	G12E (GGG→GAG)	dam ←	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,515,879	0	C	T	98.1%	162.2 / NA	53	G12E (GGG→GAG)	dam	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (26/25);  minor base G (0/1);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTCTTC  >  NC_000913/3515634‑3516118                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                aCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTCCTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCTCCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATTCTTGCCCTCTGcc                                                                                                                                                                                                                                            <  1:279063/251‑1 (MQ=255)           aatcTCGTCGCAATGTGTTTCGGCAACAAACAGATCGCGTGCGGCCGGTACGTACTCAGCAGTACGCATCTTCACAATGTTATAGAGACTGATCATGTCGCGATTGATAACGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCTACTCGGCACCTACAAAAGGATCAACCAGACAGTCTCCCGTGGGCAAATGCCGTTTAATATCATCAATCAGGGGATACGTGCCCTCTGTCCACtgctaaa                                                                                                                                                                                                                                  <  2:413769/248‑7 (MQ=255)                    gcAAGTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGGACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATAGAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCGGCCCACTTCAAAAAATCGCGAt                                                                                                                                                                                                                         <  2:392480/251‑1 (MQ=255)                     cAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGAtt                                                                                                                                                                                                                        <  1:103789/251‑1 (MQ=255)                            tttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTttcttc                                                                                                                                                                                                                 >  1:339274/1‑251 (MQ=255)                                      aaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGact                                                                                                                                                                                                        >  2:311294/1‑250 (MQ=255)                                              cgcgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAAtt                                                                                                                                                                                                >  1:100747/1‑250 (MQ=255)                                                cgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTac                                                                                                                                                                                              >  1:11548/1‑250 (MQ=255)                                                 gTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTTAACCAGACAGTCGCCCTTGTGCAAATGCCGTTTAATATCATCAAGCAGTTGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTacac                                                                                                                                                                                            <  1:379749/251‑1 (MQ=255)                                                 gTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTacac                                                                                                                                                                                            <  2:232279/251‑1 (MQ=255)                                                    cGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCtt                                                                                                                                                                                         >  1:177554/1‑251 (MQ=255)                                                                acTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGCTACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTgg                                                                                                                                                                             >  1:386654/1‑251 (MQ=255)                                                                              aTCTTCACAATGTTATAGAGACTGATCAGGTCGCTTTTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACACGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATg                                                                                                                                                                                                              >  1:243241/1‑204 (MQ=255)                                                                              aTCTTCACAATGTTATAGAGACTGATCAGGTCGCTTTTGATAGCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAGTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGTAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCGTCATg                                                                                                                                                                                                              <  2:243241/204‑1 (MQ=255)                                                                              aTCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAATCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGCCTAACTAATTAAACCTTCTCCGGCTGTGGAGAAAGCTCCGAcc                                                                                                                                                               >  2:373599/1‑250 (MQ=255)                                                                                     caATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATcc                                                                                                                                                         >  1:81420/1‑250 (MQ=255)                                                                                        tGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACtt                                                                                                                                                                                                                                       <  2:396229/169‑1 (MQ=255)                                                                                        tGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACtt                                                                                                                                                                                                                                       >  1:396229/1‑169 (MQ=255)                                                                                              agagACTGATCAGGGCTCTATTGAGATCGGCAAGGATAGAACGAGAAAAGGCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCGTGGGCAAATGCCGTTTAAGATCATCAAGCAGGGGATACGTGCCCTCTGCCCTCTTCAAAAAAGCGCGATTTTTCGTCATGCTTACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACATGATCCATCGCTTCAg                                                                                                                                               <  2:362275/251‑1 (MQ=255)                                                                                                             cGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACg                                                                                                                                >  1:299828/1‑251 (MQ=255)                                                                                                             cGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCAGCCAGCGCTTCATGCAACTCCCGCAACg                                                                                                                                <  2:272903/251‑1 (MQ=255)                                                                                                                 atttatatCGGCAAGGATATAACGAGAAAAGTCGTTGTTGAGAAACACCGACCCGGCCCCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAATCAGGGGATACTTGCCCTCTGCCCTCTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACTCCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGAt                                                                                                                              <  2:267809/245‑1 (MQ=255)                                                                                                                   tGATATCGGCAAGGATATAACGAGAAAAGTAGGTGTTGAGAAACACCGACCCGGCCCCTACACAAGGCCCAACCAGACAGTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCAGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTCCACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTCCAGGCAACTCCCGTAACGATaaa                                                                                                                           <  2:340069/250‑1 (MQ=255)                                                                                                                    gATATCGGCAAGGATATAACGCGAAAAGTCGGTGTTGAGAAACACCGACCCGGCCCCTACAAAAGGCTCAATCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATCTTTGCCCGCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAAcc                                                                                                                         <  1:4212/251‑1 (MQ=255)                                                                                                                       atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCAtt                                                                                                                      >  1:363360/1‑251 (MQ=255)                                                                                                                       atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCAtt                                                                                                                      <  1:289318/251‑1 (MQ=255)                                                                                                                              aGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCggg                                                                                                                 >  2:50064/1‑249 (MQ=255)                                                                                                                               ggATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGgaga                                                                                                              <  1:311294/251‑1 (MQ=255)                                                                                                                                      aCGAGAAAAGTCGGTGTTGAGAAACACCTACCCGGCACCTACAAAAGGCTCAACCAGTCATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAATCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCATCTCCCGTAACGATAAACCATTGCGGGAGATATCCtt                                                                                                       <  2:143461/251‑1 (MQ=255)                                                                                                                                          gaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCGACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCAGACAGCATCCCGAGCTTCAGGCACCTCCCGCCACGATAACACATTGCGGGAGATAACCTTGAtt                                                                                                   >  2:227289/1‑251 (MQ=255)                                                                                                                                                     gtTGAGAAACACCGACCCGGCACCTACACAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTc                                                                                              >  2:115753/1‑245 (MQ=255)                                                                                                                                                     gtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATc                                                                                         <  2:97439/250‑1 (MQ=255)                                                                                                                                                       tGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGc                                                                                      >  1:206779/1‑251 (MQ=255)                                                                                                                                                            aaCACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTa                                                                                 <  2:81420/251‑1 (MQ=255)                                                                                                                                                             acacCGACCCGGCACCTACAAAGGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGAAACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCTTCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATATCCTTGATTCCTTCAGATCGGCCTGTAc                                                                                <  2:358737/251‑1 (MQ=255)                                                                                                                                                                     ccGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAACAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGc                                                                         >  2:303360/1‑250 (MQ=255)                                                                                                                                                                          cacgaCAAAAGGCTCAACCATACATTCGCCCTTTGGCAAATGCCGTTTAATATCATCAAGTGGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGTCTAACTAATTACACCTACTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCGGTACCTGACGCAGCGGt                                                                   <  2:407357/247‑1 (MQ=255)                                                                                                                                                                              aCAAAAGGCTCAACCAGACATTCGCCCGTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTGCAGGCAACTCCCGCAACGATAAACCAGTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTc                                                                <  2:199139/250‑1 (MQ=255)                                                                                                                                                                              aCAAAAGGCTCAACCAGACATTCGCCCGTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTGGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCTACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATATCCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTc                                                                <  2:368217/250‑1 (MQ=255)                                                                                                                                                                                              gACATTCGCCCTTGGGCAACTGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGTCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGc                                               >  1:229183/1‑251 (MQ=255)                                                                                                                                                                                                      cccTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTTCTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACAt                                       <  2:269298/251‑1 (MQ=255)                                                                                                                                                                                                       ccTTGGGCAAATGCCGTTTAATATCATCAAGCTGGGGATATTGGCCCCCTGCCCACTTCAAAAAAGCGCGATGTTTCTGCATTCGGACTAACTAATTACACCGTCACCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTGCAGGCAACTTCCGCAGCGATAAACCATGGCGGTTGATTACTTTGATTACTTCAGATCGGCTGGGTCCTTACGCGGCGTTTTCGGCCACTTGTTTTTGGCGTGGACATc                                      <  2:263478/251‑1 (MQ=255)                                                                                                                                                                                                           gggCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATctgct                                  >  1:162938/1‑251 (MQ=255)                                                                                                                                                                                                                  tGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTACTGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATg                           >  2:215350/1‑251 (MQ=255)                                                                                                                                                                                                                           atatCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATCCCg                  >  1:10288/1‑251 (MQ=255)                                                                                                                                                                                                                           atatCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCg                  <  2:31598/251‑1 (MQ=255)                                                                                                                                                                                                                           atatCATCAAGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCg                  >  1:361442/1‑251 (MQ=255)                                                                                                                                                                                                                             atcatcAAGCAGGGAATAATTGCCCTCTGCCCACTTCAAAAAAGCGCGCTTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGATCCGACAGCATCCAGCGCTTCAGGCAACTCACGCAACGATAAACCATTGCGGGAGATCACCTTGATTACTGCAGATCGGACTGTACCTGACGCAGCGGTTTTGACCACGGGTTTTTGGCCTGGAGATCTGCTGGCAATGGAGATTCCGCt                >  2:235420/1‑251 (MQ=255)                                                                                                                                                                                                                               catcaAGAAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCttt              <  2:101859/251‑1 (MQ=255)                                                                                                                                                                                                                                    aGCAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATCACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCg         >  1:34379/1‑251 (MQ=255)                                                                                                                                                                                                                                      cAGGGGATACTTGCCCTCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTCCACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCAAGCGCTTCAGGCAACTCCCGCAACGATAAACAATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTTGCAATGTAGATACCGCTTTTTTctcc       >  2:36490/1‑248 (MQ=255)                                                                                                                                                                                                                                        gggtaaaCTTGCCCTCTTCCCCCTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCctc     <  2:324350/245‑1 (MQ=255)                                                                                                                                                                                                                                           gATACTTGCCCTCTGCACACTTCAAAAAAGCGCGCTTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGCCATCTGCTGGCAATGTAGATACCGCTTTTTTTGCAtcttc  >  1:124763/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                ACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTCTCCGGCTGTGGAGAAAGCTCCGACAGCATCCAGCGCTTCAGGCAACTCCCGCAACGATAAACCATTGCGGGAGATAACCTTGATTACTTCAGATCGGCCTGTACCTGACGCAGCGGTTTCGCCCACGGGTTTTTGGCCTGGACATCTGCTGGCAATGTAGATACCGCTTTTTTCGCCTCTTC  >  NC_000913/3515634‑3516118

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