Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	473,813	C→A	14.7%	A283E (GCG→GAG)	amtB →	ammonium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	473,813	0	C	A	14.7%	188.5 / 17.7	67	A283E (GCG→GAG)	amtB	ammonium transporter
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (5/5);  total (33/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 = 8.95e-01

GTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCGGGTGGATGA  >  NC_000913/473667‑473960                                                                                                                                                   |                                                                                                                                                    gttTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCgt                                                                                                                                                   <  1:200158/149‑1 (MQ=255)     ggCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGtt                                                                                                                                                <  1:181250/148‑1 (MQ=255)       cTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTg                                                                                                                                             >  3:10273/1‑149 (MQ=255)        tttAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTgg                                                                                                                                            >  3:209369/1‑149 (MQ=255)        tttAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTgg                                                                                                                                            >  3:340360/1‑149 (MQ=255)              gCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGAt                                                                                                                                      <  3:296021/149‑1 (MQ=255)                 gggTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCCCTGCTGGGGGCGTGTTCTGGCGCGATTGc                                                                                                                                   >  4:113700/1‑149 (MQ=255)                 gggTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGc                                                                                                                                   >  4:380160/1‑149 (MQ=255)                   gTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCg                                                                                                                                 <  4:191635/149‑1 (MQ=255)                     cAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGt                                                                                                                               <  4:209369/149‑1 (MQ=255)                      aGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTc                                                                                                                              >  3:301758/1‑149 (MQ=255)                        cGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTg                                                                                                                            <  4:340360/149‑1 (MQ=255)                          ggCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGt                                                                                                                          <  8:151546/149‑1 (MQ=255)                           gCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTc                                                                                                                         >  7:7673/1‑149 (MQ=255)                               ggCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCg                                                                                                                     >  7:123056/1‑149 (MQ=255)                                gCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGt                                                                                                                    <  6:163973/149‑1 (MQ=255)                                  gAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGa                                                                                                                  >  8:77089/1‑149 (MQ=255)                                        aaTCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCa                                                                                                             <  4:10273/148‑1 (MQ=255)                                                cACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGAtt                                                                                                                                     >  1:134844/1‑116 (MQ=255)                                                cACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGAtt                                                                                                                                     <  2:134844/116‑1 (MQ=255)                                                     gCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACAt                                                                                               >  8:243207/1‑149 (MQ=255)                                                       aTTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTg                                                                                             <  4:154011/149‑1 (MQ=255)                                                        tttGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTgg                                                                                            >  6:168111/1‑149 (MQ=255)                                                        tttGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATATTCGGTGAATGGGCGCTGCGTGGTAAGCCATCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTTACTCCAGCCTGCGGCTCCATTgg                                                                                            >  5:12184/1‑149 (MQ=255)                                                          tgtgAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTggg                                                                                           <  6:56534/148‑1 (MQ=255)                                                           gtgAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGt                                                                                         >  3:356981/1‑149 (MQ=255)                                                              aaTACTTGGGGCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGGCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTgg                                                                                      <  3:203551/149‑1 (MQ=255)                                                              aaTACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCTCTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTgg                                                                                      >  1:41040/1‑149 (MQ=255)                                                              aaTACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTgg                                                                                      >  1:258017/1‑149 (MQ=255)                                                                tACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTggcg                                                                                    >  6:153337/1‑149 (MQ=255)                                                                  cTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGTGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTggcggc                                                                                  <  8:7673/149‑1 (MQ=255)                                                                   tgtgGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGgcggcg                                                                                 >  8:174851/1‑149 (MQ=255)                                                                     tgGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCGGCGGGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGgcggcg                                                                                 <  3:54229/147‑1 (MQ=255)                                                                       gTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCATGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGcgcg                                                                               >  6:57425/1‑147 (MQ=255)                                                                       gTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCATGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGcgcg                                                                               <  5:57425/147‑1 (MQ=255)                                                                          gCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGGCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGa                                                                           <  2:109445/148‑1 (MQ=255)                                                                          gCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCg                                                                                                      >  6:42647/1‑121 (MQ=255)                                                                             aCGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTAt                                                                       <  8:123056/149‑1 (MQ=255)                                                                               ggcggcggcAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCg                                                                     <  4:301758/149‑1 (MQ=255)                                                                                        aaTTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGTAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAg                                                            >  5:129256/1‑149 (MQ=255)                                                                                             ttGGCTGGAACTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGt                                                       <  7:206359/149‑1 (MQ=255)                                                                                                gCTGGATCTTCGGTGCATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGGGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCt                                                     <  7:77089/148‑1 (MQ=255)                                                                                                    gATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGt                                                                                                                          >  8:121612/1‑75 (MQ=255)                                                                                                    gATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGt                                                                                                                          <  7:121612/75‑1 (MQ=255)                                                                                                        ttCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCtt                                            >  2:218285/1‑149 (MQ=255)                                                                                                         tCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTtg                                           >  4:340938/1‑149 (MQ=255)                                                                                                           ggTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTtgtg                                         <  5:54038/149‑1 (MQ=255)                                                                                                              gAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTgggg                                      >  1:222736/1‑149 (MQ=255)                                                                                                               actGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTATGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGc                                     <  3:294393/147‑1 (MQ=255)                                                                                                                    gcgcTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTAc                                <  5:153337/149‑1 (MQ=255)                                                                                                                        tGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATg                            >  1:195337/1‑149 (MQ=255)                                                                                                                        tGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTATGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATg                            >  7:170540/1‑149 (MQ=255)                                                                                                                              gTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTgggg                                                                                          <  3:172106/81‑1 (MQ=255)                                                                                                                              gTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTgggg                                                                                          >  4:172106/1‑81 (MQ=255)                                                                                                                              gTAAGCCTTCACTGCTGGGGGCGTGTGCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCaaa                      <  1:372190/149‑1 (MQ=255)                                                                                                                               tAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAAc                     <  2:41040/149‑1 (MQ=255)                                                                                                                               tAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAAc                     >  8:278886/1‑149 (MQ=255)                                                                                                                                 aGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGc                   >  2:282217/1‑149 (MQ=255)                                                                                                                                 aGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACg                    <  4:356981/148‑1 (MQ=255)                                                                                                                                   cattCACTGCGGGGGGAGTGTTCTGGCGCGATTGACGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCtt                 <  6:129256/147‑1 (MQ=255)                                                                                                                                     ttCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTtgc               >  3:126872/1‑149 (MQ=255)                                                                                                                                        aCTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTtgctgc            >  7:102726/1‑149 (MQ=255)                                                                                                                                            ctgGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCggg         <  1:282217/148‑1 (MQ=255)                                                                                                                                            ctgGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCggg         <  7:243207/148‑1 (MQ=255)                                                                                                                                            ctgGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGGGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCgg          <  3:380160/147‑1 (MQ=255)                                                                                                                                               ggggCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCGGGTGga     >  2:399358/1‑149 (MQ=255)                                                                                                                                                  gCGTGTGATGGCTGGATTGCCTGTAGGGTCGGCGGGTCGCCAGCCTGTGGTTAGATTGGGGGTGGAGGCGCGTTTATTATAGGCGTGGTATCTGGGCTGGCGGGCTTGTGGGGCGGTACCATGCTGAAACGCTTGCTGCGGGTGgaaga  <  6:12184/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCGGGTGGATGA  >  NC_000913/473667‑473960

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