BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,908,486	G→C	5.4%	A394G (GCC→GGC)	pyrG ←	CTP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,908,486	0	G	C	5.4%	388.8 / 10.8	112	A394G (GCC→GGC)	pyrG	CTP synthetase
Reads supporting (aligned to +/- strand): ref base G (55/51); new base C (4/2); total (59/53)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 114 total aligned reads displayed.

TGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGT  >  NC_000913/2908338‑2908633

tgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTgg                                                                                                                                                     <  3:510097/149‑1 (MQ=255)

 gcgcCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTgg                                                                                                                                                     <  2:677336/148‑1 (MQ=255)

    ccGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCaa                                                                                                                                                  <  2:462075/148‑1 (MQ=255)

     cGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACa                                                                                                                                                <  3:297203/149‑1 (MQ=255)

         aCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATgg                                                                                                                                             >  2:290329/1‑148 (MQ=255)

           gCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCg                                                                                                                                           <  2:538125/148‑1 (MQ=255)

                 tACCGCCGAGATCGCTCTTCTCGCTACGAGCTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTaa                                                                                                                                    >  3:721966/1‑149 (MQ=255)

                  aCCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTaa                                                                                                                                    >  2:570374/1‑148 (MQ=255)

                   ccgccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATc                                                                                                                                  <  3:25735/149‑1 (MQ=255)

                    cgccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCa                                                                                                                                 >  4:476845/1‑149 (MQ=255)

                     gccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCa                                                                                                                                 <  1:570374/148‑1 (MQ=255)

                      ccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCaa                                                                                                                                <  1:13331/148‑1 (MQ=255)

                         agaTCGCTCTTCTCGCTATGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTaa                                                                                                                            >  3:413358/1‑149 (MQ=255)

                          gaTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAAc                                                                                                                           <  4:170042/149‑1 (MQ=255)

                             cGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGcc                                                                                                                        >  4:304378/1‑149 (MQ=255)

                                tcttctCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCAcc                                                                                                                     >  3:910415/1‑149 (MQ=255)

                                 cttctCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGCTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCAcc                                                                                                                     >  1:659227/1‑148 (MQ=255)

                                   tctcGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTg                                                                                                                   <  2:468173/148‑1 (MQ=255)

                                       gCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCCGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATAc                                                                                                              <  3:801699/149‑1 (MQ=255)

                                         tACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATAccc                                                                                                            >  4:109953/1‑149 (MQ=255)

                                         tACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATAccc                                                                                                            >  4:628274/1‑149 (MQ=255)

                                          aCGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACTTGGCGAGCGTAATCAATTAACGCCACCTGCATAccc                                                                                                            >  2:377681/1‑148 (MQ=255)

                                           cGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAg                                                                                                          <  4:289219/149‑1 (MQ=255)

                                            gAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAgg                                                                                                         >  3:323513/1‑149 (MQ=255)

                                            gAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAgg                                                                                                         >  3:599984/1‑149 (MQ=255)

                                             aaCTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAgg                                                                                                         >  2:791894/1‑148 (MQ=255)

                                             aaCTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAgg                                                                                                         <  1:100485/148‑1 (MQ=255)

                                                 tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaaa                                                                                                     >  1:481317/1‑148 (MQ=255)

                                                 tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaaa                                                                                                     >  1:209303/1‑148 (MQ=255)

                                                 tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaaa                                                                                                     >  2:917311/1‑148 (MQ=255)

                                                 tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaaa                                                                                                     >  1:806419/1‑148 (MQ=255)

                                                 tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAAt                                                                                                    >  4:508728/1‑149 (MQ=255)

                                                  cAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAAt                                                                                                    >  1:53196/1‑148 (MQ=255)

                                                  cAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAAt                                                                                                    >  2:575348/1‑148 (MQ=255)

                                                  cAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATg                                                                                                   <  4:599984/149‑1 (MQ=255)

                                                    aCGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGc                                                                                                  >  4:874402/1‑148 (MQ=255)

                                                     cGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGccc                                                                                                <  4:244838/149‑1 (MQ=255)

                                                     cGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGcc                                                                                                 <  1:124510/148‑1 (MQ=255)

                                                            gTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                          <  1:43919/148‑1 (MQ=255)

                                                            gTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                          >  2:614036/1‑148 (MQ=255)

                                                            gTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCC‑CCTGCATACCCAGGCAAATGCCCAGATAAgg                                                                                        >  3:78207/1‑149 (MQ=255)

                                                             ttttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAgg                                                                                        >  3:742342/1‑149 (MQ=255)

                                                             ttttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAgg                                                                                        >  3:357436/1‑149 (MQ=255)

                                                             ttttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                         >  1:85171/1‑148 (MQ=255)

                                                             ttttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                         >  1:258624/1‑148 (MQ=255)

                                                                  aTCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAatat                                                                                   >  3:145718/1‑149 (MQ=255)

                                                                   tCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCGTACCCAGGCAAATGCCCAGATAAGGAatat                                                                                   <  1:245034/148‑1 (MQ=255)

                                                                   tCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAatat                                                                                   <  1:786698/148‑1 (MQ=255)

                                                                      cgcCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGt                                                                                <  1:575348/148‑1 (MQ=255)

                                                                       gcCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGtt                                                                               <  1:735264/148‑1 (MQ=255)

                                                                       gcCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTtc                                                                              >  4:107236/1‑149 (MQ=255)

                                                                          aCTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTtctc                                                                            >  1:674436/1‑148 (MQ=255)

                                                                           cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCa                                                                           >  2:41375/1‑148 (MQ=255)

                                                                           cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCa                                                                           >  1:442887/1‑148 (MQ=255)

                                                                                      gcgcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcg                                                                >  2:439415/1‑148 (MQ=255)

                                                                                      gcgcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcg                                                                >  1:776552/1‑148 (MQ=255)

                                                                                      gcgcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcg                                                                <  2:806419/148‑1 (MQ=255)

                                                                                        gcCACAACCGGGTACTTGCAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCGCGCGCAAAAcgcg                                                              <  1:729021/148‑1 (MQ=255)

                                                                                         ccACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgcgc                                                             >  2:759923/1‑148 (MQ=255)

                                                                                           acaACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCgg                                                           >  2:218592/1‑148 (MQ=255)

                                                                                                   gTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCa                                                  <  3:508728/149‑1 (MQ=255)

                                                                                                    tACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAATGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCa                                                  <  2:504630/148‑1 (MQ=255)

                                                                                                    tACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCAt                                                 >  3:475782/1‑149 (MQ=255)

                                                                                                     aCTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCAt                                                 <  2:674436/148‑1 (MQ=255)

                                                                                                      cTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGc                                               <  3:870392/149‑1 (MQ=255)

                                                                                                       ttACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGcc                                              <  4:145718/149‑1 (MQ=255)

                                                                                                       ttACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGcc                                              >  3:861505/1‑149 (MQ=255)

                                                                                                        tACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGcc                                              >  1:351135/1‑148 (MQ=255)

                                                                                                         aCAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCtt                                            <  3:476845/149‑1 (MQ=255)

                                                                                                          cAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGTCtt                                            <  2:175621/148‑1 (MQ=255)

                                                                                                            gTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCt                                          <  1:452522/148‑1 (MQ=255)

                                                                                                             tCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATTCCTTCTAc                                        >  4:438868/1‑149 (MQ=255)

                                                                                                             tCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTAc                                        <  4:798288/149‑1 (MQ=255)

                                                                                                                ggCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAGGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGcc                                     >  3:622183/1‑149 (MQ=255)

                                                                                                                  cacaAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCAc                                   <  3:107236/149‑1 (MQ=255)

                                                                                                                    caAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcg                                                                >  2:18413/1‑118 (MQ=255)

                                                                                                                    caAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcg                                                                <  1:18413/118‑1 (MQ=255)

                                                                                                                    caAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCtt                                            >  1:481928/1‑138 (MQ=255)

                                                                                                                    caAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACg                                  >  1:524243/1‑148 (MQ=255)

                                                                                                                    caAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACg                                  >  1:625667/1‑148 (MQ=255)

                                                                                                                            gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATa                               <  4:127324/143‑1 (MQ=255)

                                                                                                                            gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATa                               >  3:127324/1‑143 (MQ=255)

                                                                                                                            gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGaa                         >  4:907723/1‑149 (MQ=255)

                                                                                                                                   tGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGcc                   >  1:923418/1‑148 (MQ=255)

                                                                                                                                      cGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGt               <  4:357436/149‑1 (MQ=255)

                                                                                                                                        ttCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa              <  2:94048/148‑1 (MQ=255)

                                                                                                                                        ttCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTAc             >  4:176675/1‑149 (MQ=255)

                                                                                                                                         tcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg            >  3:118505/4‑149 (MQ=255)

                                                                                                                                         tcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg            >  3:350004/4‑149 (MQ=255)

                                                                                                                                         tcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAACCATGCCTTCTACGCCACGATAGcc                            >  3:696651/4‑133 (MQ=255)

                                                                                                                                         tcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAACCATGCCTTCTACGCCACGATAGcc                            <  4:696651/130‑1 (MQ=255)

                                                                                                                                          ctcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg            <  1:439415/148‑1 (MQ=255)

                                                                                                                                          ctcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg            <  1:319593/148‑1 (MQ=255)

                                                                                                                                          ctcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg            <  1:233727/148‑1 (MQ=255)

                                                                                                                                          cgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg            >  2:384585/3‑148 (MQ=255)

                                                                                                                                              aTGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGa        <  2:231660/148‑1 (MQ=255)

                                                                                                                                                gTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGAtt      >  2:378198/1‑148 (MQ=255)

                                                                                                                                                 ttGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc    <  3:589384/149‑1 (MQ=255)

                                                                                                                                                  tGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc    <  1:430169/148‑1 (MQ=255)

                                                                                                                                                  tGGCAACATGGCGAGCGTAATCAATTAACGCC‑CCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAAGCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGt  <  4:78207/149‑1 (MQ=255)

TGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGT  >  NC_000913/2908338‑2908633

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