Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,846,254	A→T	5.1%	F236Y (TTT→TAT) ‡	topB ←	DNA topoisomerase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,846,254	0	A	T	5.1%	251.2 / 2.9	78	F236Y (TTT→TAT) ‡	topB	DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base T (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGC  >  NC_000913/1846105‑1846402                                                                                                                                                      |                                                                                                                                                     cTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTaa                                                                                                                                                      >  2:1011278/1‑150 (MQ=255)   gATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCTGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAAc                                                                                                                                                    >  1:1040140/1‑150 (MQ=255)       cccGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCt                                                                                                                                                >  2:1287352/1‑150 (MQ=255)           tttATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATc                                                                                                                                            <  1:81345/150‑1 (MQ=255)             tATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTTTACCGCtatctgg                                                                                                                                          <  1:1197487/150‑7 (MQ=255)              aTCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTTTACCGCtatctggc                                                                                                                                         <  1:611083/150‑8 (MQ=255)               tCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCa                                                                                                                                        >  2:1240930/1‑150 (MQ=255)                 gTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAgg                                                                                                                                      >  2:803442/1‑150 (MQ=255)                 gTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAgg                                                                                                                                      <  2:1142873/150‑1 (MQ=255)                          ggTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGat                                                                                                                             >  1:1090257/1‑150 (MQ=255)                          ggTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGat                                                                                                                             >  1:1233881/1‑150 (MQ=255)                          ggTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACAGGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGat                                                                                                                             >  2:1034564/1‑150 (MQ=255)                               cAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATAtgtg                                                                                                                        >  1:1615485/1‑150 (MQ=255)                                     ccGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGa                                                                                                                  >  2:801475/1‑150 (MQ=255)                                       ggTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACt                                                                                                                <  2:1236606/150‑1 (MQ=255)                                          tGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCa                                                                                                             <  1:492332/150‑1 (MQ=255)                                           gACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCaa                                                                                                            <  2:326825/150‑1 (MQ=255)                                            aCCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCaaa                                                                                                           >  1:355895/1‑150 (MQ=255)                                              cACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAaaga                                                                                                         <  2:1446758/150‑1 (MQ=255)                                              cACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAaaga                                                                                                         <  2:1440971/150‑1 (MQ=255)                                                 tAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGt                                                                                                      <  2:42506/150‑1 (MQ=255)                                                 tAATGCGGTTAACCACATGCTCCCCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGt                                                                                                      <  1:1048971/150‑1 (MQ=255)                                                       ggTTAACCACATGCTCTGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCg                                                                                                <  2:1040140/150‑1 (MQ=255)                                                       ggTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCg                                                                                                <  1:1318665/150‑1 (MQ=255)                                                            aCCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACg                                                                                           >  1:444538/1‑150 (MQ=255)                                                            aCCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACg                                                                                           >  2:797920/1‑150 (MQ=255)                                                             ccACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGa                                                                                          >  2:214718/1‑150 (MQ=255)                                                              cacaTGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGaa                                                                                         <  2:1333343/150‑1 (MQ=255)                                                              cacaTGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGaa                                                                                         >  1:424231/1‑150 (MQ=255)                                                               acaTGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAg                                                                                        >  1:212185/1‑150 (MQ=255)                                                               acaTGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAg                                                                                        >  2:1385349/1‑150 (MQ=255)                                                                  tGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGttt                                                                                     >  2:82953/1‑150 (MQ=255)                                                                    cTCCGCCAGTTGACGATGTAACAAGCGACCTTCTTCTTCCTGGTACGGTTTACTCGCTTCGCTCGGTTGCCATATAGCGGTAAACCGCTCATCGGCAGTTGTCACGATATGTGCTTTGACTTCTAAGAAGTCTTTCGCCACGAAGTTTTc                                                                                   <  2:1606641/150‑1 (MQ=255)                                                                    cTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTc                                                                                   >  2:818954/1‑150 (MQ=255)                                                                    cTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTc                                                                                   >  1:614278/1‑150 (MQ=255)                                                                    cTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTc                                                                                   <  2:808953/150‑1 (MQ=255)                                                                    cTCCGACAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCCCACGCTTCGCTCCGTTGCCAGATAGCGCTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTc                                                                                   >  2:144614/1‑150 (MQ=255)                                                                      ccGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCaa                                                                                 <  1:1423956/150‑1 (MQ=255)                                                                         ccAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAAtct                                                                              <  1:1240930/150‑1 (MQ=255)                                                                          cAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAAtctc                                                                             <  1:1011278/150‑1 (MQ=255)                                                                             tGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTc                                                                          <  1:1467602/150‑1 (MQ=255)                                                                                 cGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATcg                                                                      >  2:66960/1‑150 (MQ=255)                                                                                 cGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATcg                                                                      >  1:1305326/1‑150 (MQ=255)                                                                                 cGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGATTTCAATCTCTTCATcg                                                                      >  2:699692/1‑150 (MQ=255)                                                                                  gATTTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGATTTGACTTCAAAGAAGTCTTTCGCCACTAAGATTTCAATCTCTTCATcgc                                                                     >  1:847819/1‑150 (MQ=255)                                                                                     gTAACAAGCGCCCTTCTTCATCCTGGTACTGTTCAGACGCTTCGCTCGGTTGCCAGATAGGGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCgcggc                                                                  <  2:1006723/150‑1 (MQ=255)                                                                                       aaCAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGgcgc                                                                <  2:77966/150‑1 (MQ=255)                                                                                        aCAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAgg                                                                                                                                      >  1:1244043/1‑79 (MQ=255)                                                                                        aCAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAgg                                                                                                                                      <  2:1244043/79‑1 (MQ=255)                                                                                          aaGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGcacc                                                             >  2:115203/1‑150 (MQ=255)                                                                                           aGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGcacca                                                            >  2:383113/1‑150 (MQ=255)                                                                                            gcgcCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGcaccac                                                           <  1:433784/150‑1 (MQ=255)                                                                                             cgcCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCaccacc                                                          >  2:827199/1‑150 (MQ=255)                                                                                             cgcCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCaccacc                                                          >  1:1045101/1‑150 (MQ=255)                                                                                              gcCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCAccacca                                                         <  1:246643/150‑1 (MQ=255)                                                                                                ccTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGc                                                       >  1:145680/1‑150 (MQ=255)                                                                                                  ttcttcATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGccc                                                     <  1:818954/150‑1 (MQ=255)                                                                                                   tcttcATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCa                                                    >  2:626467/1‑150 (MQ=255)                                                                                                    cttcATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTc                                                                                                 <  1:670731/104‑1 (MQ=255)                                                                                                    cttcATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTc                                                                                                 >  2:670731/1‑104 (MQ=255)                                                                                                         tCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCTGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACg                                              >  1:114239/1‑150 (MQ=255)                                                                                                           cTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACggg                                            >  1:1126476/1‑150 (MQ=255)                                                                                                           cTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACggg                                            >  1:448304/1‑150 (MQ=255)                                                                                                            tGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGc                                           >  2:1382832/1‑150 (MQ=255)                                                                                                             ggTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCg                                          <  1:214718/150‑1 (MQ=255)                                                                                                               tACGGTTCACACGCTTCGGTCGGTTGCCTGATAGCGGTAATCCGCTCATCGGGAGTTGTCACTATATGTGCTTTGACTTCAAAGAAGTCTTTCGGCTCGAAGTTTTCAATCTCTTCATCGCGGGGGAGCACCAGCCCATGCACGGGCGTc                                        <  2:1071605/150‑1 (MQ=255)                                                                                                                  ggTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGAGGCGAACCACCAGCCCAAGCACGGGCGTCTGc                                     >  2:715289/1‑150 (MQ=255)                                                                                                                         acGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTc                              <  1:1287352/150‑1 (MQ=255)                                                                                                                          cGCTTCGCTCGGTTTCCACATATCGGTAAACCGGCCATTGGCATGGGTCACGTTTTGTTCTTTGTCTTTAAACAAGTCCTTCGGCACGAATTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTcc                             <  2:1608016/150‑1 (MQ=255)                                                                                                                            cTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCa                           <  2:1126476/150‑1 (MQ=255)                                                                                                                                    ggTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGt                   <  2:588314/150‑1 (MQ=255)                                                                                                                                       tGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACg                >  2:1349945/1‑150 (MQ=255)                                                                                                                                         ccgatgaGCGGTTTACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGcc              >  1:1105341/7‑150 (MQ=255)                                                                                                                                         ccAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGcc              <  1:1385349/150‑1 (MQ=255)                                                                                                                                         ccAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGcc              <  1:383113/150‑1 (MQ=255)                                                                                                                                           gatgaGCGGTTTACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCt            >  1:351823/5‑150 (MQ=255)                                                                                                                                                  ggTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAAcc     <  1:1382832/150‑1 (MQ=255)                                                                                                                                                   gTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCg    >  1:1110760/1‑150 (MQ=255)                                                                                                                                                   gTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCg    <  1:1349945/150‑1 (MQ=255)                                                                                                                                                     aaaCCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCCCGGAAAGTACGCCCTGATAACCGGc  >  1:838773/1‑150 (MQ=255)                                                                                                                                                     aaaCCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGc  <  2:1161739/150‑1 (MQ=255)                                                                                                                                                     aaaCCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGc  >  1:9485/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CTGATTCCCGTTTATCGTTATAGCTGGTGACAATAGCCGGTTGACCACTAATGCGGTTAACCACATGCTCCGCCAGTGGACGATGTAACAAGCGCCCTTCTTCATCCTGGTACGGTTCACACGCTTCGCTCGGTTGCCAGATAGCGGTAAACCGCTCATCGGCAGGTGTCACGATATGTGCTTTGACTTCAAAGAAGTCTTTCGCCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGCACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGC  >  NC_000913/1846105‑1846402

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