Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,895,485	G→C	6.5%	W227C (TGG→TGC)	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,485	0	G	C	6.5%	193.7 / 2.4	60	W227C (TGG→TGC)	pabB	aminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base C (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 = 1.00e+00

CGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAA  >  NC_000913/1895340‑1895626                                                                                                                                                  |                                                                                                                                              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                                                                                   >  4:93816/1‑133 (MQ=255)                                                                          tCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCtt                                                                                   <  3:93816/133‑1 (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                                         <  1:254241/149‑1 (MQ=255)                                                                                                    tGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTAtt                                         <  6:55885/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                                                                                                      >  6:331038/1‑58 (MQ=255)                                                                                                                                  aTTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCAtt                                                                                                      <  5:331038/58‑1 (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:148643/147‑1 (MQ=255)                                                                                                                                           aggAATGCCTGCCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTaaa  <  2:153329/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)                                                                                                                                                  |                                                                                                                                              CGCTCACTTCCGACTGGCAATCCAATATGACCCGCGAGCAGTACGGCGAAAAATTTCGCCAGGTACAGGAATATCTGCACAGCGGTGATTGCTATCAGGTGAATCTCGCCCAACGTTTTCATGCGACCTATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAA  >  NC_000913/1895340‑1895626

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: ‑