BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,933,724	G→A	6.5%	C125Y (TGT→TAT) ‡	rbsD →	D‑ribose pyranase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,933,724	0	G	A	6.5%	379.0 / 9.8	122	C125Y (TGT→TAT) ‡	rbsD	D‑ribose pyranase
Reads supporting (aligned to +/- strand): ref base G (62/52); new base A (4/4); total (66/56)
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.52e-01

Only 100 of 125 total aligned reads displayed.

CCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTG  >  NC_000913/3933575‑3933873

ccATAATGCGCAACTCCACTAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATg                                                                                                                                                       <  2:340936/150‑1 (MQ=255)

   tAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTtc                                                                                                                                                    <  2:1302343/150‑1 (MQ=255)

       ccGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACCCCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCg                                                                                                                                                >  1:755087/1‑150 (MQ=255)

         gCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTa                                                                                                                                              <  1:331391/150‑1 (MQ=255)

           aaCTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATg                                                                                                                                            <  1:1684511/150‑1 (MQ=255)

             cTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGACCAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCg                                                                                                                                          >  2:958844/1‑150 (MQ=255)

                taCGAAACGTTGCTCACTCAACTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCATACGGAGAACATTCTCCG‑CTGCGAat                                                                                                                                       <  2:1530573/149‑1 (MQ=255)

                cACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAat                                                                                                                                       >  2:1138343/1‑150 (MQ=255)

                 aCGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCATGTAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAata                                                                                                                                      <  1:503763/150‑1 (MQ=255)

                  cGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAatat                                                                                                                                     >  1:1683053/1‑150 (MQ=255)

                   gAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCATACGGAGAACATTCTCCG‑CTGCGAattac                                                                                                                                    <  1:365082/150‑4 (MQ=38)

                   gAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATc                                                                                                                                    >  1:325499/1‑150 (MQ=255)

                          tGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTctct                                                                                                                             >  1:1549090/1‑150 (MQ=255)

                             tcactcacCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTAGACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCtgtg                                                                                                                          >  2:813856/1‑150 (MQ=255)

                             tcactcacCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCtgtg                                                                                                                          <  2:44929/150‑1 (MQ=255)

                                ctcacCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTg                                                                                                                       <  1:732787/150‑1 (MQ=255)

                                 tcacCTTGAGCATCTGAAAAAACACCAGGTAAATACCATTGATATTCGTTACACCACTCATGATCAATTCAAACAACAAACCGCAGATAGTCAGGCTGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTTCTgg                                                                                                                      <  1:921698/150‑1 (MQ=255)

                                  cacCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGc                                                                                                                     >  1:1230903/1‑150 (MQ=255)

                                  cacCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGc                                                                                                                     >  1:1117645/1‑150 (MQ=255)

                                   acCTTGAGCAGGTGCAAAAACACCAGGGAAATCCCATTGAAATTCGTTACACCACGCATTAACAATTCAAACTACATACCGCAGAAAGTCAGGCGGTAATTCGCT‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCg                                                                                                                    <  1:272690/150‑1 (MQ=255)

                                    ccTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGt                                                                                                                   <  2:1072889/150‑1 (MQ=255)

                                    ccTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGt                                                                                                                   >  1:802705/1‑150 (MQ=255)

                                      ttGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGa                                                                                                                 >  2:1430551/1‑150 (MQ=255)

                                           cagcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCGCGCATGAACAATTCAACCAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTc                                                                                                            >  2:1092069/1‑150 (MQ=255)

                                           cagcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTc                                                                                                            >  1:1566419/1‑150 (MQ=255)

                                            agcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCt                                                                                                           >  2:872590/1‑150 (MQ=255)

                                            agcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGATGGCGTGACGTTCt                                                                                                           >  2:878502/1‑150 (MQ=255)

                                              cTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGa                                                                                                         >  2:1573113/1‑150 (MQ=255)

                                                gCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCACTGTGCTGGCGTGCCGTTCTGAgg                                                                                                       >  2:476118/1‑150 (MQ=255)

                                                gCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTAAGGCGGTAATTCGCA‑TCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAgg                                                                                                       >  2:1370866/1‑150 (MQ=255)

                                                   aaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCg                                                                                                    >  1:1343432/1‑150 (MQ=255)

                                                   aaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCg                                                                                                    >  2:1302500/1‑150 (MQ=255)

                                                   aaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCg                                                                                                    <  2:497328/150‑1 (MQ=255)

                                                   aaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCg                                                                                                    <  1:1633265/150‑1 (MQ=255)

                                                   aaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCg                                                                                                    >  1:315035/1‑150 (MQ=255)

                                                      acacCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGG‑AATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCAt                                                                                                >  1:319286/1‑150 (MQ=255)

                                                          cAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGCGGCCGTCATGGa                                                                                             >  1:828270/1‑150 (MQ=255)

                                                          cAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGa                                                                                             >  1:1365869/1‑150 (MQ=255)

                                                            gggAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAg                                                                                           <  2:1477566/150‑1 (MQ=255)

                                                                aaTACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAACCCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCAtt                                                                                       <  2:262921/150‑1 (MQ=255)

                                                                aaTACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCAtt                                                                                       >  2:1707132/1‑150 (MQ=255)

                                                                  tACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTAc                                                                                     <  1:642144/150‑1 (MQ=255)

                                                                     cATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTc                                                                                  <  2:321180/150‑1 (MQ=255)

                                                                     cATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTc                                                                                  <  2:58116/150‑1 (MQ=255)

                                                                       ttGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAg                                                                                <  2:366308/150‑1 (MQ=255)

                                                                         gAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCt                                                                              <  1:1102360/150‑1 (MQ=255)

                                                                         gAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCt                                                                              <  1:1405152/150‑1 (MQ=255)

                                                                          aaaTTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCtt                                                                             <  1:690830/150‑1 (MQ=255)

                                                                            aTTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCTGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCTTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCATATTaa                                                                           >  1:1070959/1‑150 (MQ=255)

                                                                              tCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAg                                                                         >  2:371146/1‑150 (MQ=255)

                                                                               cGTTACACCACGCATGAACAATTCATACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAgg                                                                        <  2:776266/150‑1 (MQ=255)

                                                                               cGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAgg                                                                        >  1:854219/1‑150 (MQ=255)

                                                                                gTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCAGTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTCCTTCAGCTTAAAGGc                                                                       >  1:122211/1‑150 (MQ=255)

                                                                                gTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGATCTCCGTATGCGAATATCATTATCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGc                                                                       >  1:685515/1‑150 (MQ=255)

                                                                                  tACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCAt                                                                     <  1:1326163/150‑1 (MQ=255)

                                                                                  tACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGCCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCAt                                                                     <  1:1662149/150‑1 (MQ=255)

                                                                                   acacCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATc                                                                    <  2:516738/150‑1 (MQ=255)

                                                                                     accacGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGa                                                                  >  1:499940/1‑150 (MQ=255)

                                                                                     accacGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGa                                                                  >  1:506619/1‑150 (MQ=255)

                                                                                     accacGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGa                                                                  >  2:1335703/1‑150 (MQ=255)

                                                                                       cacGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATa                                                                >  2:650576/1‑150 (MQ=255)

                                                                                         cGCATGAACAATTCAAACCACAACCCGCAGAAAGTGAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGGATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGAATCGATaaa                                                              <  1:89470/150‑1 (MQ=255)

                                                                                           cATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGc                                                            >  2:1436030/1‑150 (MQ=255)

                                                                                           cATGAACAAGTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCGGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGc                                                            >  2:73525/1‑150 (MQ=255)

                                                                                              gAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCtt                                                         >  1:1432667/1‑150 (MQ=255)

                                                                                                aCAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTcc                                                       >  2:958618/1‑150 (MQ=255)

                                                                                                 cAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGATGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTccc                                                      >  2:440795/1‑150 (MQ=255)

                                                                                                 cAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTccc                                                      <  1:1218850/150‑1 (MQ=255)

                                                                                                    ttCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCggg                                                   <  2:1128307/150‑1 (MQ=255)

                                                                                                      cAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGATTATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCg                                                 >  2:1648093/1‑150 (MQ=255)

                                                                                                      cAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCg                                                 >  1:442929/1‑150 (MQ=255)

                                                                                                      cAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCg                                                 >  1:518359/1‑150 (MQ=255)

                                                                                                      cAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTACCCGGGCg                                                 >  2:259678/1‑150 (MQ=255)

                                                                                                       aaaCAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGt                                                >  2:666428/1‑150 (MQ=255)

                                                                                                          caacaaACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCGGTATGCGAATATCATTCTATGTGCTGGAGTGACGTTATTAGTCCGTCAATGAAGCATGACTTCAGAATAAAGTCATCGATAAAGCCTTCCCGGGCGTaaa                                             >  2:340206/1‑150 (MQ=255)

                                                                                                          caacaaACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTaaa                                             >  2:464923/1‑150 (MQ=255)

                                                                                                            acaaACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAg                                           >  1:636417/1‑150 (MQ=255)

                                                                                                               aaCCGCAGAAAGTCTGGCGGTATTTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGGTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGccc                                        <  1:869448/150‑1 (MQ=255)

                                                                                                                 ccGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGCCGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTATAAGCCctc                                      >  1:1056957/1‑150 (MQ=255)

                                                                                                                      gAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCggg                                 <  1:1362134/150‑1 (MQ=255)

                                                                                                                      gAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCggg                                 <  2:325499/150‑1 (MQ=255)

                                                                                                                      gAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCggg                                 <  1:651195/150‑1 (MQ=255)

                                                                                                                          gTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCa                             >  1:1559349/1‑150 (MQ=255)

                                                                                                                            cAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGc                           <  1:930741/150‑1 (MQ=255)

                                                                                                                              ggCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGt                         >  2:1638389/1‑150 (MQ=255)

                                                                                                                               gctgTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGtt                        <  2:1727235/147‑1 (MQ=255)

                                                                                                                                cGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTACAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTa                       <  2:1305955/150‑1 (MQ=255)

                                                                                                                                     aTTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGt                  >  1:1400109/1‑150 (MQ=255)

                                                                                                                                      ttCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTc                 <  1:1174229/150‑1 (MQ=255)

                                                                                                                                      ttCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGGGACGGTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTc                 <  2:78960/150‑1 (MQ=255)

                                                                                                                                       tCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGAATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGTGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCt                >  1:286229/1‑150 (MQ=255)

                                                                                                                                        cGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTa               <  1:79049/150‑1 (MQ=255)

                                                                                                                                        cGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTa               >  2:70218/1‑150 (MQ=255)

                                                                                                                                         gcatacGGAGAACATTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTa               >  1:1464893/6‑150 (MQ=255)

                                                                                                                                         gaatacGGAGAACATACTCCGTAAGCGAATATCATTCTCTGTGCTGGCGTGACCTTCTCAGGCCGTCATGGAAGCATTACTTCTGCTTAAAGGCATCGATAAAGCCTTCCCGGTCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTcga               >  2:1522110/6‑148 (MQ=255)

                                                                                                                                          cA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATc             <  1:1707132/150‑1 (MQ=255)

                                                                                                                                            tacGGAGAACATTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATcc            >  2:472161/3‑150 (MQ=255)

                                                                                                                                                 agaATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGcc       <  1:650576/150‑1 (MQ=255)

                                                                                                                                                      gTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTg  >  2:443108/1‑150 (MQ=255)

                                                                                                                                                      gTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAAGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTg  <  2:1567772/150‑1 (MQ=255)

CCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTG  >  NC_000913/3933575‑3933873

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: `‑`