Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,143	T→C	10.2%	intergenic (‑10/+53)	argZ ← / ← argY	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,143	0	T	C	10.2%	150.1 / 5.3	67	intergenic (‑10/+53)	argZ/argY	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base C (5/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACT  >  NC_000913/2817999‑2818287                                                                                                                                                 |                                                                                                                                                 ccTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttac                                                                                                                                              >  1:1077546/1‑148 (MQ=255)  cTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttatt                                                                                                                                             >  2:578254/1‑149 (MQ=255)  cTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCAGAGAATGTGGTGCATCCGGGAGGATTCGAACCTCCGCCCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatttact                                                                                                                                             >  2:948245/1‑147 (MQ=255)     cGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCGGAGCTCCGGATGCAGCGGGGAACTTACTTTa                                                                                                                                          >  2:1427292/1‑149 (MQ=14)     cGGTCAATTCCGTGACCTTTTCGATACCTTCAGTGCTGAATATTTCGGAGTATATGGGGCATCCGGGAGGATGCGAACCTCCGACCGCTCGGTTCGGAGCCGCTTTCTCTATCCAGCTGAGTTACGGGTGCATCGGGACttttatttta                                                                                                                                          <  2:531351/149‑1 (MQ=14)      ggTTAATTCCGGCACATATTCGATCCCTTCATAGCGGAAGATCTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTAc                                                                                                                                         <  2:941130/149‑1 (MQ=14)       gTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCCCCCGAGTCCGAGCCGGGTACTCTATCCAGCTGAGCGACGGCTGCAGCGGGAAACTTAGTTTAcg                                                                                                                                        >  2:1015292/1‑148 (MQ=14)           acgcggCTAGCCCTTCGATGCCTTCATCGCTGAAGGTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATTCATCGGGAAATTTATTTTActgctg                                                                                                                                    <  1:872118/144‑1 (MQ=14)           aTTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTTCATCCGAGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATTCATCGCGAAATTTATTTTActgctg                                                                                                                                    >  1:967441/1‑149 (MQ=255)             tCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGGACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGat                                                                                                                                  >  2:730863/1‑149 (MQ=255)               cGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGattt                                                                                                                                >  2:1441509/1‑147 (MQ=255)               cGCTACCTCTTCAATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGattt                                                                                                                                <  2:1047972/149‑3 (MQ=14)                gCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATAtt                                                                                                                               >  1:973342/1‑149 (MQ=255)                    cctcctCGAGACCTTCATTGCTGCGGCTGTCGGAGAATTTGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCTGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATa                                                                                                                           <  2:951338/149‑1 (MQ=14)                     ctctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACGGCTCGGTTCGTAGCCGAGTACTCTCTCCAGCTGAGCTACGGATGCGTGGGGAAACTGATGTTACTGCAGATATTGTTAc                                                                                                                          >  2:80513/1‑149 (MQ=14)                     ctctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAc                                                                                                                          >  1:575180/1‑149 (MQ=255)                      tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATAttttgacc                                                                                                                         >  2:700266/1‑143 (MQ=255)                      tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAcc                                                                                                                         >  1:1331525/1‑149 (MQ=255)                        ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATTTTTATATCGc                                                                                                                       >  1:1094920/1‑149 (MQ=255)                          cGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                                                                                                                     <  1:795451/149‑1 (MQ=255)                            aTACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGCAACTCCGACCGCTCGGTTCGTCGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAGATTTATTTTGGATCTGACATTGACAACgctgct                                                                                                                   >  1:300679/1‑149 (MQ=14)                               ccTTCATTTCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                                                                                                                <  1:731117/149‑1 (MQ=255)                               ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                                                                                                                <  2:937725/149‑1 (MQ=255)                                 ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGAAGCAGCGGGAAATTTATTTTACTGCTGATATTGTTACCGCTCCTAAAGc                                                                                                              >  1:693688/1‑149 (MQ=255)                                   cATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCAGAGTCCTCAATCCAGCTAAGCTACGGATGCATCGGGAAAATTACTGTACTGCCGATCTTTCTACCCCCCCTAAAGccg                                                                                                            >  2:1246356/1‑148 (MQ=11)                                     ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCata                                                                                                          <  2:1300242/149‑1 (MQ=255)                                     ctGCTGACGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCata                                                                                                          <  1:243102/148‑1 (MQ=255)                                        cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCCTCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATct                                                                                                       >  1:467902/1‑148 (MQ=255)                                         tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATCGATACCGCTACTAAAGCCATATCaa                                                                                                      <  1:814465/149‑1 (MQ=255)                                          gAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg                                                                                                     <  1:578254/149‑1 (MQ=255)                                            aGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTa                                                                                                   <  1:165315/149‑1 (MQ=255)                                                   ggAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATg                                                                                            <  1:7605/149‑1 (MQ=255)                                                   ggAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATg                                                                                            <  1:419470/149‑1 (MQ=255)                                                     agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGAGGGt                                                                                          >  1:1158717/1‑149 (MQ=255)                                                      gaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                         >  1:1078476/1‑149 (MQ=255)                                                      gaATATGCTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                         <  1:1134149/149‑1 (MQ=255)                                                               ccATCCGCGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGgga                                                                                <  1:751912/148‑1 (MQ=255)                                                                cATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTGCGGATGCCTCGGGGAACTTAGTTTACTTCAGATATTGATACAGCGCCTAAAGACACATCTAGAAAGAGAGGGGGCATCCGggag                                                                               >  2:1389269/1‑149 (MQ=11)                                                                    cGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCGACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGAtt                                                                           <  2:1326863/149‑1 (MQ=255)                                                                    cGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGCACTCTATCCAGCTGAGCTACGGATGCAACGGGAAATTTATTTTCCTGCTTATATTGATACCGCTACTAAAGCCCTATCAAGTAAGAGATGGTGCATCCGGGAGGAtt                                                                           <  1:1066519/149‑1 (MQ=14)                                                                     gggAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAattta                                                                                                                                               >  2:387920/1‑80 (MQ=34)                                                                     gggAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAattta                                                                                                                                               <  1:387920/80‑1 (MQ=34)                                                                          gaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTa                                                                                                   <  1:530084/119‑1 (MQ=255)                                                                          gaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTa                                                                                                   >  2:530084/1‑119 (MQ=255)                                                                               gAACCTCCGTCCGCTTGGGTCGTAGCCGAGGACTCTCTCCAGCTGTGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                <  2:575180/149‑1 (MQ=17)                                                                               gAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                <  1:103309/149‑1 (MQ=255)                                                                               gAACCTCCGACCGCTCGGTTCGTAGCCGAGTACCCTATCCAGCTGAGCTACGGCTGCATCGGGAAATTTATTTTACTGCTGATATTGAGATACCTATTATAAAGATATAAAGAAAGAGATGGTGAGTCCGGGAGGATGCGGACTTCGGa                                                                >  2:798765/1‑149 (MQ=14)                                                                               gAACCGCCGACCGCTCGGTTCGTAGCCCAGTACTCTACCCAGCTGAGCTACGGATGCATCGGAAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                <  1:674230/149‑1 (MQ=255)                                                                                  ccTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCg                                                             >  1:1178323/1‑149 (MQ=255)                                                                                    cccGACCGCGCGGTCCGTAGCCCCTTACCCTTCGCACTTCAGTTACGGATGCAGCGGGAAATGTCATTTACGGCAGGTGTTGAGTCCGTTACTAAATCCAGATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCt                                                           <  2:408116/148‑1 (MQ=12)                                                                                        aCCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGt                                                       <  1:746026/149‑1 (MQ=255)                                                                                          cGCTCGGTTCGTAGCCGCGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGCGTTCGAACCTCCGACCGCTCGGTTc                                                     <  1:870443/149‑1 (MQ=255)                                                                                            cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGCGCTACGGATGCATCGGGAACTTTATTTTACTGCTGATTTTGATACCGCTCCTACACCCCTATCAGAAAAAAGAGGGTGCGTCCGGGAGGATACAACCTTCCGACCCCCCggcccct                                                   >  2:28216/1‑144 (MQ=9)                                                                                            cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAAAGATGGTGCGTCCGGGAGGATTCGCACCTCCGACCGCTCGGTTCGt                                                   >  1:1139276/1‑149 (MQ=255)                                                                                                gtccggaGCCGAGTACTCTACCCAGCTGAGCTACGTATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                               <  1:560433/143‑1 (MQ=255)                                                                                                gTTCGTGGCGGGGTACCTTATCAAGTTGGGCTACGGATGCTTCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                               <  2:319244/149‑1 (MQ=18)                                                                                                    ctAGCCCCGTAATCTACCCCGCTGCGCTCCGGCTGCCTCTAGACATTTCTTTTCCGGCTGCTCTTGGATCCCCTACTAGATCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGt                                           <  2:1331525/148‑1 (MQ=12)                                                                                                      aGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCACCCGGGGGGATACAAACCTCCACCCGCTCGGTTCGGAGACGAGTAc                                         >  2:546039/1‑149 (MQ=17)                                                                                                        ccGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGGGGGTTCGACCCGCCGACCGGTCGGTTCGTAGCCGGGTActc                                       >  1:116269/1‑149 (MQ=18)                                                                                                               ctctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGACCCTCCGCCCGCTCGGTTCGTGGCCGAGTGCTCTATcccg                                >  1:115504/1‑147 (MQ=255)                                                                                                                tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGTTGCATCCGGGAGAGTTCGAACCTCCGACCGCGCGGTTCGTAGCGGGGTACTCTATCCAGc                               >  2:932321/1‑149 (MQ=21)                                                                                                                tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTCCGTCGCCGCGTACTCTAccccac                               >  2:413401/1‑146 (MQ=18)                                                                                                                tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCCAGCGCCCGGTTCGTCGCCGCGTCCTCTATCCacc                               >  2:19041/1‑147 (MQ=18)                                                                                                                      cAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTCAGAGATGGTGCATCCGGGAGGATTCGAAGCTCTGACCGCTAGGCTCGTCGCCGAGTACTCTATTCAGCTGTGCt                         >  1:661095/1‑149 (MQ=18)                                                                                                                       aGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTGTTCAGCTGGGcaa                        >  1:198585/1‑147 (MQ=255)                                                                                                                         cTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTCCTCTATCCAGCTGAGCTACg                      >  1:701012/1‑149 (MQ=255)                                                                                                                           gggCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                    <  1:121870/147‑1 (MQ=255)                                                                                                                                   gATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCGCGCTGAGCTACGGATGCAACgg            >  1:1122037/1‑149 (MQ=255)                                                                                                                                      gCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGCCCGCCCGGTTCGTAGCCGAGGCCTCTATCCAGCTGAGCTCCGGAGGCACCgggga         >  2:839169/1‑147 (MQ=25)                                                                                                                                       cATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaa        <  1:1401647/149‑1 (MQ=255)                                                                                                                                         tCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAact      >  1:396653/1‑149 (MQ=255)                                                                                                                                           gggAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAactta    >  2:909825/1‑149 (MQ=255)                                                                                                                                             gAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTAGACGCTTCGCGCGCCCGGCTCGTAGCCGTGTACTCTATCCGGTTGAGCTCCGGGTGCCTGGGGGTatttact  >  2:361481/1‑149 (MQ=17)                                                                                                                                                 |                                                                                                                                                 CCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACT  >  NC_000913/2817999‑2818287

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