Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,879	1	.	C	96.6%	213.9 / ‑0.2	59	coding (671/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (26/31);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTC  >  NC_000913/3815743‑3816012                                                                                                                                           |                                                                                                                                       aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑ccc                                                                                                                                     >  2:2216631/1‑140 (MQ=255)   aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑ccccc                                                                                                                                   <  1:544408/140‑2 (MQ=255)    aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑ccccct                                                                                                                                  <  2:955915/140‑3 (MQ=255)    aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑ccccct                                                                                                                                  <  2:1806465/140‑3 (MQ=255)    aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑ccccct                                                                                                                                  >  1:462040/1‑138 (MQ=255)         gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGc                                                                                                                             >  2:1485443/1‑140 (MQ=255)          cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGcc                                                                                                                            >  1:1687349/1‑140 (MQ=255)             tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCaga                                                                                                                         <  1:227221/140‑1 (MQ=255)              gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCagag                                                                                                                        <  1:1485443/140‑1 (MQ=255)                 atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCca                                                                                                                     <  1:761683/140‑1 (MQ=255)                   ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCag                                                                                                                          <  2:372448/133‑1 (MQ=255)                   ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCag                                                                                                                          >  1:372448/1‑133 (MQ=255)                      ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAg                                                                                                                <  1:1061419/140‑1 (MQ=255)                      ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAg                                                                                                                <  1:537547/140‑1 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATg                                                                                                             >  1:291234/1‑140 (MQ=255)                          tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATg                                                                                                             >  2:382605/1‑140 (MQ=255)                              ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  1:1054358/140‑1 (MQ=255)                              ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                        <  1:1712913/140‑1 (MQ=255)                               tCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGc                                                                                                       <  1:750606/140‑1 (MQ=255)                                       cATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                                                                               >  1:2038162/1‑140 (MQ=255)                                           ttACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                                                                           <  2:1687349/140‑1 (MQ=255)                                           ttACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                                                                           <  2:842618/140‑1 (MQ=255)                                             aCTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGACCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                                                                                         >  1:1950924/1‑140 (MQ=255)                                                  tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                    <  2:750075/140‑1 (MQ=255)                                                  tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                    <  2:2038162/140‑1 (MQ=255)                                                    tACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                                                  >  1:87028/1‑140 (MQ=255)                                                    tACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                                                  >  1:2074639/1‑140 (MQ=255)                                                      cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                                                <  2:291234/140‑1 (MQ=255)                                                      cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                                                <  2:378238/140‑1 (MQ=255)                                                      cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                                                >  2:501793/1‑140 (MQ=255)                                                      cAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                                                >  2:922755/1‑140 (MQ=255)                                                       aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                                                                               >  2:835603/1‑140 (MQ=255)                                                              aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                         <  1:382605/140‑1 (MQ=255)                                                                   aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                    <  2:137044/140‑1 (MQ=255)                                                                    ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGACCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                   <  2:1950924/140‑1 (MQ=255)                                                                      cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                 >  1:1840649/1‑140 (MQ=255)                                                                      cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                 <  2:2009925/140‑1 (MQ=255)                                                                       gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGc                                                                >  1:1751239/1‑140 (MQ=255)                                                                        aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCACGCCTTCTGCCGTCCCCTGCa                                                               >  2:2176213/1‑140 (MQ=255)                                                                        aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TAGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa                                                               >  1:568814/1‑140 (MQ=255)                                                                         cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCAc                                                              >  2:2075143/1‑140 (MQ=255)                                                                                    ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CCCCATCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATgc                                                   >  2:1810398/1‑140 (MQ=255)                                                                                      ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                 >  2:780733/1‑124 (MQ=255)                                                                                      ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                 <  1:780733/124‑1 (MQ=255)                                                                                        aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGccc                                               >  2:1875977/1‑140 (MQ=255)                                                                                          aaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TAGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGt                                             <  2:568814/140‑1 (MQ=255)                                                                                               aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                        >  1:898556/1‑140 (MQ=255)                                                                                               aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                        >  1:2210824/1‑140 (MQ=255)                                                                                                gTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                       <  1:835603/140‑1 (MQ=255)                                                                                                gTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                       <  2:2074639/140‑1 (MQ=255)                                                                                                   tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                    >  1:1705885/1‑140 (MQ=255)                                                                                                    gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                   <  1:1164294/140‑1 (MQ=255)                                                                                                      cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                 <  1:2075143/140‑1 (MQ=255)                                                                                                      cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CCCCATCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                 <  1:1810398/140‑1 (MQ=255)                                                                                                       aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                <  2:87028/140‑1 (MQ=255)                                                                                                         cgccgcCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                 >  1:1977933/1‑137 (MQ=255)                                                                                                         cgccgcCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatca                                 <  2:1977933/137‑1 (MQ=255)                                                                                                         cgccgcCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                   >  1:574329/1‑135 (MQ=255)                                                                                                         cgccgcCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                   <  2:574329/135‑1 (MQ=255)                                                                                                         cgccgcCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                    >  1:729592/1‑134 (MQ=255)                                                                                                         cgccgcCTTCTGCGTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                    <  2:729592/134‑1 (MQ=255)                                                                                                                      gTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                 <  2:1705885/140‑1 (MQ=255)                                                                                                                      gTCGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                 <  1:454253/140‑1 (MQ=255)                                                                                                                        cGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                           >  2:1858603/1‑128 (MQ=255)                                                                                                                        cGCTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                           <  1:1858603/128‑1 (MQ=255)                                                                                                                          cTACAATGGATTCGATTCCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct             >  1:71738/1‑140 (MQ=255)                                                                                                                                  gattcgattcCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga     <  1:1875977/140‑1 (MQ=255)                                                                                                                                    ttcgattcCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGt   >  2:152951/1‑140 (MQ=255)                                                                                                                                     tcgattcCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc  >  1:1511758/1‑140 (MQ=255)                                                                                                                                     tcgattcCCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc  >  1:812454/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                       ATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTC  >  NC_000913/3815743‑3816012

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