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	94.8%	191.6 / 2.9	58	coding (671/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (24/31);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGA  >  NC_000913/3815743‑3816009                                                                                                                                            |                                                                                                                                    aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCTAGCACCGCCTTCTGCGCCGCTACAATGGATTCGA‑T‑T‑c‑cc                                                                                                                                 >  2:591517/1‑140 (MQ=255)    aaCTGGCGCTGATATTGTTTCATGCTTTCCCTCCTCACCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑c‑cccct                                                                                                                              <  1:591517/140‑3 (MQ=255)         gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGc                                                                                                                          <  2:86792/140‑1 (MQ=255)         gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGc                                                                                                                          <  2:391170/140‑1 (MQ=255)          cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGcc                                                                                                                         >  2:1072000/1‑140 (MQ=255)            cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCag                                                                                                                       >  2:158158/1‑140 (MQ=255)                tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGcc                                                                                                                   <  2:624072/140‑1 (MQ=255)                      ttCATGCCTTCGCTCTTCATCTCACCTTTCTACAGACCAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAg                                                                                                             <  1:434206/140‑1 (MQ=255)                      ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAg                                                                                                             <  1:32457/140‑1 (MQ=255)                      ttCATGCCTTCGCTCCTCATCTCACTTTTCTTCAGCCAAAAAAAAGGCCACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAg                                                                                                             <  2:590517/140‑1 (MQ=255)                       tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGa                                                                                                            >  2:1085975/1‑140 (MQ=255)                        cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGAt                                                                                                           <  1:1085975/140‑1 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGT‑T‑T‑CGCCCTCGGGCCAGAGCCAACAAGATg                                                                                                          >  1:513465/1‑140 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATg                                                                                                          >  2:146772/1‑140 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATg                                                                                                          >  1:635108/1‑140 (MQ=255)                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATg                                                                                                          >  1:793031/1‑140 (MQ=255)                              ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                     <  1:1542421/140‑1 (MQ=255)                                     ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                                                                              >  1:950293/1‑140 (MQ=255)                                          cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                     <  2:791784/128‑1 (MQ=255)                                          cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                     >  1:791784/1‑128 (MQ=255)                                                 ttCTATAGTCAAGAAAAAGGCGACTCATCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                                                  <  2:1424474/140‑1 (MQ=255)                                                 ttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTCCCAGCCCCGCCTTCTGCGGCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                                                  >  2:1017277/1‑140 (MQ=255)                                                  tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                 >  2:230483/1‑140 (MQ=255)                                                  tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                 >  1:77694/1‑140 (MQ=255)                                                  tCTACAGACAAAAAAAACGCGACTCATCAGTGGCCTTAACAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                 <  2:513465/140‑1 (MQ=255)                                                   cTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                >  1:172197/1‑140 (MQ=255)                                                    tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                                                               >  1:1424474/1‑140 (MQ=255)                                                      cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGGCAAGAAGTTGAGTAGCCCTTCATGGGTGAACGGCTCGcc                                                                             >  2:596538/1‑140 (MQ=255)                                                        gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                           <  2:838969/140‑1 (MQ=255)                                                        gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                           <  2:370890/140‑1 (MQ=255)                                                                aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                                                                   <  2:635108/140‑1 (MQ=255)                                                                  aGGCGACTCATCAGTGGCCTTACAACTCACTTTGCCATCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                 <  1:743438/140‑1 (MQ=255)                                                                  aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCTCCTTCTGCCGTccc                                                                 <  2:950293/140‑1 (MQ=255)                                                                     cGACTCATCAGTCGCTTTAAAAGTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                              <  1:596538/140‑1 (MQ=255)                                                                     cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCGACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                              <  1:29700/140‑1 (MQ=255)                                                                     cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACACGATGAGTATCTCGTTACGGGTGAACGGCTCGACTTCTGCAGTCCCCTg                                                              >  2:743438/1‑140 (MQ=255)                                                                     cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                              <  2:138695/140‑1 (MQ=255)                                                                                  gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                              >  2:860413/1‑127 (MQ=255)                                                                                  gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                              <  1:860413/127‑1 (MQ=255)                                                                                  gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                 >  2:1294619/1‑140 (MQ=255)                                                                                   ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGGACGCCTCGCCTTCTGCAGTCCTCTGCACTTCAATGATgc                                                >  2:397254/1‑140 (MQ=255)                                                                                   ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTACATGGGTGCACGGCTCCCCTTCTGCCGTCCCCTGCACTTCAatgatgt                                                >  2:629934/1‑139 (MQ=255)                                                                                     ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                <  2:244394/122‑1 (MQ=255)                                                                                     ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                                                >  1:244394/1‑122 (MQ=255)                                                                                     ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc                                              >  2:315911/1‑140 (MQ=255)                                                                                     ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGACGGCCCCTGAACTTCAATGATgcgc                                              >  1:399311/1‑140 (MQ=255)                                                                                            tCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCtt                                       <  2:172197/140‑1 (MQ=255)                                                                                            tCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCtt                                       <  2:77694/140‑1 (MQ=255)                                                                                              aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCg                                     <  1:1294619/140‑1 (MQ=255)                                                                                                tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                   <  1:158158/140‑1 (MQ=255)                                                                                                 ttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtc                                  <  1:315911/140‑1 (MQ=255)                                                                                                  tGCCAGCGCCGCCTTCTTCGTCTCTACAATGCATTCGA‑TAC‑C‑CCCTCGTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                 <  2:399311/140‑1 (MQ=255)                                                                                                   gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                <  2:793031/140‑1 (MQ=255)                                                                                                      aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                             <  2:833680/140‑1 (MQ=255)                                                                                                          ccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCAGCCGTCCCCTGCCCTTGAATGATGCGCCCGTGTTCGGTCATCACTACg                         >  2:756494/1‑140 (MQ=255)                                                                                                           cgccTTCTGCGTCGCTTCAATGCATTCGATT‑T‑C‑TCCTCGCGCCAGAGCGAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                        <  2:746632/140‑1 (MQ=255)                                                                                                                   gCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGGGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGTGCCCGTCTTCGGTCATCACTACGTTCATGTCg                <  1:629934/140‑1 (MQ=255)                                                                                                                   gCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                <  1:357284/140‑1 (MQ=255)                                                                                                                   gCGTCGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                <  1:151777/140‑1 (MQ=255)                                                                                                                       cGCTACAATGGATTCGA‑T‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGCAGCTCTTCATGGGTGAACGGCTCGCCTTCTGTCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct            >  2:1477760/1‑140 (MQ=255)                                                                                                                                 gattcga‑t‑tcC‑CCCTTGGGCCACAGTCAAGAAGATGATTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga  <  1:756494/140‑1 (MQ=255)                                                                                                                                   ttcga‑t‑tcC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa    >  1:222518/1‑136 (MQ=255)                                                                                                                                   ttcga‑t‑tcC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCa    <  2:222518/136‑1 (MQ=255)                                                                                                                                            |                                                                                                                                    ATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑C‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGA  >  NC_000913/3815743‑3816009

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