Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,845,901	G→A	5.2%	H354Y (CAT→TAT)	topB ←	DNA topoisomerase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,845,901	0	G	A	5.2%	244.9 / 2.4	77	H354Y (CAT→TAT)	topB	DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base A (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

TAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAA  >  NC_000913/1845753‑1846049                                                                                                                                                     |                                                                                                                                                     taGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATg                                                                                                                                                      <  4:30862/149‑1 (MQ=255)      cttcgcttcgTTCTCCGTCAGGTTGATCGCATAAATCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGCTCCGGTGCATGAacac                                                                                                                                                 >  3:266491/1‑149 (MQ=255)       ttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACt                                                                                                                                                >  4:345603/1‑149 (MQ=255)         cgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACt                                                                                                                                                <  2:271775/147‑1 (MQ=255)           cttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATg                                                                                                                                            <  8:112934/149‑1 (MQ=255)               gTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGTAATGATTGAGTGGTGCGCATCGAACTTTTTTTCATACCAACAGAGGTTGCGTATATCTGGAGCTACCACTGCGTGCGGCAACAGATCCGTGGCATGCACACTGATGGcct                                                                                                                                        >  6:236318/1‑147 (MQ=255)                  ctcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtc                                                                                                                                      <  2:142044/148‑1 (MQ=255)                   tcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtca                                                                                                                                     >  8:415957/1‑148 (MQ=255)                           gTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGa                                                                                                                                              >  5:211700/1‑131 (MQ=255)                           gTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGa                                                                                                                                              <  6:211700/131‑1 (MQ=255)                            ttGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTg                                                                                                                           <  5:307731/149‑1 (MQ=255)                             tGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTgg                                                                                                                          >  4:240152/1‑149 (MQ=255)                                 cGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTc                                                                                                                      >  1:362783/1‑149 (MQ=255)                                   cAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCg                                                                                                                    <  7:166701/149‑1 (MQ=255)                                   cAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCg                                                                                                                    <  4:171659/149‑1 (MQ=255)                                     gAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGc                                                                                                                  <  7:75604/149‑1 (MQ=255)                                         tCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa                                                                                                              >  3:17082/1‑149 (MQ=255)                                           cGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAAt                                                                                                             <  7:174901/148‑1 (MQ=255)                                                 gTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGttcttc                                                                                                      <  8:54083/149‑1 (MQ=255)                                                 gTCGGAATGATGGCGTGGTGCGCATCGAACTTTTTGTCATCCCAACATCTGTTTAGTATATCTGGATATAACACTGGCTGCGGCAACAGATCCGGTTCATGAACACTGATGGCATTCATCTACTCGTGGCTTCATGCAAAATGttctta                                                                                                      <  8:102222/149‑2 (MQ=255)                                                   cGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTg                                                                                                    >  7:147438/1‑149 (MQ=255)                                                     gAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTgg                                                                                                   >  7:131259/1‑148 (MQ=255)                                                             gCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGc                                                                                           <  8:93085/148‑1 (MQ=255)                                                              cGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         <  8:153916/149‑1 (MQ=255)                                                               gtggtgCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         >  7:75081/1‑148 (MQ=255)                                                               ggggggcgcATCGACCTTTTTGTCATCCCAACAGCGGTAGCGTATATCTGGAGCTCCCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc                                                                                        <  5:328899/144‑1 (MQ=255)                                                                   tgCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTTCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATc                                                                                    <  3:329855/149‑1 (MQ=255)                                                                    gcgcATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAAt                                                                                     >  1:281830/1‑147 (MQ=255)                                                                     cgcATCGAACTTTTTTTCATACCAACAGCGGTTGCGTATATCTGGATCTACCACTGTCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCAATCATCACCGCGTGGCGTCCTGCAAAATGTTCTTCTGGCAAATAGCGACAATCa                                                                                   <  4:266491/148‑1 (MQ=255)                                                                       cATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa                                                                                                              <  2:281473/119‑1 (MQ=255)                                                                       cATCGACCTTTTTGTCATCCCAACAGCGGTTGCGAATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa                                                                                                              >  1:281473/1‑119 (MQ=255)                                                                              cTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGAt                                                                         <  8:108477/149‑1 (MQ=255)                                                                              cTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGAt                                                                         <  3:345603/149‑1 (MQ=255)                                                                                   tGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGt                                                                     <  7:2004/148‑1 (MQ=255)                                                                                    gTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGa                                                                   >  5:341569/1‑149 (MQ=255)                                                                                     tCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGa                                                                   <  8:131259/148‑1 (MQ=255)                                                                                        tCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAg                                                               >  3:278222/1‑149 (MQ=255)                                                                                         cccAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGc                                                              >  7:401076/1‑149 (MQ=255)                                                                                           cAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatc                                                                                                                                   >  8:319796/1‑78 (MQ=255)                                                                                           cAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcatc                                                                                                                                   <  7:319796/78‑1 (MQ=255)                                                                                           cAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCtt                                                            >  6:180052/1‑149 (MQ=255)                                                                                            aaCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGGACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTtg                                                           >  7:381765/1‑149 (MQ=255)                                                                                            aaCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         >  1:138742/1‑119 (MQ=255)                                                                                            aaCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGa                                                                                         <  2:138742/119‑1 (MQ=255)                                                                                             aCAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTtgt                                                          <  8:401076/149‑1 (MQ=255)                                                                                             aCAGCGGTTGCGTAAATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTtgt                                                          <  6:327586/149‑1 (MQ=255)                                                                                               aGCGGTTGCGAATATCTGGATCTACCACTGGCTGCGGCAACAGATCAGGTGCATGAACACTTATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGAGACAATCAGAACGCGGATAAGTGATTAGCTTGTGc                                                        <  6:297109/149‑1 (MQ=255)                                                                                                    ttGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTc                                                   <  7:415957/149‑1 (MQ=255)                                                                                                      gCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCg                                                  <  2:362783/148‑1 (MQ=255)                                                                                                      gCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCg                                                  <  8:75081/148‑1 (MQ=255)                                                                                                       cGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTa                                                <  3:240152/149‑1 (MQ=255)                                                                                                         tataTCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTa                                                >  7:55939/1‑147 (MQ=255)                                                                                                             tCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGttt                                          >  4:360517/1‑149 (MQ=255)                                                                                                                 gATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTgg                                      >  4:83143/1‑149 (MQ=255)                                                                                                                     tACCACTGTCTGCGGCAACAGATCCGGTGCATTAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATATCGACAATCATAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGa                                  <  5:236318/149‑1 (MQ=255)                                                                                                                     tACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGc                                                              <  1:409977/121‑1 (MQ=255)                                                                                                                     tACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGc                                                              >  2:409977/1‑121 (MQ=255)                                                                                                                       ccACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGata                                >  8:220718/1‑149 (MQ=255)                                                                                                                       ccACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGat                                 <  6:341569/148‑1 (MQ=255)                                                                                                                       ccACTGGCGGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGa                                  <  1:314269/147‑1 (MQ=255)                                                                                                                        cACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGatat                               >  3:151187/1‑149 (MQ=255)                                                                                                                         aCTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGatat                               >  2:7091/1‑148 (MQ=255)                                                                                                                             gCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAg                           >  2:23968/1‑148 (MQ=255)                                                                                                                              cTGCGGCAACAGATCCGGTGCATGGACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCa                         <  8:381765/149‑1 (MQ=255)                                                                                                                              cTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCa                         <  7:330841/149‑1 (MQ=255)                                                                                                                               tGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCAc                        >  4:27459/1‑149 (MQ=255)                                                                                                                                     aaCAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTg                  >  2:62991/1‑149 (MQ=255)                                                                                                                                      aCAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTgc                 <  3:27459/149‑1 (MQ=255)                                                                                                                                        aGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTgcgc               >  4:139080/1‑149 (MQ=255)                                                                                                                                          aTCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGat                                 <  4:243902/129‑1 (MQ=255)                                                                                                                                          aTCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGat                                 >  3:243902/1‑129 (MQ=255)                                                                                                                                          aTCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCAc             >  3:192383/1‑149 (MQ=255)                                                                                                                                           tCCGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACt            >  7:163899/1‑149 (MQ=255)                                                                                                                                            ccGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc                                                                                        >  7:372448/1‑72 (MQ=255)                                                                                                                                            ccGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGAc                                                                                        <  8:372448/72‑1 (MQ=255)                                                                                                                                             cGGTGCATAAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCa          <  8:163899/149‑1 (MQ=255)                                                                                                                                                tGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGAc       <  1:7091/149‑1 (MQ=255)                                                                                                                                                   aTGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCa     >  8:23813/1‑148 (MQ=255)                                                                                                                                                     gAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCaaaa  >  2:88528/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAATGTTCTTCTGGCAAATAGCGACAATCAGAACGCGGATAAGTGATTAGCTTGTGCGTTTCGTACAGTTTCTGGCAGATATCAAGCACGTTCTGCGCACTCAGACCAAAA  >  NC_000913/1845753‑1846049

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