BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,995,722	T→G	5.1%	V47G (GTA→GGA)	yigA →	DUF484 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,722	0	T	G	5.1%	356.5 / 7.2	118	V47G (GTA→GGA)	yigA	DUF484 family protein
Reads supporting (aligned to +/- strand): ref base T (65/47); new base G (3/3); total (68/50)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 118 total aligned reads displayed.

TATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTC  >  NC_000913/3995574‑3995868

tATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTa                                                                                                                                                   >  2:76736/1‑150 (MQ=255)

     aTCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCgg                                                                                                                                              <  2:612360/150‑1 (MQ=255)

       cTATGAAGCACCCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTCTCCGTAATGCGCGCGCAGTAGAAGAGATACGTGTGCCTCATCCGGTACGCGGCa                                                                                                                                            >  1:300093/1‑150 (MQ=255)

         aTGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCGTACCGGATGCGGCAca                                                                                                                                          <  2:1892905/150‑2 (MQ=255)

          tGAAGCACCCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCGTACCGGATGCGGCAcac                                                                                                                                         <  1:998615/150‑3 (MQ=255)

          tGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCg                                                                                                                                         >  1:321134/1‑150 (MQ=255)

          tGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCg                                                                                                                                         >  2:594148/1‑150 (MQ=255)

          tGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGGCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCg                                                                                                                                         <  1:1642098/150‑1 (MQ=255)

          tGAAGCAACCAGCGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCACTAGAAGCGAAACGTCTGCCGCATCCGGTACGCGCCAccc                                                                                                                                         >  1:1327905/1‑149 (MQ=255)

           gAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGt                                                                                                                                        >  1:593428/1‑150 (MQ=255)

            aaGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCGTACCGGATGCGGCAcacgt                                                                                                                                       <  1:1709631/150‑5 (MQ=255)

             aGCAACCAGGGGAGGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGttt                                                                                                                                      >  2:570352/1‑150 (MQ=255)

             aGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGttt                                                                                                                                      >  2:1321291/1‑150 (MQ=255)

              gCAACCCGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGCCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTc                                                                                                                                     <  1:1713529/150‑1 (MQ=255)

                  ccAGGGGAAGACCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTg                                                                                                                                 <  2:583705/150‑1 (MQ=255)

                   cAGGGGAAGAACTGCAGGCAACACTCAAGGAGCTTGATGACCGGGCGGTTGTCGATTATCTCATTAAAAATCCTGAGTTTTTTAACCTTAATGCGCGCGCAGTAGAAGCGATACGTGAGACGCATCCGGTACGCGGCACCGTTTCGTTgg                                                                                                                                >  1:1546535/1‑150 (MQ=255)

                   cAGGGGAAGAACTGCAGGAAACACTCACGGAGCGTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTACCGCATCCGGTACGCGGCACCGTTTCGTTgg                                                                                                                                <  2:524593/150‑1 (MQ=255)

                         aagacCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGGGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGt                                                                                                                          <  1:628853/150‑1 (MQ=255)

                         aactaCTGCAGGAAACACTCACGGAGGTTGAAGACCGGGCGGTTGTCGATTATCTGATTAAAACTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGt                                                                                                                          <  1:1633502/146‑1 (MQ=255)

                          agaaCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTg                                                                                                                         >  2:1734696/1‑150 (MQ=255)

                           gaaCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTgg                                                                                                                        <  1:409896/150‑1 (MQ=255)

                             aCTGCAGGAACCACTCACGGCGCTTGATGACCGGGCGGTTGTCGAGTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGca                                                                                                                      <  2:1006665/150‑1 (MQ=255)

                               tGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGcaca                                                                                                                    >  2:45398/1‑150 (MQ=255)

                                gCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACAt                                                                                                                   >  2:783847/1‑150 (MQ=255)

                                gCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACAt                                                                                                                   <  1:1304232/150‑1 (MQ=255)

                                gCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACAt                                                                                                                   >  2:1894560/1‑150 (MQ=255)

                                 cAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATg                                                                                                                  <  1:594148/150‑1 (MQ=255)

                                   ggAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGc                                                                                                                >  2:1749240/1‑150 (MQ=255)

                                      aaCACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCTGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCcg                                                                                                             <  2:512444/150‑1 (MQ=255)

                                        cacTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCcgcg                                                                                                           >  1:1496114/1‑150 (MQ=255)

                                          cTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCTTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGTTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCa                                                                                                         >  1:169467/1‑150 (MQ=255)

                                          cTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCa                                                                                                         <  1:168862/150‑1 (MQ=255)

                                             aCGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGt                                                                                                      >  1:861909/1‑150 (MQ=255)

                                               ggAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGAACATGGCCCGCTCACGTaa                                                                                                    >  1:439158/1‑150 (MQ=255)

                                                gAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAAt                                                                                                   <  1:1614579/150‑1 (MQ=255)

                                                gAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAAt                                                                                                   <  1:744332/150‑1 (MQ=255)

                                                  gCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCa                                                                                                 >  1:1008055/1‑150 (MQ=255)

                                                    ttGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCata                                                                                               >  1:22246/1‑150 (MQ=255)

                                                       atgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTc                                                                                            <  1:534220/150‑1 (MQ=255)

                                                         gaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCAt                                                                                          <  2:1531012/150‑1 (MQ=255)

                                                              ggCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCt                                                                                     <  1:959193/150‑1 (MQ=255)

                                                               gCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTg                                                                                    <  1:182598/150‑1 (MQ=255)

                                                               gCGGTTGTCGAGTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTg                                                                                    <  1:156541/150‑1 (MQ=255)

                                                                  gTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGaa                                                                                 >  1:48649/1‑150 (MQ=255)

                                                                   ttGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATTATATTCATGTTCTGGAag                                                                                <  1:76736/150‑1 (MQ=255)

                                                                       cGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGaa                                                                            <  1:1900317/150‑1 (MQ=255)

                                                                       cGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGaa                                                                            >  1:1202423/1‑150 (MQ=255)

                                                                          ttATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACAt                                                                         >  1:645918/1‑150 (MQ=255)

                                                                              cTGATTAAAAATCCTGAGTATTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACAGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAACAGAACATGgcg                                                                     >  2:290042/1‑150 (MQ=255)

                                                                               tGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgcgc                                                                    >  1:1720201/1‑150 (MQ=255)

                                                                               tGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgcgc                                                                    >  1:1343882/1‑150 (MQ=255)

                                                                                gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCt                                                                   >  2:1058988/1‑150 (MQ=255)

                                                                                gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCt                                                                   >  1:683407/1‑150 (MQ=255)

                                                                                gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCt                                                                   <  1:784889/150‑1 (MQ=255)

                                                                                gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGCAAGAGAACATGgcgat                                                                   >  1:1779576/1‑148 (MQ=255)

                                                                                  ttAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGGTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCtct                                                                 >  2:1716350/1‑148 (MQ=255)

                                                                                     aaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGa                                                              <  2:1564059/150‑1 (MQ=255)

                                                                                     aaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGa                                                              >  2:26435/1‑150 (MQ=255)

                                                                                     aaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGa                                                              >  2:1317026/1‑150 (MQ=255)

                                                                                         tCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGa                                                          <  1:386864/150‑1 (MQ=255)

                                                                                           cTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGATGAGAACATGGCGCTGTTGATGGAAc                                                        >  1:1463905/1‑150 (MQ=255)

                                                                                           cTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAAc                                                        <  2:780914/150‑1 (MQ=255)

                                                                                            tGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACa                                                       <  2:562072/150‑1 (MQ=255)

                                                                                             gAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAg                                                      <  2:17918/150‑1 (MQ=255)

                                                                                               gTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGc                                                    <  2:594648/150‑1 (MQ=255)

                                                                                               gTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGGGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGc                                                    <  1:370511/150‑1 (MQ=255)

                                                                                                  ttttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTAt                                                 >  1:429267/1‑150 (MQ=255)

                                                                                                   tttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATc                                                >  2:1489239/1‑150 (MQ=255)

                                                                                                    ttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCg                                               <  2:1494027/150‑1 (MQ=255)

                                                                                                         cgtcaggcgcgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAAc                                          <  1:1383278/144‑1 (MQ=255)

                                                                                                         cgaaaTGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAAc                                          <  1:1447907/147‑1 (MQ=255)

                                                                                                         cGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAAc                                          <  2:189366/150‑1 (MQ=255)

                                                                                                                cgcgcgcAGTAGAAGCGATCCGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGcc                                   <  1:1833089/150‑1 (MQ=255)

                                                                                                                cgcgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGcc                                   <  1:866839/150‑1 (MQ=255)

                                                                                                                 gcgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGATCAGGCTATCGCCAACGAAGGCCt                                  >  2:143784/1‑150 (MQ=255)

                                                                                                                  cgcgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTg                                 >  2:1325285/1‑150 (MQ=255)

                                                                                                                    cgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGtt                               >  2:1514557/1‑150 (MQ=255)

                                                                                                                     gcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGGTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGtgt                              >  1:1085956/1‑148 (MQ=255)

                                                                                                                       aGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCTAGTGTCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTa                            <  2:774677/150‑1 (MQ=255)

                                                                                                                           gAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGc                        >  1:1231192/1‑150 (MQ=255)

                                                                                                                           gAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGc                        >  1:806266/1‑150 (MQ=255)

                                                                                                                            aaGCGATACGTGTTCCGCATCCGGTACGCGTCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTTATCATATTCATGTTCTGGAAGAGAAAATGGCGTTGTTGTTGTAACAGGCTTTCGCCAACGAAGGCCTGTTTTATCGcc                       <  2:1794889/150‑1 (MQ=255)

                                                                                                                            aaGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGcc                       >  1:925180/1‑150 (MQ=255)

                                                                                                                            aaGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGcc                       <  1:1342132/150‑1 (MQ=255)

                                                                                                                            aaGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGcc                       <  1:1448250/150‑1 (MQ=255)

                                                                                                                             aGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCt                      <  1:570352/150‑1 (MQ=255)

                                                                                                                               cGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTAc                    >  2:859670/1‑150 (MQ=255)

                                                                                                                                gATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTAct                   >  1:1254853/1‑150 (MQ=255)

                                                                                                                                     cggtgCCGCGTACCGGATGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTAcc              >  1:1385799/3‑150 (MQ=255)

                                                                                                                                      tgtgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt             >  1:707070/1‑150 (MQ=255)

                                                                                                                                      ggtgCCGCGTACCGGATGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt             >  1:146567/2‑150 (MQ=255)

                                                                                                                                        tgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGc           <  1:519063/150‑1 (MQ=255)

                                                                                                                                         gCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCa          >  2:650019/1‑150 (MQ=255)

                                                                                                                                         gCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCAGCGCACGTAATCATATTCATGTTCTGGAAGAGAACTTGGCGGTGTTCATGGAACAGGCTATCGCCAACGCAGGGCTGTTTTATCGCATACTCTACCTGcc          >  2:1862371/1‑149 (MQ=255)

                                                                                                                                         gCCGCATCCGGTACGCGGCACCGTTTCGATGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCa          >  2:667897/1‑150 (MQ=255)

                                                                                                                                            gCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTATGGAAGAGAACATTGCGCTGGTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATAGCCTACTCTACCTGCAgcg       >  2:1455030/1‑150 (MQ=255)

                                                                                                                                            gCATCCGGTACGCGGCACCGGTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGATGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCCACGAAGGCCTGTGTTTTCGCCTACTCTACCTgcggcg       >  1:302472/1‑150 (MQ=255)

                                                                                                                                                 cGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCCCGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCAACCTGCAGCGCAGtc  >  1:157823/1‑150 (MQ=255)

                                                                                                                                                 cGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGAGGTAACAGTCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGtc  <  1:1862696/150‑1 (MQ=255)

                                                                                                                                                 cGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGATCTGGAAGAGAACATTGCGCTGTAGATGGAACAGGCTATCGCCAACGAAGCCCTGATTAATCGCCTACTCTACCTGCCGCGCAGta  >  1:1167181/1‑149 (MQ=255)

TATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTC  >  NC_000913/3995574‑3995868

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`