Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	187,510	A→G	5.3%	Y381H (TAC→CAC)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	187,510	0	A	G	5.3%	230.9 / 2.8	75	Y381H (TAC→CAC)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base G (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

GGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGC  >  NC_000913/187367‑187658                                                                                                                                                |                                                                                                                                                     ggggTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACa                                                                                                                                                 <  4:274774/149‑1 (MQ=255)  ggttgccGCAGAAGGATCAACAACAGTTTTCGTGCTTCCGGAATCTGACACAGCGGGTGGTGCAGCTGGCGACGGGCATGGCGTAACTGGCGCCGCGGTGTGGAGTAAATGCCGGTGCTCGCACTGTTGTGCACCATGGTGTAGAACAt                                                                                                                                                <  6:45211/143‑1 (MQ=255)       aCGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGc                                                                                                                                            >  1:230013/1‑148 (MQ=255)         gCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACCGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGc                                                                                                                                            >  3:243020/1‑146 (MQ=255)           aGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGAt                                                                                                                                       <  7:320407/149‑1 (MQ=255)                     aaaaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGc                                                                                                                             >  4:13712/1‑149 (MQ=255)                       aaCAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTGGAACATACGCAAGATGGCTTCCGGCTg                                                                                                                           >  6:313940/1‑149 (MQ=255)                         cAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGgc                                                                                                                         <  7:203358/149‑1 (MQ=255)                          aGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGgcg                                                                                                                        >  4:117702/1‑149 (MQ=255)                           gTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGgcg                                                                                                                        >  7:67321/1‑148 (MQ=255)                              ttCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGGGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGgcgctt                                                                                                                     >  2:78271/1‑146 (MQ=255)                                  cgatctaggAATATTTAAAAGAGGTTGTTGCAGATTGAGACGGGAATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATa                                                                                                                <  6:219136/142‑1 (MQ=255)                                       ccGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTGGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAgtgt                                                                                                            <  5:313940/148‑1 (MQ=255)                                        cGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGtt                                                                                                           <  8:288272/148‑1 (MQ=255)                                        cGGAATATTACAAAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGtt                                                                                                           <  7:75544/148‑1 (MQ=255)                                         ggAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAACTGCCGGTGATCGCACTGTTGTGCCCCATGGTGTAGAACATACTCAAGATTGCTTCCGGCTGTCGCATAAATatagatct                                                                                                         >  5:219136/1‑142 (MQ=255)                                         ggAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATa                                                                                                                <  8:323242/142‑1 (MQ=255)                                         ggAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATa                                                                                                                >  7:323242/1‑142 (MQ=255)                                          gAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTtca                                                                                                        <  5:54390/149‑1 (MQ=255)                                                   acaGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAgg                                                                                               >  8:2533/1‑149 (MQ=255)                                                   acaGCGGTAGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGAGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCAAGGTGTAGAACATACTCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAgg                                                                                               >  3:224603/1‑149 (MQ=255)                                                    caGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGt                                                                                              <  7:75724/149‑1 (MQ=255)                                                     aGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTc                                                                                             >  7:337295/1‑149 (MQ=255)                                                     aGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTc                                                                                             <  7:188061/149‑1 (MQ=255)                                                      gCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTc                                                                                             >  7:103160/1‑148 (MQ=255)                                                        ggTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGAt                                                                                          <  1:50066/149‑1 (MQ=255)                                                        ggTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGAt                                                                                          <  2:273373/149‑1 (MQ=255)                                                            gttgCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAg                                                                                       <  6:153939/148‑1 (MQ=255)                                                            gttgCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATCGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAg                                                                                       <  6:280058/148‑1 (MQ=255)                                                               gCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTtca                                                                                                        >  3:231411/1‑128 (MQ=255)                                                               gCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTtca                                                                                                        <  4:231411/128‑1 (MQ=255)                                                                 aGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTGGAACATACGCAAGATGGCTTc                                                                                                                                 >  1:259582/1‑101 (MQ=255)                                                                 aGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTGGAACATACGCAAGATGGCTTc                                                                                                                                 <  2:259582/101‑1 (MQ=255)                                                                 aGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTtcatc                                                                                                      >  1:198276/1‑128 (MQ=255)                                                                 aGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTtcatc                                                                                                      <  2:198276/128‑1 (MQ=255)                                                                 aGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTAcc                                                                                 >  6:26142/1‑149 (MQ=255)                                                                         aCGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCt                                                                         <  2:78999/149‑1 (MQ=255)                                                                             gCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGaaa                                                                     >  3:114387/1‑149 (MQ=255)                                                                               aTGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACt                                                                   >  6:325838/1‑149 (MQ=255)                                                                                       aCTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCAtcgtc                                                            <  5:131148/148‑1 (MQ=255)                                                                                           gcgcAGCGTGGTGGAGTTAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGAGCATAAATAGTTTTTCATCTCGCATGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGAtt                                                        <  4:224603/148‑1 (MQ=255)                                                                                           gcgcAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTg                                                       >  3:121714/1‑149 (MQ=255)                                                                                            cgcAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTGATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTg                                                       >  6:258851/1‑148 (MQ=255)                                                                                                   tggtggAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGt                                               >  3:86253/1‑149 (MQ=255)                                                                                                       ggAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGttttt                                           >  2:45714/1‑149 (MQ=255)                                                                                                        gAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTtc                                          >  4:99964/1‑149 (MQ=255)                                                                                                        gAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTtc                                          <  1:153124/149‑1 (MQ=255)                                                                                                        gAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTtc                                          <  3:13712/149‑1 (MQ=255)                                                                                                          gTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGaaa                                                                     >  5:322585/1‑120 (MQ=255)                                                                                                          gTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGaaa                                                                     <  6:322585/120‑1 (MQ=255)                                                                                                           tAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTtcgtc                                       >  2:41714/1‑149 (MQ=255)                                                                                                             aaTGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTtcgtcg                                      <  8:337295/148‑1 (MQ=255)                                                                                                              aTGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCGAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTtc                                          >  8:49363/1‑143 (MQ=255)                                                                                                              aTGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCGAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTtc                                          <  7:49363/143‑1 (MQ=255)                                                                                                               tGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTcggcg                                   <  3:117702/149‑1 (MQ=255)                                                                                                                 ccGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAATATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTcggcg                                   <  1:312636/147‑1 (MQ=255)                                                                                                                 ccGGTGATCGCACTGTTGTCCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGAACCGAGTATCTGAAACTCATCGTCGATTGGACGAGGTTTTTCGTcggcg                                   >  2:312636/1‑147 (MQ=255)                                                                                                                   ggTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGaa                               <  2:230013/149‑1 (MQ=255)                                                                                                                    gTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGt                                                                                              <  6:143634/85‑1 (MQ=255)                                                                                                                    gTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGt                                                                                              >  5:143634/1‑85 (MQ=255)                                                                                                                        tCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGc                           >  2:259927/1‑148 (MQ=255)                                                                                                                          gCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGa                         <  7:2533/148‑1 (MQ=255)                                                                                                                          gCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAg                        >  6:318051/1‑149 (MQ=255)                                                                                                                           cACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGGCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAgg                       <  5:26142/149‑1 (MQ=255)                                                                                                                             cTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGAt                     <  4:114387/149‑1 (MQ=255)                                                                                                                              tgttgtGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAACCTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATg                    >  5:17864/1‑149 (MQ=255)                                                                                                                                ttgtGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGc                  >  2:201299/1‑149 (MQ=255)                                                                                                                                 tgtgCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCt                 >  8:284717/1‑149 (MQ=255)                                                                                                                                   tgCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCtt                >  8:309212/1‑148 (MQ=255)                                                                                                                                     cACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCAt             <  5:325838/149‑1 (MQ=255)                                                                                                                                      aCCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATc            <  4:35559/149‑1 (MQ=255)                                                                                                                                       ccATGGTGTAGAACATACGCAAGATGGCGTCCGGCTGGCGCATAAATAGTGTTTACTCACGCTTGTCGATTAGCGTACCGAGTAGCTGAAACTCATCGTCGATTGTACGTGGTTTTGAGTCGGCGGGAAGGGCGAGTATGGCTTCATc            >  4:258758/1‑148 (MQ=255)                                                                                                                                           ggTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGaa         >  7:58011/1‑147 (MQ=255)                                                                                                                                              gtAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGTGGGAAGGGCGAGGATGGCTTCATCGAACAGtt    >  7:75406/1‑149 (MQ=255)                                                                                                                                              gtAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGACGGGCGAGGATGGCTTCAGCGACCAGtt    >  4:160948/1‑149 (MQ=255)                                                                                                                                                aGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTgc  >  7:329837/1‑149 (MQ=255)                                                                                                                                                aGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTG‑AACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTg   >  2:302517/1‑147 (MQ=255)                                                                                                                                                |                                                                                                                                                     GGGGTGACGCAGAATGCTCAAAAACAGTTTTCGTGCTTCCGGAATATTACACAGCGGTTGTTGCAGATGGCGACGGGCATGGCGTAACTGGCGCAGCGTGGTGGAGTAAATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGC  >  NC_000913/187367‑187658

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