Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,845,795	C→A	5.5%	R389L (CGG→CTG)	topB ←	DNA topoisomerase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,845,795	0	C	A	5.5%	275.1 / 3.5	82	R389L (CGG→CTG)	topB	DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base C (35/42);  new base A (3/2);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TTAGCGACAAATTTGCCTTTGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAA  >  NC_000913/1845651‑1845941                                                                                                                                                 |                                                                                                                                                   ttAGCGACAAATTTGCCTTTGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTg                                                                                                                                                 >  6:82832/1‑148 (MQ=255)       aCAAATTTGCCTTTGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCg                                                                                                                                           >  1:235264/1‑148 (MQ=255)         aaaTTTGCCTTTGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCgg                                                                                                                                          >  1:309902/1‑147 (MQ=255)            tttGCCTTTGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAatg                                                                                                                                      <  1:70890/148‑1 (MQ=255)              tGCCTTTGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCTGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAatga                                                                                                                                     <  3:53503/147‑1 (MQ=255)                    tGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgt                                                                                                                               <  5:382987/147‑1 (MQ=255)                    tGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgt                                                                                                                               <  7:83536/147‑1 (MQ=255)                     ggCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgtg                                                                                                                              <  2:259818/147‑1 (MQ=255)                      gCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgtgg                                                                                                                             <  3:284420/147‑1 (MQ=255)                          tGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTgcg                                                                                                                         >  1:78420/1‑147 (MQ=255)                                gTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCg                                                                                                                                                   <  6:48147/115‑1 (MQ=255)                                gTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCg                                                                                                                                                   >  5:48147/1‑115 (MQ=255)                                gTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgt                                                                                                                               <  4:283940/135‑1 (MQ=255)                                gTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgt                                                                                                                               >  3:283940/1‑135 (MQ=255)                                 ttCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGAcc                                                                                                                 >  3:306735/1‑148 (MQ=255)                                 ttCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGAcc                                                                                                                 <  6:39786/148‑1 (MQ=255)                                  tCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCt                                                                                                                >  2:103962/1‑148 (MQ=255)                                  tCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCAGCGAccg                                                                                                                >  6:277332/1‑147 (MQ=255)                                    gATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCttt                                                                                                              <  2:309902/148‑1 (MQ=255)                                    gATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCttt                                                                                                              <  7:379713/148‑1 (MQ=255)                                    gATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCtt                                                                                                               <  2:235264/147‑1 (MQ=255)                                         cacacTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCAGTCAGGTTGATCGCAGAACTCAGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGt                                                                                                          >  2:90068/1‑147 (MQ=255)                                          acacTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACATTCGCTTCGTTCTCCGTCAGGTTGATAGCAGAACTCAGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCa                                                                                                        >  6:210409/1‑148 (MQ=255)                                             cTTGCGGACCACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATcc                                                                                                     >  5:322803/1‑148 (MQ=255)                                                gCGGAACACCGCATCCGGGCAGAATTGCATCAGATGCTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCa                                                                                                   <  7:38103/147‑1 (MQ=255)                                                 cGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCaa                                                                                                  >  6:26317/1‑147 (MQ=255)                                                    aaCACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGc                                                                                              >  5:311619/1‑148 (MQ=255)                                                     acacCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGTACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCg                                                                                             >  8:394685/1‑148 (MQ=255)                                                       acCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGGCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGt                                                                                           <  2:378848/148‑1 (MQ=255)                                                           cATCCGGGCAGAATTGCATCAGTTACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGTACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCg                                                                                       <  7:394685/148‑1 (MQ=255)                                                            aTCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGt                                                                                      <  2:78420/148‑1 (MQ=255)                                                            aTCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGt                                                                                      <  6:220747/148‑1 (MQ=255)                                                            aTCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGt                                                                                      <  8:22979/148‑1 (MQ=255)                                                               cGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGtata                                                                                   <  1:323849/148‑1 (MQ=255)                                                                 ggCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATc                                                                                 <  6:322803/148‑1 (MQ=255)                                                                 ggCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATc                                                                                 <  2:156709/148‑1 (MQ=255)                                                                      aaTTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGAt                                                                            >  8:169654/1‑148 (MQ=255)                                                                       aTTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCTGAACTCCGTGCGGTAGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATc                                                                           <  6:310093/148‑1 (MQ=255)                                                                           cATCAGATACTGACTGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTTATGGCAGTACTCCGGGCGGGCGGAATGATGGCGTGGGGCGCAGCGACCTTTTTGGCATCCCAACAGCGGTTGCGTATATCTGGATCTAcc                                                                       <  5:82832/148‑1 (MQ=255)                                                                              cAGATACTGTTGGGCTTTCAGGTTATTGACCTTCGGTTTGTTGTGCGTGAGGTTTATCGCTGATCTCCGTGCGGTCGGAATGATGTCGTGGTGCGCATGGACCTTTTTGTCATCCCAACAGGGGTTGCGtttt                                                                                   <  8:155992/133‑4 (MQ=255)                                                                               aGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGCTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTatat                                                                                  >  7:155992/1‑133 (MQ=255)                                                                                 aTACTGACGGGCAATCAGGTGATAGACCTTCGCTTCGTTATCCGTCAGGGTGATCGCAGAACTACGTGCGGTCGGAATGAAGGCGTGGTTCTCAGCGACAGTTTTTTCATCCCCCCAGCGGTGGCGGATATCTGGCTCTCCCACTGGc                                                                 >  6:15595/1‑148 (MQ=255)                                                                                       aCGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCgg                                                            <  1:38975/147‑1 (MQ=255)                                                                                       aCGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCg                                                             <  7:345656/146‑1 (MQ=255)                                                                                        cGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGc                                                           >  7:15590/1‑147 (MQ=255)                                                                                         gggCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgtggtg                                                                                                                           <  3:161080/82‑1 (MQ=255)                                                                                         gggCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCgtggtg                                                                                                                           >  4:161080/1‑82 (MQ=255)                                                                                          ggCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTa                                                                         <  2:152381/131‑1 (MQ=255)                                                                                          ggCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTa                                                                         >  1:152381/1‑131 (MQ=255)                                                                                          ggCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCAGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCaa                                                         >  1:396059/1‑147 (MQ=255)                                                                                               tCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGtata                                                                                   <  2:392342/116‑1 (MQ=255)                                                                                               tCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGtata                                                                                   >  1:392342/1‑116 (MQ=255)                                                                                                cAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCgg                                                                                            <  7:81899/106‑1 (MQ=255)                                                                                                cAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCgg                                                                                            >  8:81899/1‑106 (MQ=255)                                                                                                cAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                  >  4:97192/1‑148 (MQ=255)                                                                                                  ggTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCAGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCgg                                                <  2:396059/148‑1 (MQ=255)                                                                                                    ttATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTg                                              >  7:359015/1‑148 (MQ=255)                                                                                                     tataGACCTTCGCTTCGTACTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACATTTTTGTCATCCCAACAGCGGTTGCGTATATCTGCATCTACCACTGGCTGCGGCAACAGATCCtgtg                                              >  3:106611/1‑147 (MQ=255)                                                                                                       taGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCa                                            <  3:97192/147‑1 (MQ=255)                                                                                                       taGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCAt                                           >  2:234883/1‑148 (MQ=255)                                                                                                         gACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATg                                          >  5:111329/1‑147 (MQ=255)                                                                                                            cttcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCAGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAaca                                      <  1:90068/148‑1 (MQ=255)                                                                                                              tcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                  <  6:207715/134‑1 (MQ=255)                                                                                                              tcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                  <  5:276389/134‑1 (MQ=255)                                                                                                              tcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                  >  5:207715/1‑134 (MQ=255)                                                                                                              tcgcttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATcc                                                  >  6:276389/1‑134 (MQ=255)                                                                                                                 cttcgTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGAt                                 <  5:210409/148‑1 (MQ=255)                                                                                                                      ttCTCCGTCAGGTTGATCGCTTAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCAt                            <  5:280606/148‑1 (MQ=255)                                                                                                                      ttCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCAt                            >  2:105921/1‑148 (MQ=255)                                                                                                                      ttCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCAt                            <  1:103962/148‑1 (MQ=255)                                                                                                                       tctcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGAt                                 <  2:20824/142‑1 (MQ=255)                                                                                                                       tctcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGAt                                 >  1:20824/1‑142 (MQ=255)                                                                                                                        ctcCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCAtt                           <  8:359015/147‑1 (MQ=255)                                                                                                                          ccGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtcat                        >  8:267441/1‑148 (MQ=255)                                                                                                                           cGTGAGGTTGTTCGCATAACTGCGTGGGGTAGGATTGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCAAACAGCGGTTGCGTATATGTGGATCTACCACTGGCTGCGGGAACAGATCCGGTGCATGAACACTGATGGCATtcat                        <  4:106611/147‑1 (MQ=255)                                                                                                                           cGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtc                          <  1:200802/145‑1 (MQ=255)                                                                                                                           cGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATtc                          <  1:105921/145‑1 (MQ=255)                                                                                                                              cAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCAcc                    >  3:12853/1‑148 (MQ=255)                                                                                                                                       cGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGt           >  1:285138/1‑148 (MQ=255)                                                                                                                                         cAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTcc         <  4:306735/148‑1 (MQ=255)                                                                                                                                             aCTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGc      >  1:36135/1‑147 (MQ=255)                                                                                                                                             aCTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGGTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGc      >  3:274024/1‑147 (MQ=255)                                                                                                                                                ccGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCaaaa  >  6:141151/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                   TTAGCGACAAATTTGCCTTTGGCAATGTCCAGTTCGATAACACACTTGCGGAACACCGCATCCGGGCAGAATTGCATCAGATACTGACGGGCAATCAGGTTATAGACCTTCGCTTCGTTCTCCGTCAGGTTGATCGCAGAACTCCGTGCGGTCGGAATGATGGCGTGGTGCGCATCGACCTTTTTGTCATCCCAACAGCGGTTGCGTATATCTGGATCTACCACTGGCTGCGGCAACAGATCCGGTGCATGAACACTGATGGCATTCATCACCGCGTGGCGTCCGGCAAAA  >  NC_000913/1845651‑1845941

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