Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	473,176	A→G	5.5%	T71A (ACG→GCG)	amtB →	ammonium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	473,176	0	A	G	5.5%	233.5 / 3.4	72	T71A (ACG→GCG)	amtB	ammonium transporter
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base G (2/2);  total (36/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 = 5.89e-01

GCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATGGGCAGCATTTATCAGT  >  NC_000913/473032‑473323                                                                                                                                                 |                                                                                                                                                    gCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGGTGCTGACGCAggtgacggt                                                                                                                                                 >  3:353409/1‑149 (MQ=255) gCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAggtgacggt                                                                                                                                                 <  4:48279/149‑1 (MQ=255) acACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAggtgacggt                                                                                                                                                 <  3:66764/148‑1 (MQ=255)    ccTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAggtgacggtgac                                                                                                                                              <  1:6846/149‑1 (MQ=255)     cTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACa                                                                                                                                             >  3:148244/1‑149 (MQ=255)     cTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACa                                                                                                                                             <  3:333037/149‑1 (MQ=255)        cGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACAttt                                                                                                                                          <  1:231885/149‑1 (MQ=255)                  tAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCt                                                                                                                                <  3:196869/149‑1 (MQ=255)                      gCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTAt                                                                                                                            >  2:70854/1‑149 (MQ=255)                         gACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGGCGGTGACATTTGCACTGGTCTGTATtct                                                                                                                         >  3:167474/1‑149 (MQ=255)                         gACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGGCGGTGACATTTGCACTGGTCTGTATtct                                                                                                                         >  1:328434/1‑149 (MQ=255)                         gACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCt                                                                                                                                <  4:678647/142‑1 (MQ=255)                         gACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCt                                                                                                                                >  3:678647/1‑142 (MQ=255)                         gACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATtct                                                                                                                         >  4:511000/1‑149 (MQ=255)                             aTGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTgg                                                                                                                     <  2:601079/149‑1 (MQ=255)                                cGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGggtg                                                                                                                  <  2:537133/149‑1 (MQ=255)                                 gTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTAGGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGggtgg                                                                                                                 >  4:518951/1‑149 (MQ=255)                                 gTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGggtgg                                                                                                                 >  2:206207/1‑149 (MQ=255)                                 gTTTATGATGATTTGTACTGCGCTGGTGCTGTCTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGTGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGggtgg                                                                                                                 <  2:593179/149‑1 (MQ=255)                                      tgatgaTTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTAc                                                                                                            <  3:511000/149‑1 (MQ=255)                                        atgatTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACgg                                                                                                          >  3:758256/1‑149 (MQ=255)                                         tgatTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGt                                                                                                         >  3:375970/1‑149 (MQ=255)                                            tttGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTAc                                                                                                      <  3:212575/149‑1 (MQ=255)                                                      gcTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATtct                                                                                                                         >  1:168982/1‑120 (MQ=255)                                                      gcTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATtct                                                                                                                         <  2:168982/120‑1 (MQ=255)                                                         ggTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTg                                                                                         <  1:70854/149‑1 (MQ=255)                                                           tGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGt                                                                                       <  4:353409/149‑1 (MQ=255)                                                           tGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGt                                                                                       >  1:368773/1‑149 (MQ=255)                                                           tGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGt                                                                                       >  2:311360/1‑149 (MQ=255)                                                              tGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAg                                                                                    >  1:330796/1‑149 (MQ=255)                                                              tGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAg                                                                                    <  3:662153/149‑1 (MQ=255)                                                              tGTTTATGACTATCCCGGGGATTGCCATGTTTTACGGTGGGTTGATTCGCGGCAAAAGCGTGCTGTCGAGGCTGAGGGAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCGGGCGTTTGGTGAg                                                                                    <  2:46372/149‑1 (MQ=255)                                                               gTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAgg                                                                                   <  2:529259/149‑1 (MQ=255)                                                                tttATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAg                                                                                    >  4:234607/1‑147 (MQ=255)                                                                tttATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAg                                                                                    <  3:234607/147‑1 (MQ=255)                                                                   aTGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTa                                                                                                                             >  4:575211/1‑103 (MQ=255)                                                                   aTGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTa                                                                                                                             <  3:575211/103‑1 (MQ=255)                                                                   aTGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGggt                                                                                                                   >  4:567614/1‑113 (MQ=255)                                                                   aTGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGggt                                                                                                                   <  3:567614/113‑1 (MQ=255)                                                                     gACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGcaaca                                                                             >  4:79220/1‑149 (MQ=255)                                                                     gACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGcaaca                                                                             >  3:487127/1‑149 (MQ=255)                                                                     gACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGcaaca                                                                             >  2:719082/1‑149 (MQ=255)                                                                        tATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAg                                                                                    <  4:599604/139‑1 (MQ=255)                                                                        tATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAg                                                                                    >  3:599604/1‑139 (MQ=255)                                                                        tATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAAct                                                                          >  3:501678/1‑149 (MQ=255)                                                                        tATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGAGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAAct                                                                          >  2:164118/1‑149 (MQ=255)                                                                              ggggATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCg                                                                    >  1:621886/1‑149 (MQ=255)                                                                                  aTTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACttcttc                                                                     >  3:321496/1‑144 (MQ=255)                                                                                  aTTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACttcttc                                                                     <  4:321496/144‑1 (MQ=255)                                                                                  aTTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCaa                                                                <  1:229916/149‑1 (MQ=255)                                                                                  aTTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGATGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCaa                                                                >  2:662825/1‑149 (MQ=255)                                                                                     gCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACAt                                                             >  1:216125/1‑149 (MQ=255)                                                                                       ccTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTa                                                           >  4:411734/1‑149 (MQ=255)                                                                                       ccTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTa                                                           >  2:694137/1‑149 (MQ=255)                                                                                            tttACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTgg                                                      <  1:719082/149‑1 (MQ=255)                                                                                             ttACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGt                                                     <  2:259056/149‑1 (MQ=255)                                                                                                  gTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATg                                                <  2:511407/149‑1 (MQ=255)                                                                                                   tGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGGCGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGc                                               <  4:167474/149‑1 (MQ=255)                                                                                                   tGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGGCGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGc                                               <  2:328434/149‑1 (MQ=255)                                                                                                      gTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGa                                            >  2:713888/1‑149 (MQ=255)                                                                                                           ttCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAAc                                       <  3:205012/149‑1 (MQ=255)                                                                                                            tCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACa                                      <  1:196527/149‑1 (MQ=255)                                                                                                              gcgGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATc                                    >  4:693024/1‑149 (MQ=255)                                                                                                                 gCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGaa                                 <  2:65174/149‑1 (MQ=255)                                                                                                                    aaaaCGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTg                              <  4:487127/149‑1 (MQ=255)                                                                                                                     aaaCGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGcaacaa                                                                            >  2:728437/1‑102 (MQ=255)                                                                                                                     aaaCGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGcaacaa                                                                            <  1:728437/102‑1 (MQ=255)                                                                                                                     aaaCGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGa                             <  1:90964/149‑1 (MQ=255)                                                                                                                           gCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGAcggcgg                       <  1:311360/149‑1 (MQ=255)                                                                                                                            cTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGt                      <  2:330796/149‑1 (MQ=255)                                                                                                                                 gATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGaa                                 <  1:217895/133‑1 (MQ=255)                                                                                                                                 gATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGaa                                 >  2:217895/1‑133 (MQ=255)                                                                                                                                 gATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATgg                 >  4:312258/1‑149 (MQ=255)                                                                                                                                        aCGCAGGTGACGGTGACATTTGCACTGGTCTGTCTTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATGGGCAGCAt          <  4:501678/149‑1 (MQ=255)                                                                                                                                            aGGTGACGGTGACATTTGCATTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATGGGCAGCATTTAt      >  4:409726/1‑149 (MQ=255)                                                                                                                                            aGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATGGGCAGCATTTAt      >  1:114365/1‑149 (MQ=255)                                                                                                                                              gtgacggtgacATTTGCACTGGTCTGTATTCTCTGGGTGGTTTGCGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATGGGCAGCATTTATCa    <  4:342759/149‑1 (MQ=255)                                                                                                                                                gacggtgacATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATGGGCAGCATTTATCAGt  >  3:272591/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGGTGCTGTTTATGACTATTCCGGGGATTGCCCTGTTTTACGGTGGGTTGATTCGCGGCAAAAACGTGCTGTCGATGCTGACGCAGGTGACGGTGACATTTGCACTGGTCTGTATTCTCTGGGTGGTTTACGGTTACTCGCTGGCGTTTGGTGAGGGCAACAACTTCTTCGGCAACATTAACTGGTTGATGCTGAAAAACATCGAACTGACGGCGGTGATGGGCAGCATTTATCAGT  >  NC_000913/473032‑473323

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