BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,831,768	A→G	5.3%	T438A (ACA→GCA)	yicH →	AsmA famly protein YicH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,831,768	0	A	G	5.3%	890.9 / 24.2	284	T438A (ACA→GCA)	yicH	AsmA famly protein YicH
Reads supporting (aligned to +/- strand): ref base A (128/141); new base G (11/4); total (139/145)
Fisher's exact test for biased strand distribution p-value = 6.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 285 total aligned reads displayed.

GCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACAC  >  NC_000913/3831619‑3831915

   ccAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAt                                                                                                                                                  >  1:1049305/1‑150 (MQ=255)

   ccAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAt                                                                                                                                                  <  2:1779495/150‑1 (MQ=255)

    cAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAtt                                                                                                                                                 >  1:1138698/1‑150 (MQ=255)

    cAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAtt                                                                                                                                                 <  2:1558928/150‑1 (MQ=255)

     aGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTc                                                                                                                                                <  2:1272759/150‑1 (MQ=255)

           aaTCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattc                                                                                                                                          <  1:350260/150‑6 (MQ=255)

           aaTCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCCGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgt                                                                                                                                          >  1:1409619/1‑150 (MQ=255)

            aTCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtg                                                                                                                                         <  1:1795734/150‑1 (MQ=255)

               tGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                                                                      >  1:394288/1‑150 (MQ=255)

               tGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCCGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                                                                      <  2:412520/150‑1 (MQ=255)

                 aTCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGAt                                                                                                                                    <  1:1683198/150‑1 (MQ=255)

                   cATTGATATCGCCCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCAGCCGCCGCGACATTCAATCGTGTTGATGt                                                                                                                                  >  2:1966704/1‑150 (MQ=255)

                     ttGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtc                                                                                                                                <  1:1524760/150‑1 (MQ=255)

                     ttGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtc                                                                                                                                >  2:643658/1‑150 (MQ=255)

                       gATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgt                                                                                                                              >  2:1735579/1‑150 (MQ=255)

                           tCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGcc                                                                                                                          <  2:809531/150‑1 (MQ=255)

                            cGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGccc                                                                                                                         >  2:900477/1‑150 (MQ=255)

                             gACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                                                                                        >  2:848537/1‑150 (MQ=255)

                             gACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                                                                                        >  1:1770110/1‑150 (MQ=255)

                              aCCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGt                                                                                                                       >  2:1820689/1‑150 (MQ=255)

                                 cTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                                                                                    <  1:92487/150‑1 (MQ=255)

                                    aCTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                                                                                 <  1:1293042/150‑1 (MQ=255)

                                     cTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgc                                                                                                                <  2:1410712/150‑1 (MQ=255)

                                       tCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGccc                                                                                                                         >  1:463760/1‑139 (MQ=255)

                                       tCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgcgc                                                                                                              <  1:1822383/150‑1 (MQ=255)

                                        cccGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCt                                                                                                             >  2:185750/1‑150 (MQ=255)

                                          cGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                                                                      <  2:472100/123‑1 (MQ=255)

                                           gTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                                                                                          <  2:1924066/150‑1 (MQ=255)

                                             ggCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                                                                        >  2:1300613/1‑150 (MQ=255)

                                                 gCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCaa                                                                                                    >  2:554878/1‑150 (MQ=255)

                                                  cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGGTCAGCGCCagc                                                                                                   <  2:263486/150‑3 (MQ=255)

                                                  cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                                                                                   <  1:382139/150‑1 (MQ=255)

                                                    cACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                                                                                 >  2:1008421/1‑150 (MQ=255)

                                                     aCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCt                                                                                                             <  2:1449880/137‑1 (MQ=255)

                                                     aCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                                                                <  2:1373932/150‑1 (MQ=255)

                                                     aCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagc                                                                                                >  1:839957/1‑150 (MQ=255)

                                                      ccGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                                                                               >  2:773682/1‑150 (MQ=255)

                                                      ccGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCCGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                                                                               <  1:1505698/150‑1 (MQ=255)

                                                            tCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                                                                         <  2:1222530/150‑1 (MQ=255)

                                                             cGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGt                                                                                        <  2:1157469/150‑1 (MQ=255)

                                                                         tGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGa                                                                                      <  1:696405/140‑1 (MQ=255)

                                                                           ccAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGCAc                                                                          >  2:711057/1‑150 (MQ=255)

                                                                           ccAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                                                                          >  1:943062/1‑150 (MQ=255)

                                                                                cTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCCGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                                                                     <  2:1082009/150‑1 (MQ=255)

                                                                                 tGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                                                                    >  1:1249035/1‑150 (MQ=255)

                                                                                  gACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc                                                                   >  2:1668505/1‑150 (MQ=255)

                                                                                   aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa                                                                  >  1:309134/1‑150 (MQ=255)

                                                                                    cGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAt                                                                 <  2:1039245/150‑1 (MQ=255)

                                                                                      cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAttt                                                               <  1:1325360/150‑1 (MQ=255)

                                                                                      cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACCGTGAATATCAGCGAACTGAGTGCAttt                                                               >  2:1831222/1‑150 (MQ=255)

                                                                                        ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                             <  2:257303/150‑1 (MQ=255)

                                                                                        ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGGCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc                                                             <  2:327518/150‑1 (MQ=255)

                                                                                          ttACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTg                                                           <  2:22488/150‑1 (MQ=255)

                                                                                           tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGa                                                          <  2:924020/150‑1 (MQ=255)

                                                                                                aTCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAg                                                     <  1:1668505/150‑1 (MQ=255)

                                                                                                 tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAgg                                                    <  1:1983808/150‑1 (MQ=255)

                                                                                                  cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGc                                                   >  2:1190754/1‑150 (MQ=255)

                                                                                                   aTAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCa                                                  >  1:241003/1‑150 (MQ=255)

                                                                                                     aaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAtt                                                >  2:1974803/1‑150 (MQ=255)

                                                                                                     aaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGACGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCAtt                                                >  2:591701/1‑150 (MQ=255)

                                                                                                      aaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTc                                               <  1:1190754/150‑1 (MQ=255)

                                                                                                        tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCTTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCTAACTGAGTGCATTTACTGAAAAAGGCATTCTg                                             <  1:586579/150‑1 (MQ=255)

                                                                                                        tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGGGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTg                                             <  2:1936534/150‑1 (MQ=255)

                                                                                                         ggggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTgg                                            <  1:1706549/150‑1 (MQ=255)

                                                                                                           ggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg                                                                      >  1:896183/1‑122 (MQ=255)

                                                                                                           ggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGaa                                          >  2:870432/1‑150 (MQ=255)

                                                                                                            gCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAg                                         <  2:1596738/150‑1 (MQ=255)

                                                                                                                 tGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCAcc                                    <  1:1618172/150‑1 (MQ=255)

                                                                                                                    aGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGcc                                 <  1:554878/150‑1 (MQ=255)

                                                                                                                    aGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGcc                                 <  1:1592051/150‑1 (MQ=255)

                                                                                                                      tGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAg                               <  2:1520077/150‑1 (MQ=255)

                                                                                                                       ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgt                              <  1:1441913/150‑1 (MQ=255)

                                                                                                                       ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgt                              >  1:508804/1‑150 (MQ=255)

                                                                                                                       ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCAGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgt                              >  2:1474140/1‑150 (MQ=255)

                                                                                                                        gCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtg                             >  2:573943/1‑150 (MQ=255)

                                                                                                                         cTCGGCGAATCTGAATGCTGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtgt                            >  2:1579430/1‑150 (MQ=255)

                                                                                                                             gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTca                        <  1:1582694/150‑1 (MQ=255)

                                                                                                                              cGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTcac                       <  1:1920872/150‑1 (MQ=255)

                                                                                                                              cGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTcac                       <  2:499981/150‑1 (MQ=255)

                                                                                                                               gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGATGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTcaca                      <  1:1948270/150‑1 (MQ=255)

                                                                                                                                cgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACaa                     >  1:663841/5‑150 (MQ=255)

                                                                                                                                cgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACaa                     >  1:1469299/5‑150 (MQ=255)

                                                                                                                                aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACaa                     >  2:1757220/1‑150 (MQ=255)

                                                                                                                                  tCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAAc                   <  2:748694/150‑1 (MQ=255)

                                                                                                                                  tCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAAc                   <  2:450628/150‑1 (MQ=255)

                                                                                                                                  attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAAc                   >  1:439869/3‑150 (MQ=255)

                                                                                                                                  attGAATGTCGCGGCGGCGGCATTAAATCCAGTTGATGTTCGTCGCCAGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTCCATTTACTGAAAAAGGCAATCTGGAAGCCCCCGCCAGTGTTTCACAAAc                   >  2:1618342/3‑150 (MQ=255)

                                                                                                                                   gtGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACg                  <  1:199579/149‑1 (MQ=255)

                                                                                                                                   cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACg                  <  1:1218128/150‑1 (MQ=255)

                                                                                                                                   cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACg                  <  1:773682/150‑1 (MQ=255)

                                                                                                                                     gAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGcc                <  2:369144/150‑1 (MQ=255)

                                                                                                                                     gAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGcc                <  2:285323/150‑1 (MQ=255)

                                                                                                                                         gccgccgccGCGACATTCAATCGTGTTGATGTTCGTCGACCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACaa            >  2:1325658/1‑150 (MQ=255)

                                                                                                                                          ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAgtg                             >  2:1785493/1‑132 (MQ=255)

                                                                                                                                           cgccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACg          <  1:219229/150‑1 (MQ=255)

                                                                                                                                           cgccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACg          >  1:457713/1‑150 (MQ=255)

                                                                                                                                            gccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGt         <  1:1015893/150‑1 (MQ=255)

                                                                                                                                              cgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGATCGGAGTGGATTTACTGATAAAGGCTTTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCa       >  2:623864/1‑150 (MQ=255)

                                                                                                                                                 cgcgACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGAATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGac    >  2:1710499/1‑150 (MQ=255)

                                                                                                                                                  gcgACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGaca   <  1:47326/150‑1 (MQ=255)

GCGCCAGCCGCAATCTGATCATTGATATCGACCCTGACTTCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACCGCCAGTGTTTCACAAACGCCACAACGTCAGACAC  >  NC_000913/3831619‑3831915

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