Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,099,793	G→A	G312G (GGC→GGT)	ansB ←	L‑asparaginase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,099,793	0	G	A	97.1%	229.1 / NA	72	G312G (GGC→GGT)	ansB	L‑asparaginase 2
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (40/28);  minor base T (0/2);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAGCCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGGATGCGTTAGCGTAGTTATAAACAATGCCGACTTTCGGCAGTTCATTCA  >  NC_000913/3099549‑3100040                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        ttCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTAt                                                                                                                                                                                                                                                   >  2:187873/2‑251 (MQ=255)              tCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATcc                                                                                                                                                                                                                                       >  1:375424/1‑250 (MQ=255)                 cccGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACtt                                                                                                                                                                                                                                   >  2:549096/1‑251 (MQ=255)                    gAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCgg                                                                                                                                                                                                                                >  1:503474/1‑251 (MQ=255)                         aCCTCACGCGGAATAATATATTCCGGTGATTATCATTAGGAAAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGCATTTCGCATCATCCCCTTCGGCATcc                                                                                                                                                                                                                           >  2:646472/1‑251 (MQ=255)                                       aatTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACAGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGcc                                                                                                                                                                                                             >  1:642203/1‑251 (MQ=255)                                        atTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGccc                                                                                                                                                                                                            >  1:287568/1‑251 (MQ=255)                                             ttCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGGGGTAGCGCCCGTCgg                                                                                                                                                                                                       >  1:474197/1‑251 (MQ=255)                                                     aTTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCggg                                                                                                                                                                                               >  2:230974/1‑251 (MQ=255)                                                        aTCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGga                                                                                                                                                                                              >  1:252245/1‑249 (MQ=255)                                                        aTCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGga                                                                                                                                                                                              >  1:466319/1‑249 (MQ=255)                                                             tAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGACCCTGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGGTCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGACAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCTCCCGTAGGTACGCGGGACGAACGc                                                                                                                                                                                       >  1:522200/1‑251 (MQ=255)                                                                ggTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACtt                                                                                                                                                                                                                                   >  1:435092/1‑204 (MQ=255)                                                                ggTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACtt                                                                                                                                                                                                                                   <  2:435092/204‑1 (MQ=255)                                                                    aTTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGAGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTg                                                                                                                                                                                >  2:322730/1‑251 (MQ=255)                                                                            tGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAATGCACGACTGCAGTACCg                                                                                                                                                                        >  1:273071/1‑251 (MQ=255)                                                                             ggCAATGGACGTGAGTCTTAAGTGAAAA‑GCCCCGGCGCGGTACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGATTCAGCGTGCCAGAGGCGTCGAAGCCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGTAAGAACGCCCTACTGCATTACCGGt                                                                                                                                                                      <  2:55476/251‑1 (MQ=255)                                                                               cAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGtt                                                                                                                                                                     >  2:472457/1‑251 (MQ=255)                                                                                 aTGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATACTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTAACGGgttt                                                                                                                                                                   >  2:2123/1‑251 (MQ=255)                                                                                    gACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTcgc                                                                                                                                                                <  1:72983/251‑1 (MQ=255)                                                                                     aCGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCTGGGCGAGGCGATTAGTACTGAGTGAAGATCTGCTGGATCTGCTGCGGATCTTTTGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGATGAAACCGTATTTCGCATCATCCACTTCGGCAGCCTGAGTGGTAGCGCACGTCGGTACGCGGGAAGAACGCGCGACTGCAGTACCGGTTTTcgcg                                                                                                                                                               <  2:566755/251‑1 (MQ=255)                                                                                       gTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCgcggc                                                                                                                                                             >  1:358741/1‑251 (MQ=255)                                                                                         gAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGcggcgg                                                                                                                                                           <  1:575795/251‑1 (MQ=255)                                                                                               gAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCg                                                                                                                                                        <  1:345025/248‑1 (MQ=255)                                                                                                aaGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAg                                                                                                                                                    <  1:472457/251‑1 (MQ=255)                                                                                                      aaaaGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTc                                                                                                                                              <  1:36570/251‑1 (MQ=255)                                                                                                      aaaaGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTc                                                                                                                                              <  2:167812/251‑1 (MQ=255)                                                                                                       aaaGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGAGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCTCGAGAGTGTCg                                                                                                                                             >  2:430064/1‑251 (MQ=255)                                                                                                        aaGCCCAGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCgg                                                                                                                                            >  1:529874/1‑250 (MQ=255)                                                                                                              cGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGCTCTGCTGCGGCTCTTTGGTTTGCGTCCGCGCCCGTTGCCGCCGCCCGCGCGCTTTTTGCGGGTTCCGCGTGCCCGCGGCGCCGCCCCCGTCTTTCGCCTCCTCCCCTTCGGCCTCCTGCGTGGTCGCGCCCGTCGGTCCGCGGGCCGCCCGCCCGCCTGCCGTCCCGGTTTTCGCTGCGGTCGCCCGCgttgcgccccccg                                                                                                                                      >  1:15584/1‑240 (MQ=255)                                                                                                                    gATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTtat                                                                                                                                >  1:34648/1‑251 (MQ=255)                                                                                                                            ggcgaggcgaTTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACCGGTTg                                                                                                                        >  2:501339/1‑251 (MQ=255)                                                                                                                            ggcgaggcgaTTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTg                                                                                                                        <  1:316750/251‑1 (MQ=255)                                                                                                                            ggcgaggcgaTTAGTACTGATTGAAGATCTGCCGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTg                                                                                                                        <  1:144824/251‑1 (MQ=255)                                                                                                                              cgaggcgaTTAGTACTGATAGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGGAGAAAGCACGACTGCAGTACCGGTTTTCGAGGCGGTCGCCAGCGTGTGGCACACAGATTTTTACAGGtttcc                                                                                                                      >  2:70565/1‑248 (MQ=255)                                                                                                                                        aGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACCCAGATTTCTAAAGGTTGCCGTTACCcac                                                                                                             >  1:231613/1‑250 (MQ=255)                                                                                                                                          tACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACAcc                                                                                                          >  1:72174/1‑251 (MQ=255)                                                                                                                                            cTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGGTTCGGACACAGATTTATACAGGTTTCCCTTAcccccccccc                                                                                                        >  2:18428/1‑244 (MQ=255)                                                                                                                                             tGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAgc                                                                                                       >  1:40963/1‑251 (MQ=255)                                                                                                                                              gATTGAAGAGCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAgc                                                                                                       <  2:240275/250‑1 (MQ=255)                                                                                                                                                  gAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTGTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCTTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATGTATACAGGTTGCCGTTACCCACACCATCGCTa                                                                                                   <  2:237019/250‑1 (MQ=255)                                                                                                                                                        ctgctggatctgctgCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCCAACACAGATTTATACAGGTTTCCGTTACCCACACCAGCGCTAACGNTGc                                                                                            >  2:599262/1‑251 (MQ=255)                                                                                                                                                         tgctggatctgctgCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACAAAGATTTATACAGGTTGCCGGTACCCACACCAGCGCTAACGaaggc                                                                                           >  2:210592/1‑247 (MQ=255)                                                                                                                                                          gctggatttgcttCGGATCTTTTGTGTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTCTGCGGGTTTAGCGTGCCAGAGGCGACGAATCCGTATTTCTCAACATCCATTTCGGCTTCCTGACTGGTAGCGCCCGTCGGTACGCGGGAATAACGCACGACTGCAGTACCGGTTTTCTCGTCGGTGGATAGCGTGTCTAACGTAGAATTATACAGGTTGCATTTACCCCCACCATCGCGAACGATg                                                                                             <  2:647439/248‑1 (MQ=255)                                                                                                                                                                  tgctgcGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCCGCGTGCCGGACACAGATTTATACCGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGcc                                                                                  >  2:529844/1‑251 (MQ=255)                                                                                                                                                                        ggATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATACACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGAATGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGGGTCGCAAACCGATTTATACAGGAGGACGGTAACACCACACGCGCTAAAGATGCCATCATAGCCCGccgt                                                                            >  2:190900/1‑248 (MQ=255)                                                                                                                                                                               tGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAg                                                                      <  1:501339/250‑1 (MQ=255)                                                                                                                                                                               tGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCGTTT‑GCGGGTTCATCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGAGGCGGTCGCCAGCGTGTCGAACACAGATTTATATAGGTTGCCGTTACCCGCACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGt                                                                     <  1:415369/250‑1 (MQ=255)                                                                                                                                                                                   ttGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTCTACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAACCCGCCTCTACCAATGCtt                                                                 >  1:635062/1‑251 (MQ=255)                                                                                                                                                                                     gCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTa                                                               <  1:390863/251‑1 (MQ=255)                                                                                                                                                                                        tCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTc                                                                                                                                              <  1:238312/169‑1 (MQ=255)                                                                                                                                                                                        tCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTc                                                                                                                                              >  2:238312/1‑169 (MQ=255)                                                                                                                                                                                        tCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTa                                                               <  1:442107/248‑1 (MQ=255)                                                                                                                                                                                        tCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGGTAACCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTa                                                               >  2:442107/1‑248 (MQ=255)                                                                                                                                                                                         cAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGcc                                                            <  1:505444/250‑1 (MQ=255)                                                                                                                                                                                          agagCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACTCCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCgg                                                          <  1:257088/251‑1 (MQ=255)                                                                                                                                                                                                aGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATc                                                    <  2:137919/251‑1 (MQ=255)                                                                                                                                                                                                aGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATc                                                    <  2:318748/251‑1 (MQ=255)                                                                                                                                                                                                    gcagcaGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCgaat                                                <  1:413663/251‑4 (MQ=255)                                                                                                                                                                                                      agcagAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTTCCGTTACCCACAACAGCGCCAACGATGCCATCTAAGCCCGCATATACCAGTGCTTTAACCCGAAGATccaatgt                                              >  2:501143/1‑245 (MQ=255)                                                                                                                                                                                                      agcagAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGAGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCCGCGTTTCGAACCCAGATTTATACAGGTTGCCGTTACCCACAACCGCGCTAACGGTGCCATACTAACCCGCATCTACCAGTGCTTTAGcccgacaatcgaatgc                                              >  2:34857/1‑237 (MQ=255)                                                                                                                                                                                                        cagAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCCGCGCTCACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCgg                                            >  2:334126/1‑250 (MQ=255)                                                                                                                                                                                                            aCGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGc                                        <  2:472569/251‑1 (MQ=255)                                                                                                                                                                                                               cgcgcTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGCGCTTTAGCCGGAAGATCGAATGCGTTAGCGt                                      >  2:447547/1‑250 (MQ=255)                                                                                                                                                                                                                  gcTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGCGTAGtt                                  <  2:505202/251‑1 (MQ=255)                                                                                                                                                                                                                     ttttGCGGTTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACCCCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGCGTAGTta                                 <  1:124012/249‑1 (MQ=255)                                                                                                                                                                                                                       ttGCGGGTTCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGCGTAGTTATaa                              >  1:482637/1‑250 (MQ=255)                                                                                                                                                                                                                          cGGGTTCTGCGTGCCAGAGTCGACGAATCCGTATGTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTTTCGAACACAGATTTATTCAGGTTGCCGTTACCCACACCATCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGGTCGAATGCGTTAGCGTAGTTATAAACaa                          <  1:409278/251‑1 (MQ=255)                                                                                                                                                                                                                              ttCAGCGTGCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGCTGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGAGTAGTTATAAACAATGcc                      >  1:472074/1‑251 (MQ=255)                                                                                                                                                                                                                                cgncGTGCCAGAGGCGACGAAACCGTATTTCGCATAATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGGTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGCGTAGTTATAAACAATGCCg                     <  2:594412/247‑1 (MQ=255)                                                                                                                                                                                                                                      gCCAGAGGCGACGAAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGCGTAGTTATAAACAATGCCGACTTTc               <  2:468587/250‑1 (MQ=255)                                                                                                                                                                                                                                               gacgaAACCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTTCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATCTATACAGGTTGCCGTTACCCACGCCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGAATGCGTTAGCGTAGTTATAAACAATGCCGACTTTCGGCCGttcat     <  2:474197/251‑1 (MQ=255)                                                                                                                                                                                                                                                  gaANCCGTATTTCGCATCATCCACTTCGTCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAACATCGAATGCGTTAGCGTAGTTATAAACAATGCCGACTTTCGGCAGttcattca  <  2:587165/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        ATCCATTGCTTTTTCACCCGAAAAACCTCACGCGGAATAATTTATTCCGGTGATTATCATTAGGGTAATTAAAATTGGCAATGGACGTGAGTCTGAAGTGAAAAAGCCCCGGCACGATACCGGGGCGAGGCGATTAGTACTGATTGAAGATCTGCTGGATCTGCTGCGGATCTTTGGTTTGCGTCAGAGCCAGTTGCAGCAGAACGCGCGCTTTTTGCGGGTTCAGCGTGCCAGAGGCGACGAAGCCGTATTTCGCATCATCCACTTCGGCATCCTGAGTGGTAGCGCCCGTCGGTACGCGGGAAGAACGCACGACTGCAGTACCGGTTTTCGCGGCGGTCGCCAGCGTGTCGAACACAGATTTATACAGGTTGCCGTTACCCACACCAGCGCTAACGATGCCATCATAGCCCGCATCTACCAGTGCTTTAGCCGGAAGATCGGATGCGTTAGCGTAGTTATAAACAATGCCGACTTTCGGCAGTTCATTCA  >  NC_000913/3099549‑3100040

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