Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,895,487	A→T	6.7%	Q228L (CAG→CTG)	pabB →	aminodeoxychorismate synthase, subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,895,487	0	A	T	6.7%	188.8 / 3.6	60	Q228L (CAG→CTG)	pabB	aminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base T (2/2);  total (26/34)
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.18e-01

CGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGC  >  NC_000913/1895340‑1895633                                                                                                                                                    |                                                                                                                                                   cGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcag                                                                                                                                                   <  8:44956/149‑1 (MQ=255) cGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcag                                                                                                                                                   <  7:287961/149‑1 (MQ=255)  gCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcagg                                                                                                                                                  >  2:170065/1‑149 (MQ=255)   cTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcagg                                                                                                                                                  >  8:28970/1‑148 (MQ=255)    tCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATggcaggca                                                                                                                                                >  2:53202/1‑149 (MQ=255)       cTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTc                                                                                                                                             <  8:44342/149‑1 (MQ=255)       cTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTc                                                                                                                                             >  2:309716/1‑149 (MQ=255)       cTTCCGACTGGCAATACAATATGACCCGGGAGCAGTACGGTGAAAATTTTCGCCAGGTACAGGAATATCTGCAGAGCGGTTATTTCTATCAGGATTATCTCGCCCATCGTTTTCCTTAGACCTAATCTGCCGATGAATGGCAGGCATTc                                                                                                                                             <  6:122403/149‑1 (MQ=255)           cGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCtt                                                                                                                                          <  2:321200/148‑1 (MQ=255)            gACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCa                                                                                                                                        <  5:98608/149‑1 (MQ=255)               tGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCt                                                                                                                                     >  1:277951/1‑149 (MQ=255)                 gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                    >  5:38339/1‑148 (MQ=255)                 gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                    >  4:158164/1‑148 (MQ=255)                 gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                    >  2:231723/1‑148 (MQ=255)                 gCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCtt                                                                                                                                    >  6:82898/1‑148 (MQ=255)                   aaTCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAAt                                                                                                                                 >  3:283227/1‑149 (MQ=255)                    aTCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAAt                                                                                                                                 >  4:80322/1‑148 (MQ=255)                     tCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGGGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCa                                                                                                                               <  5:241178/149‑1 (MQ=255)                       cAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAg                                                                                                                              <  8:12563/148‑1 (MQ=255)                         atatGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGc                                                                                                                            >  6:69333/1‑148 (MQ=255)                            tGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCaa                                                                                                                         <  7:173850/148‑1 (MQ=255)                             gACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAAcc                                                                                                                       >  5:55885/1‑149 (MQ=255)                                 cgcgAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACcgcgc                                                                                                                   <  1:53202/149‑1 (MQ=255)                                            cGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgc                                                                                                         <  1:309716/148‑1 (MQ=255)                                             ggCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAgcgc                                                                                                       <  6:254464/149‑1 (MQ=255)                                                     aTTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTAc                                                                                               >  7:281805/1‑149 (MQ=255)                                                         cGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTc                                                                                            <  1:55093/148‑1 (MQ=255)                                                               gTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGaa                                                                                       >  1:11923/1‑147 (MQ=255)                                                               gTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACa                                                                                     <  7:28970/149‑1 (MQ=255)                                                                       atatCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCaa                                                                             <  3:80322/149‑1 (MQ=255)                                                                          tCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCtt                                                                                          <  3:93816/133‑1 (MQ=255)                                                                          tCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCtt                                                                                          >  4:93816/1‑133 (MQ=255)                                                                          tCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAAttt                                                                          >  2:136848/1‑149 (MQ=255)                                                                              cacaGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCTACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTaa                                                                       <  4:101430/148‑1 (MQ=255)                                                                              cacaGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTaa                                                                       <  3:158164/148‑1 (MQ=255)                                                                              cacaGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTa                                                                        <  5:82898/147‑1 (MQ=255)                                                                                        aTTTCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATGCTGGCGATTAATGTCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGcc                                                             <  7:228807/148‑1 (MQ=255)                                                                                            cTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCagag                                                         <  8:281805/148‑1 (MQ=255)                                                                                               tCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTCCTCAGGCCCACCGCGCGCCATTTAGCGCTTGTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCgg                                                      <  3:49967/148‑1 (MQ=255)                                                                                                 aGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGtt                                                    >  3:238313/1‑148 (MQ=255)                                                                                                  ggTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTa                                                  <  1:318008/149‑1 (MQ=255)                                                                                                    tGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTAtt                                                <  6:55885/149‑1 (MQ=255)                                                                                                    tGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTAtt                                                <  1:254241/149‑1 (MQ=255)                                                                                                           cGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTGCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTtgt                                          >  6:56083/1‑148 (MQ=255)                                                                                                              ccAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGat                                       >  7:146110/1‑148 (MQ=255)                                                                                                                   gTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGt                                 <  6:56353/149‑1 (MQ=255)                                                                                                                     tttCATGAGACCTATTCTGGCGATGAATGGCAGTCATTCCTTCAGATTAATCAGGCCAACCGCGCGCCATTTAGAGCGTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCGAGAGCGGTTTATGCTTGGTGATAATAGTg                                >  6:254737/1‑148 (MQ=255)                                                                                                                            cGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCa                         <  4:283227/148‑1 (MQ=255)                                                                                                                            cGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCa                         <  8:17664/148‑1 (MQ=255)                                                                                                                             gACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCTGTTTATTCTTTGTGATAATAGTGAAATCCAGa                       <  5:159205/149‑1 (MQ=255)                                                                                                                               ccTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGAcc                     >  3:144975/1‑149 (MQ=255)                                                                                                                               ccTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAGTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCATAcc                     >  6:318338/1‑149 (MQ=255)                                                                                                                                  aTTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAtt                                                                                                             <  5:331038/58‑1 (MQ=255)                                                                                                                                  aTTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAtt                                                                                                             >  6:331038/1‑58 (MQ=255)                                                                                                                                      tGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGa              >  7:191172/1‑149 (MQ=255)                                                                                                                                      tGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGa              >  6:35539/1‑149 (MQ=255)                                                                                                                                        tgaaggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAt             >  1:153329/6‑148 (MQ=255)                                                                                                                                        tgaaggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAt             >  1:148643/6‑148 (MQ=255)                                                                                                                                        gCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGAtt            <  6:38339/149‑1 (MQ=255)                                                                                                                                          gATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaa          <  1:136848/149‑1 (MQ=255)                                                                                                                                           aggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa         <  2:153329/147‑1 (MQ=255)                                                                                                                                           aggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa         <  2:148643/147‑1 (MQ=255)                                                                                                                                           aTGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa         <  2:277951/149‑1 (MQ=255)                                                                                                                                                tGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGa                       >  7:328036/1‑130 (MQ=255)                                                                                                                                                tGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGa                       <  8:328036/130‑1 (MQ=255)                                                                                                                                                   caggcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCAGATTAAAGGCACGc  >  6:111904/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGC  >  NC_000913/1895340‑1895633

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