BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,995,724	C→T	5.1%	R48C (CGC→TGC)	yigA →	DUF484 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,724	0	C	T	5.1%	375.7 / 4.9	117	R48C (CGC→TGC)	yigA	DUF484 family protein
Reads supporting (aligned to +/- strand): ref base C (63/47); new base T (3/3); total (67/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 117 total aligned reads displayed.

ATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTC  >  NC_000913/3995579‑3995868

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)

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

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

              cAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTgg                                                                                                                                >  2:1697734/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                                                                                                                   >  1:230913/1‑150 (MQ=255)

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

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

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

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

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

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

                                    ccTCACGGAGCGTGATGACCGGGCGGGTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCgcgc                                                                                                          <  1:1653208/149‑1 (MQ=255)

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

                                     cTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCa                                                                                                         >  1:837631/1‑150 (MQ=255)

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

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

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

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

                                                tgatgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCatat                                                                                              >  2:572682/1‑150 (MQ=255)

                                                 gatgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATAtt                                                                                             <  2:1012521/150‑1 (MQ=255)

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

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

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

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

                                                            ggTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGa                                                                                  >  2:68651/1‑150 (MQ=255)

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

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

                                                              ttGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAag                                                                                <  1:1734882/150‑1 (MQ=255)

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

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

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

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

                                                                         cTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgcg                                                                     >  1:397229/1‑150 (MQ=255)

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

                                                                          tGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACTTGTGCCGCATCCGGTAAGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgcgc                                                                    >  1:651739/1‑150 (MQ=255)

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

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

                                                                           gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGCACATGGCGCt                                                                   >  2:144171/1‑150 (MQ=255)

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

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

                                                                           gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCt                                                                   >  2:1058988/1‑150 (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: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)

                                                                                          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)

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

                                                                                                                                 tgtgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGCACATGGCGCGGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt             >  2:312387/1‑150 (MQ=255)

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

                                                                                                                                 ggtgCCGCGTACCGGATGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt             >  1:554983/2‑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)

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

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

                                                                                                                                            cGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGtc  <  2:1825513/150‑1 (MQ=255)

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

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

ATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTC  >  NC_000913/3995579‑3995868

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