Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,799,867	C→A	R43L (CGA→CTA)	ygaC ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,799,867	0	C	A	98.0%	383.2 / ‑1.0	98	R43L (CGA→CTA)	ygaC	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (51/45);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTCGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCTCGTTATTATCAAC  >  NC_000913/2799719‑2800017                                                                                                                                                     |                                                                                                                                                       gATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTagc                                                                                                                                                      <  2:770651/151‑1 (MQ=255)      cGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTc                                                                                                                                                 <  3:1388328/151‑1 (MQ=255)      cGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTc                                                                                                                                                 >  3:1093972/1‑151 (MQ=255)        tAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTcgc                                                                                                                                               >  2:329691/1‑151 (MQ=255)        tAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTcgc                                                                                                                                               >  3:804802/1‑151 (MQ=255)         aaTGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTcgc                                                                                                                                               <  2:1040322/150‑1 (MQ=255)         aaTGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTcgc                                                                                                                                               <  2:106593/150‑1 (MQ=255)              tCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTa                                                                                                                                          >  3:417873/1‑150 (MQ=255)              tCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTa                                                                                                                                          >  2:532627/1‑150 (MQ=255)               cGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTaa                                                                                                                                         <  3:428638/150‑1 (MQ=255)               cGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTaa                                                                                                                                         <  2:346339/150‑1 (MQ=255)                  gtCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACAt                                                                                                                                      >  2:646989/1‑150 (MQ=255)                     gctcgcCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGAc                                                                                                                                  >  2:919302/1‑151 (MQ=255)                      ctcgcCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGAc                                                                                                                                  >  3:271917/1‑150 (MQ=255)                       tcgcCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACAt                                                                                                                                >  3:382277/1‑151 (MQ=255)                       tcgcCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACAt                                                                                                                                >  3:64754/1‑151 (MQ=255)                       tcgcCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACAt                                                                                                                                >  2:343672/1‑151 (MQ=255)                        cgcCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACAt                                                                                                                                <  3:1539736/150‑1 (MQ=255)                             gCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAAtttt                                                                                                                          >  3:890702/1‑151 (MQ=255)                             gCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAAtttt                                                                                                                          >  2:1053669/1‑151 (MQ=255)                              cTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTc                                                                                                                         >  2:538954/1‑151 (MQ=255)                                aaaTAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGTCATAATTTTccc                                                                                                                       >  3:389522/1‑151 (MQ=255)                                  aTAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTcccc                                                                                                                      >  3:912068/1‑150 (MQ=255)                                       gCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCAcggcg                                                                                                                <  3:868911/151‑1 (MQ=255)                                        cGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCAcggcgg                                                                                                               >  2:520256/1‑151 (MQ=255)                                        cGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCAcggcgg                                                                                                               <  2:1057149/151‑1 (MQ=255)                                         ggAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGt                                                                                                              <  2:1550547/151‑1 (MQ=255)                                          gAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTa                                                                                                             <  2:1137397/151‑1 (MQ=255)                                            aaTTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAAc                                                                                                           <  2:1063102/151‑1 (MQ=255)                                             aTTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAAcc                                                                                                          <  2:1161300/151‑1 (MQ=255)                                              ttGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAAcc                                                                                                          >  3:895072/1‑150 (MQ=255)                                               tGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCa                                                                                                         >  3:1531527/1‑150 (MQ=255)                                                gCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCAt                                                                                                        >  3:414831/1‑150 (MQ=255)                                                 cTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATa                                                                                                       >  3:1510819/1‑150 (MQ=255)                                                    gTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCt                                                                                                   >  2:1273320/1‑151 (MQ=255)                                                     tAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCtt                                                                                                  >  2:1114088/1‑151 (MQ=255)                                                       aTGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTc                                                                                                <  2:935801/151‑1 (MQ=255)                                                         gATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCt                                                                                               <  3:1593193/150‑1 (MQ=255)                                                          aTCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTgg                                                                                             >  2:1629358/1‑151 (MQ=255)                                                             gCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGtt                                                                                          >  2:387596/1‑151 (MQ=255)                                                              cAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGtt                                                                                          >  3:796796/1‑150 (MQ=255)                                                               aaGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTa                                                                                         <  3:1499211/150‑1 (MQ=255)                                                               aaGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTa                                                                                         >  2:1499361/1‑150 (MQ=255)                                                               aaGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTAc                                                                                        <  3:1121872/151‑1 (MQ=255)                                                                aGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTAcc                                                                                       >  3:1549469/1‑151 (MQ=255)                                                                 gTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACca                                                                                      <  3:1646049/151‑1 (MQ=255)                                                                  ttttAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACca                                                                                      >  3:895906/1‑150 (MQ=255)                                                                    ttAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACcaca                                                                                    <  3:9259/150‑1 (MQ=255)                                                                     tAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACAt                                                                                   <  2:1200382/150‑1 (MQ=255)                                                                      aaTATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCg                                                                                 <  2:750145/151‑1 (MQ=255)                                                                       atatCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCg                                                                                 >  2:573685/1‑150 (MQ=255)                                                                         atCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGAc                                                                               >  2:355018/1‑150 (MQ=255)                                                                         atCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGAc                                                                               >  2:1060636/1‑150 (MQ=255)                                                                           cGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAg                                                                             <  2:1261513/150‑1 (MQ=255)                                                                           cGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGt                                                                            <  3:182649/151‑1 (MQ=255)                                                                              aaGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTaaa                                                                         <  3:496639/151‑1 (MQ=255)                                                                              aaGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTaa                                                                          >  2:560070/1‑150 (MQ=255)                                                                                 cGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAAcc                                                                      >  2:610334/1‑151 (MQ=255)                                                                                    gTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGAc                                                                   >  2:735367/1‑151 (MQ=255)                                                                                        aGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACtttt                                                               >  3:947655/1‑151 (MQ=255)                                                                                          ccAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTc                                                             <  3:41520/151‑1 (MQ=255)                                                                                             gtgtCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAg                                                          <  3:978221/151‑1 (MQ=255)                                                                                             gtgtCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGCTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGCCTTTTTCAAg                                                          <  3:38810/151‑1 (MQ=255)                                                                                               gtCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGt                                                         <  2:251403/150‑1 (MQ=255)                                                                                                 cGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTAcg                                                      >  2:861541/1‑151 (MQ=255)                                                                                                 cGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTAc                                                       <  3:1528009/150‑1 (MQ=255)                                                                                                        cGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCa                                               <  3:1603672/151‑1 (MQ=255)                                                                                                          ttCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCAcc                                             >  3:142192/1‑151 (MQ=255)                                                                                                            cTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACACAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTc                                           >  2:1030022/1‑151 (MQ=255)                                                                                                              tttAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTc                                           <  3:532591/149‑1 (MQ=255)                                                                                                                tAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCTCGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTc                                        >  3:1272970/1‑150 (MQ=255)                                                                                                                 aaTGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTg                                      <  2:418187/151‑1 (MQ=255)                                                                                                                  aTGTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTgg                                     >  2:1036677/1‑151 (MQ=255)                                                                                                                    gTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTc                                   >  3:1614813/1‑151 (MQ=255)                                                                                                                    gTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTc                                   <  3:1522164/151‑1 (MQ=255)                                                                                                                    gTTGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTc                                   <  3:862884/151‑1 (MQ=255)                                                                                                                     ttGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCg                                  >  2:318780/1‑151 (MQ=255)                                                                                                                      tGGATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCg                                  <  4:10374/150‑1 (MQ=255)                                                                                                                        gATGAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTaa                               <  2:632640/151‑1 (MQ=255)                                                                                                                           gAATTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAAt                             >  2:86201/1‑150 (MQ=255)                                                                                                                            aaTTACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATAc                           >  2:1040627/1‑151 (MQ=255)                                                                                                                               tACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACat                         >  2:609923/1‑150 (MQ=255)                                                                                                                               tACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACat                         <  3:637939/150‑1 (MQ=255)                                                                                                                               tACCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACat                         >  3:961511/1‑150 (MQ=255)                                                                                                                                aCCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACatat                       <  2:1519644/151‑1 (MQ=255)                                                                                                                                aCCAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACatat                       <  3:72958/151‑1 (MQ=255)                                                                                                                                 ccAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACatat                       <  3:1601322/150‑1 (MQ=255)                                                                                                                                 ccAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATAtt                      >  2:1410313/1‑151 (MQ=255)                                                                                                                                 ccAGAGCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATAtt                      <  3:1098753/151‑1 (MQ=255)                                                                                                                                     agCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGa                   >  4:6672/1‑150 (MQ=255)                                                                                                                                     agCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGAt                  <  3:621258/151‑1 (MQ=255)                                                                                                                                     agCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGAt                  >  2:1257525/1‑151 (MQ=255)                                                                                                                                     agCGGTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGAt                  <  2:1220332/151‑1 (MQ=255)                                                                                                                                        ggTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCt               <  3:329665/151‑1 (MQ=255)                                                                                                                                        ggTACGTCCCATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCt               <  3:1024742/151‑1 (MQ=255)                                                                                                                                           aCGTCCCATTCGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCTCGt            >  2:425370/1‑151 (MQ=255)                                                                                                                                               cccATTCGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCTCGttatt        <  3:1306201/151‑1 (MQ=255)                                                                                                                                               cccATTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCTCGttatt        >  2:785767/1‑151 (MQ=255)                                                                                                                                                  aTTAGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCTCGTTATTATc     <  3:1048714/151‑1 (MQ=255)                                                                                                                                                     cgcgCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCTCGTTATTATCAAc  <  2:731277/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCAAGTTTTAATATCGGAAGCGGGTTCAGCCAGTGTCGAACTGCGTTCTTTTAATGTTGGATGAATTACCAGAGCGGTACGTCCCATTCGCGCTTCGCGATTAACATAGACATAATTTTCCCCACGGCGGTAACCATACGCTTTCTGGGTTACCACATCGACAGTAAAACCGACTTTTTCAAGTACGCGCGCCACCTCGTCTGGTCGTAAATACATATTTGATCCTCGTTATTATCAAC  >  NC_000913/2799719‑2800017

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