Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,511,347	A→G	9.7%	V25A (GTC→GCC)	yhfY ←	PRD domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,511,347	0	A	G	9.7%	172.0 / 8.2	69	V25A (GTC→GCC)	yhfY	PRD domain protein
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base G (4/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

GCCAGTTCTGCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGT  >  NC_000913/3511200‑3511495                                                                                                                                                     |                                                                                                                                                     gcCAGTTCTGCCAGCAACTCGTTATCCAGCGGCTCTCTTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTg‑acgac                                                                                                                                                     <  2:848808/149‑1 (MQ=255)    aGTTCTGCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTg                                                                                                                                                  <  1:797850/149‑1 (MQ=255)        cTGCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGcatc                                                                                                                                              <  1:719750/149‑1 (MQ=255)         tGCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTT‑ACGACCTGcatca                                                                                                                                             <  1:622607/149‑1 (MQ=255)          gCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGGCCTGcatcat                                                                                                                                            >  2:159627/1‑149 (MQ=255)          gCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGGCCTGcatcat                                                                                                                                            >  2:717504/1‑149 (MQ=255)          gCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGGCCTGcatcat                                                                                                                                            >  2:348311/1‑149 (MQ=255)           ccAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATg                                                                                                                                           <  1:872683/149‑1 (MQ=255)             cgCAACTCGTTATCCAGCGGCTCTATTTCTTCACCCCGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGcc                                                                                                                                         <  2:645959/148‑1 (MQ=255)             aGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGTCCTGCATCATGcc                                                                                                                                         >  1:166839/1‑149 (MQ=255)                                 ctctATTTCTTCACCGCGGCGCCTGCGCATCAGTGCCCTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCaa                                                                                                                     <  1:47947/149‑1 (MQ=255)                                  tctATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAAt                                                                                                                    >  2:813655/1‑149 (MQ=255)                                  tctAATTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCAGCAGGCCTTTGCAGATGTCCTTATCAAt                                                                                                                    >  2:417380/1‑149 (MQ=255)                                   ctATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCGTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATa                                                                                                                   >  1:608505/1‑149 (MQ=255)                                    tATTTCTTCCCCACGGCGACTGCGCATCAGTGCAGTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATaa                                                                                                                  <  2:56524/149‑1 (MQ=255)                                    tATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATaa                                                                                                                  <  1:263361/149‑1 (MQ=255)                                    tATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTG‑TCCGATAACGTGGTACGGCCTCAAACATGCCTTGGCAGATTTCCTTATCAATac                                                                                                                  >  1:887284/1‑148 (MQ=255)                                    tATTTCTTCACCACGGCGACTGCGCACCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATaa                                                                                                                  >  1:844420/1‑149 (MQ=255)                                    tATTTCTGCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATaa                                                                                                                  <  1:714191/149‑1 (MQ=255)                                      tttCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTCCAGATGTCCTTATCAATAACg                                                                                                                >  2:815763/1‑149 (MQ=255)                                      tttCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGGACGTTG‑ACGACCTGCATCATGCATTTGCAGATGTCCTTATCAATAACg                                                                                                                >  1:400261/1‑149 (MQ=255)                                        ccctcACCACGGCGACTGCGCATCAGTGCCCTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCCGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACgcc                                                                                                              <  1:83488/146‑1 (MQ=255)                                           tcACCCCGGCGACGGCGCATCAGGGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACgcctgc                                                                                                           <  2:351948/149‑1 (MQ=255)                                             accacGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACgcctgcct                                                                                                         >  1:699985/1‑149 (MQ=255)                                                acGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGc                                                                                                      <  2:438025/149‑1 (MQ=255)                                                acGGCGACTGCGCATCAGTGCACTCGACATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGc                                                                                                      >  2:898493/1‑149 (MQ=255)                                                  ggCGACTGCGCATCAGTGAACTCGCCATACGTGCCCTCCACATCGCTCCTTGCTCACTGGGCACCGGCAGATGGCACTCTGTTTCCCGTACATTG‑ACGACCTGCCTGCTGCCTCTCCACATGTCCTTAACGATGACGCCTGCCTCGCaa                                                                                                    >  2:383171/1‑149 (MQ=255)                                                   gCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCaaa                                                                                                   >  1:812077/1‑149 (MQ=255)                                                      aCTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCa                                                                                                <  1:553212/149‑1 (MQ=255)                                                        tGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAgg                                                                                              >  2:399851/1‑149 (MQ=255)                                                          cgcATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGtt                                                                                            >  2:613566/1‑149 (MQ=255)                                                            cagcAGTGCACTCGCCATAGGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACTCCTGCCTCGCAAAGCAGGTTGa                                                                                          <  2:771337/146‑1 (MQ=255)                                                             aTCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATCACGCCTGCCTCGCAAAGCAGGTTGAg                                                                                         >  1:187210/1‑149 (MQ=255)                                                                    cACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCCGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGtt                                                                                  >  2:802205/1‑149 (MQ=255)                                                                        cGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCa                                                                              <  1:548028/149‑1 (MQ=255)                                                                          ccATATGTGTCATCGCCATCGTTCCTTGCTCACGCCGCACCGGCAGCTGGCACTCTGGTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGAAGTCCTTATCAATAACGCCCGCCTCGCAAAGCAGGTTGAGTCTGGCTTCCata                                                                            >  2:693762/1‑149 (MQ=255)                                                                          ccAAATGTGTCATCGCCATCGTACCTTGCACACTGCGCACCGGCAGATGGCACTCGGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAACGCAGGTTGAGTCCGGTTTCCaca                                                                            >  2:312353/1‑147 (MQ=255)                                                                            atatGTGTCTGCGCCTTCGTTCCCTGCTCTCTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTTCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATc                                                                          <  2:770447/149‑1 (MQ=255)                                                                              atGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGGCCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAg                                                                        <  1:159627/149‑1 (MQ=255)                                                                              atGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGGCCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAg                                                                        <  1:348311/149‑1 (MQ=255)                                                                              atGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGGCCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAg                                                                        <  1:717504/149‑1 (MQ=255)                                                                                 tgtCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGTTTGAGTCTGGTTTACATATCAGGTg                                                                     >  1:75663/1‑149 (MQ=255)                                                                                     aTCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCtt                                                                 >  1:257737/1‑149 (MQ=255)                                                                                     aTCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTAAGTTA‑CCCACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCtt                                                                 <  1:891804/149‑1 (MQ=255)                                                                                        gCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCg                                                              <  2:47421/149‑1 (MQ=255)                                                                                         ccATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGt                                                             >  1:775675/1‑149 (MQ=255)                                                                                            tCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTta                                                          <  1:385671/149‑1 (MQ=255)                                                                                            tCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTta                                                          >  2:176174/1‑149 (MQ=255)                                                                                              gTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTCTCAATAACGCCTGCGTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTtaat                                                        >  1:890172/1‑149 (MQ=255)                                                                                              gTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGGACGTTG‑GCGACCTGCATGATGTCTTTGCAGATATCCTTAACAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCGTATCAGCTGCCTTTCGttgcat                                                        >  2:255458/1‑145 (MQ=255)                                                                                              gTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGCTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTCGTTTCCATATCAGATGCCTTTCGTTtaat                                                        >  2:223125/1‑149 (MQ=255)                                                                                              gTTCCTTGCTCACTGCGAACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGTATCATGCCTTTGCAGATGTCCTTATCAATAACGCCAGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTtaat                                                        >  2:673582/1‑149 (MQ=255)                                                                                                    tGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCCGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCGGGTGGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCg                                                  <  2:831048/149‑1 (MQ=255)                                                                                                    tGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCCTATCAGATGCCTTTGGTTTAATAACtag                                                  >  2:894102/1‑147 (MQ=255)                                                                                                    tGCCCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCg                                                  >  2:340930/1‑149 (MQ=255)                                                                                                          cTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCt                                            <  1:712701/149‑1 (MQ=255)                                                                                                          cTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCt                                            <  1:193317/149‑1 (MQ=255)                                                                                                                 ccGGCAGATGGCACTCTGTTTGCAGTACGTGG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTcccc                                     >  1:476937/1‑149 (MQ=255)                                                                                                                  cGGCAGATGGCACTCTGTTTCCAGTTCGGTG‑ACGACCTCCATCACTCCTTTGCAGATGCCCCTATCAAAAATGCCTGGCTCGCAAATCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCca                                    <  2:328172/149‑1 (MQ=255)                                                                                                                  cGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCca                                    >  1:364299/1‑149 (MQ=255)                                                                                                                  cGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCca                                    >  2:768356/1‑149 (MQ=255)                                                                                                                      aGATGCCACTCCGCTTCCCGGACTTTG‑ACCACCTGCATCATGCCGTCGCAGATGTCCTTATCAAAAACGCCTGCCTCGAAAAGCACGTTGACTCTGGCTTCCATATCAGATGCCTTTCCTTTAATAACTCGGGTCCTGTTCCCcactca                                <  2:281979/149‑1 (MQ=255)                                                                                                                        aTGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGGCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCcactcact                              >  1:50638/1‑149 (MQ=255)                                                                                                                          ggCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTAc                            >  2:811345/1‑149 (MQ=255)                                                                                                                            cACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCCGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCACACCCCCTAccc                          >  2:210282/1‑149 (MQ=255)                                                                                                                                    ttCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGa                  <  2:257737/149‑1 (MQ=255)                                                                                                                                      ccAGTACGTTG‑ACGACCTGCATCATGCCTTTGGAGATGTCCTTATCAATACCGGCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAGCTCGGTTCCTGTTCCCTACTCACTACCCGTTGTTGATg                <  2:643020/149‑1 (MQ=255)                                                                                                                                        aGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATggg              <  2:166839/149‑1 (MQ=255)                                                                                                                                              ttG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCa        <  1:811345/149‑1 (MQ=255)                                                                                                                                               tG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTCCCCGTTGTTGATGGGCGAACAg       >  2:23011/1‑149 (MQ=255)                                                                                                                                                  acgacCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTCTTTTTTTGTTTCCTGTTCCCCACTCACTATCGGTTGTTGATGGGCGAGCAggg     >  1:887809/1‑149 (MQ=255)                                                                                                                                                     acCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCgt  >  2:705885/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCCAGTTCTGCCAGCAACTCGTTATCCAGCGGCTCTATTTCTTCACCACGGCGACTGCGCATCAGTGCACTCGCCATATGTGTCATCGCCATCGTTCCTTGCTCACTGCGCACCGGCAGATGGCACTCTGTTTCCAGTACGTTG‑ACGACCTGCATCATGCCTTTGCAGATGTCCTTATCAATAACGCCTGCCTCGCAAAGCAGGTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGT  >  NC_000913/3511200‑3511495

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