Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,215	(G)8→9	noncoding (72/87 nt)	leuP ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,207	1	.	G	95.5%	208.2 / 3.0	66	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (34/29);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606074‑4606341                                                                                                                                         |                                                                                                                                         aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑gggggg                                                                                                                                   <  1:17046/140‑1 (MQ=255)          gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGa                                                                                                                          <  1:883568/140‑1 (MQ=255)          gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGa                                                                                                                          >  1:1245318/1‑140 (MQ=255)          gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑ggggggggggacttg                                                                                                                          >  2:1272711/1‑133 (MQ=255)          gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑ggggggggggacttg                                                                                                                          >  2:403499/1‑133 (MQ=255)           cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGaa                                                                                                                         >  1:721351/1‑140 (MQ=255)           cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑GGGGGGGGA‑CTTGAAc                                                                                                                        >  1:1369130/1‑140 (MQ=255)                ggggggACTTGAACCCCCGCGTCCGTAAAGACACTAACCCCTTAAGGTTGCGCGTCTACCAAATCCGCCACCTTCGCACCGTCAACTTAATTTTTTTGAAATTGCCCCGTTTTGTGCG‑GA‑‑GGGGGGGGA‑CTTGAAccccc                                                                                                                    <  2:1052570/140‑1 (MQ=18)                ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAAccccc                                                                                                                    <  1:198056/140‑1 (MQ=255)                  ggggACTTGAACCCCCACGTCCGTAAGGACACTAACTCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCAc                                                                                                                  <  1:142129/140‑1 (MQ=255)                   gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTAATGATGTCTCTTCGTCTGGTGC‑‑GA‑‑GGCTGCGGACCATGAGCCCCCACg                                                                                                                 >  2:501877/1‑140 (MQ=17)                    ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGGGGGGGGGA‑CTTGAACCCCCACg                                                                                                                 <  1:403499/140‑1 (MQ=255)                            aCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGGGGGGGGGA‑CTTGAACCCCCACGTCCGTAAg                                                                                                         <  1:1272711/140‑1 (MQ=255)                               cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAaca                                                                                                     <  2:721351/140‑1 (MQ=255)                                ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAacac                                                                                                    >  1:108488/1‑140 (MQ=255)                                 cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACt                                                                                                   <  1:637641/140‑1 (MQ=255)                                 cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACt                                                                                                   >  2:602231/1‑140 (MQ=255)                                  aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGATTGGTGC‑‑GAGGGGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACt                                                                                                   >  2:208836/1‑140 (MQ=255)                                   cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                 <  1:917317/140‑1 (MQ=255)                                         tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                                                           <  1:293241/140‑1 (MQ=255)                                              acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                                                      >  2:904955/1‑140 (MQ=255)                                               cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                     <  1:55507/140‑1 (MQ=255)                                               cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                     <  1:602231/140‑1 (MQ=255)                                                acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                    >  1:156433/1‑140 (MQ=255)                                                acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                    >  1:1269557/1‑140 (MQ=255)                                                  tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                  >  2:539705/1‑140 (MQ=255)                                                    acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                                                                >  1:276171/1‑140 (MQ=255)                                                       ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                             >  1:295301/1‑140 (MQ=255)                                                       ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                             >  2:696411/1‑140 (MQ=255)                                                          gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                          >  1:1140526/1‑140 (MQ=255)                                                          gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                          <  1:539705/140‑1 (MQ=255)                                                                aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGC‑CTTGAACCCCCACGTCCGTAAAAACCCTAACACCTGAAGCTCGCACGTCTACCAATTcc                                                                    >  2:376588/1‑140 (MQ=37)                                                                aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                        <  1:1251882/136‑1 (MQ=255)                                                                aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                        >  2:1251882/1‑136 (MQ=255)                                                                aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                    >  2:394719/1‑140 (MQ=255)                                                                 gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:1369130/140‑1 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  >  1:583999/1‑140 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  1:472021/140‑1 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:1140526/140‑1 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  1:374828/140‑1 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  1:228426/140‑1 (MQ=255)                                                                  cgcgTCTACCAAATCCGGCACCTTAGCACAGTCAACTTACTTTTTTTGATATCGCCTCGTTTTGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:414539/140‑1 (MQ=21)                                                                  cgcgTCTACCAAAACCGCCGCCTTCGAACAGTCAACTTACTTTTTATGATATCAACTCGATTTGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:256696/140‑1 (MQ=18)                                                                            aaTTCCGCCGCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  1:722055/140‑1 (MQ=255)                                                                            aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  2:108488/140‑1 (MQ=255)                                                                            aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  2:295301/140‑1 (MQ=255)                                                                            aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  1:696411/140‑1 (MQ=255)                                                                                ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                    <  1:895908/140‑1 (MQ=255)                                                                                   ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                 <  1:904955/140‑1 (MQ=255)                                                                                    cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGGGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                 <  1:208836/140‑1 (MQ=255)                                                                                         tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        >  2:1251833/1‑127 (MQ=255)                                                                                         tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  1:1251833/127‑1 (MQ=255)                                                                                           gCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat                                         >  1:1140171/1‑140 (MQ=255)                                                                                               aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                     <  2:276171/140‑1 (MQ=255)                                                                                                 tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGGGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATAAACg                                    >  2:603097/1‑140 (MQ=255)                                                                                                 tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                   >  2:209033/1‑140 (MQ=255)                                                                                                 tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                   >  2:590504/1‑140 (MQ=255)                                                                                                   aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                                 >  2:1150068/1‑140 (MQ=255)                                                                                                    tCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                >  1:857249/1‑140 (MQ=255)                                                                                                         cTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                           >  2:268426/1‑140 (MQ=255)                                                                                                           ttttttGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                         >  1:1161181/1‑140 (MQ=255)                                                                                                            tttttGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                        <  1:209033/140‑1 (MQ=255)                                                                                                            tttttGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                        <  2:583999/140‑1 (MQ=255)                                                                                                              tttGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGcgcg                      >  1:838848/1‑140 (MQ=255)                                                                                                              tttGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                      <  1:178677/140‑1 (MQ=255)                                                                                                               ttGATATCGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                     >  2:801404/1‑140 (MQ=255)                                                                                                                     tCGCCTTGTTTGGTGCGGGG‑‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggg                <  1:488958/140‑1 (MQ=21)                                                                                                                      cGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg              >  1:497664/1‑140 (MQ=255)                                                                                                                      cGCCTCGTTTGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg              >  2:171572/1‑140 (MQ=255)                                                                                                                             tttGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGaa      >  1:814760/1‑140 (MQ=255)                                                                                                                             tttGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGaa      >  2:1138968/1‑140 (MQ=255)                                                                                                                             tttGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGaa      >  1:1144316/1‑140 (MQ=255)                                                                                                                               tGGTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAAcc    >  2:1332360/1‑140 (MQ=255)                                                                                                                                 gTGC‑‑GAG‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAAcccc  >  2:57601/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                         ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑GGGGGGGGA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606074‑4606341

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