BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,693,826	A→C	5.4%	R146R (CGT→CGG)	uidB ←	glucuronide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,693,826	0	A	C	5.4%	245.1 / 2.2	103	R146R (CGT→CGG)	uidB	glucuronide transporter
Reads supporting (aligned to +/- strand): ref base A (45/50); new base C (4/3); total (49/54)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 104 total aligned reads displayed.

AGCACCATTCCGGCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGA  >  NC_000913/1693677‑1693974

aGCACCATTCCGGCAATCGCCAGCACATTTGTCCAGAAATGGTATACCGACACCATCGCATCAGGGCAGGAGTTCTTTATGCTAGGTCCTATCAGAAATGCCAGGCAGCCAAAGGTCAATGAAGCGGCATTCCCACGAGCCCCGCCCAGa                                                                                                                                                      >  1:579717/1‑150 (MQ=255)

 gcagcaTTCAGGCAATCTCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGGTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGAc                                                                                                                                                     <  2:170299/150‑1 (MQ=255)

  cACCATTCCGGCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACg                                                                                                                                                    >  2:1820718/1‑150 (MQ=255)

        tCCGGCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGcgcg                                                                                                                                              >  2:1561361/1‑150 (MQ=255)

         ccGGCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCGTCTGGGCGCgg                                                                                                                                             <  1:1010793/150‑1 (MQ=255)

         ccGGCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCgg                                                                                                                                             <  1:1029368/150‑1 (MQ=255)

          cGGCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCGTCTGGGCGCggc                                                                                                                                            <  1:1669516/150‑2 (MQ=255)

            gCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGAt                                                                                                                                          <  2:1714689/150‑1 (MQ=255)

                tCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATtgtg                                                                                                                                      <  2:935341/150‑1 (MQ=255)

                 cGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTgg                                                                                                                                     <  1:1234195/150‑1 (MQ=255)

                  gCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCATAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGgt                                                                                                                                    <  1:813225/150‑1 (MQ=255)

                   ccAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGgtt                                                                                                                                   >  2:434857/1‑150 (MQ=255)

                    cAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGgttg                                                                                                                                  <  1:926151/150‑1 (MQ=255)

                    cAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCACCGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGCAGCGGCACTCCCACGAGCCGCGCCCCGCCGGGCGCGGGATTGTGgtgg                                                                                                                                  >  2:884439/1‑148 (MQ=255)

                    cAGCACAATTGACCAGAAATGGTATACCGACACCATCTCTACCGGGCTGGACTTCTTAATGCTCGGTCCTAGCAGAAATCCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGgttg                                                                                                                                  <  2:1110453/150‑1 (MQ=255)

                     aGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCGCGGGCCTCTCAGAACTGCCAGGCAGACAAAGGTCAATGAAGCGGCAAGCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGgttgt                                                                                                                                 <  1:272035/150‑1 (MQ=255)

                      ccacaATTGTCCACAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTTTCAGAAATGCCAGGCACACAAAGGTCAATGAAGCGGCAATCCCCCGAGCCGCGCCCAGACGGGCGCGGGATTGTGgttgtt                                                                                                                                <  2:333369/149‑1 (MQ=255)

                         caATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTggg                                                                                                                             >  1:175070/1‑150 (MQ=255)

                                ccAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATcgc                                                                                                                      >  2:1018558/1‑150 (MQ=255)

                                ccAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATcgc                                                                                                                      >  1:1796441/1‑150 (MQ=255)

                                 cAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATcgcg                                                                                                                     >  1:672257/1‑150 (MQ=255)

                                 cAGAAATGGTATACAGACACCATCTCTTCAGGGGTGGAGTTCTTAATGCTCGGTCCTTTCAGAAATGCCAGGCAGACAAAGGTCGAAGAAGCGGCAATCCCACTATCCTCGCCAAGACGGTCGCACGATTTTGTTTGTTGGGTCCTcgcg                                                                                                                     >  2:103061/1‑150 (MQ=255)

                                  aGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCgg                                                                                                                    >  1:1445428/1‑150 (MQ=255)

                                    acaTGGAAAACCGACCCCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGGCAATGAAGCGGCAATCCCACGCGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGGCATCGCGGTa                                                                                                                  <  2:911241/148‑1 (MQ=255)

                                    aaaTGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTa                                                                                                                  <  1:201686/150‑1 (MQ=255)

                                     aaTGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAg                                                                                                                 <  1:13554/150‑1 (MQ=255)

                                      aTGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGc                                                                                                                <  2:781903/150‑1 (MQ=255)

                                           ataCCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTg                                                                                                           <  2:1042712/150‑1 (MQ=255)

                                             aCCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGACTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGaa                                                                                                         >  1:136551/1‑150 (MQ=255)

                                                    ccATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGCCGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAg                                                                                                  >  2:508658/1‑150 (MQ=255)

                                                    ccATCTCTTCCGGGCTGGAGATCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGCCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAg                                                                                                  <  1:1458094/150‑1 (MQ=255)

                                                     cATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAgg                                                                                                 <  1:378380/150‑1 (MQ=255)

                                                     cATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAgg                                                                                                 >  1:474377/1‑150 (MQ=255)

                                                      aTCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGa                                                                                                <  1:1031209/150‑1 (MQ=255)

                                                         tctTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGAGCGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAata                                                                                             <  1:370113/150‑1 (MQ=255)

                                                         tctTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAata                                                                                             <  1:370107/150‑1 (MQ=255)

                                                         tctTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAata                                                                                             >  1:133764/1‑150 (MQ=255)

                                                          ctTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAatat                                                                                            >  2:579241/1‑150 (MQ=255)

                                                          ctTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAatat                                                                                            >  2:591017/1‑150 (MQ=255)

                                                          ctTCCGGGCTGGAGTTCTTAATGATCGGTCCTAACAGAAATGCCAGGAAGACAAAGGTCAATGAAGCGGCAATCACAAGAGCCGCGCCCAGACGGGCGCTGGATTGTGGTTGTTTGGTCATCGCGGAAGCTAGTGAACCATAAGGAatat                                                                                            >  1:1738946/1‑150 (MQ=255)

                                                                 gCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCAcccg                                                                                     >  1:1379277/1‑148 (MQ=255)

                                                                   tGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGc                                                                                   <  1:1094010/150‑1 (MQ=255)

                                                                       gTTCTTACTGCTCGGGCCTATCAGAAAGGCCAGGCACACAAAGGTCAATGAAGCGGCAATCCCCCGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTa                                                                               <  2:1173801/150‑1 (MQ=255)

                                                                       gTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTa                                                                               >  1:1834809/1‑150 (MQ=255)

                                                                         tCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGCCGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGc                                                                             >  2:1378903/1‑150 (MQ=255)

                                                                         tCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGc                                                                             >  1:1732888/1‑150 (MQ=255)

                                                                         tCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGCTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGc                                                                             >  1:147532/1‑150 (MQ=255)

                                                                          cTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCa                                                                            >  2:1327538/1‑150 (MQ=255)

                                                                          cTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCa                                                                            >  1:919379/1‑150 (MQ=255)

                                                                          cTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCa                                                                            >  1:771526/1‑150 (MQ=255)

                                                                             aaTGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAg                                                                         <  1:1648114/150‑1 (MQ=255)

                                                                             aaTGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAg                                                                         <  1:1491784/150‑1 (MQ=255)

                                                                               tGCTCGGTCCTATCAGAAATGCCAGGCAGCCAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCTTAGATGGAGCATGTGCTCCATTAGGCATATTCACCATGATTTAGAAAAGcc                                                                       >  2:294331/1‑150 (MQ=255)

                                                                               tGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGcc                                                                       >  1:109960/1‑150 (MQ=255)

                                                                                gCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGATTATTCACCCGGCTGTAGCAACGccc                                                                      >  1:620797/1‑150 (MQ=255)

                                                                                 cTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGCCGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACGAGGCTGTAGCAAAGCCCg                                                                     >  2:142342/1‑150 (MQ=255)

                                                                                  tCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGa                                                                    <  1:530300/150‑1 (MQ=255)

                                                                                     gTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGc                                                                 >  2:1057725/1‑150 (MQ=255)

                                                                                        cTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCa                                                              <  2:1700372/150‑1 (MQ=255)

                                                                                          aTCAGAAATGCAAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCACAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCAtt                                                            >  2:1045158/1‑149 (MQ=255)

                                                                                           tCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGt                                                           >  2:655130/1‑150 (MQ=255)

                                                                                           tCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGGGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAACGCCCGAGGCCCATGt                                                           >  2:1156710/1‑150 (MQ=255)

                                                                                            cAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGGCCTCGCGGTAGCAAGTGAACCATAAGGAATATTAACCAGGCTGTAGCAAAGCCCGTGGCCCATGTa                                                          >  1:1644789/1‑150 (MQ=255)

                                                                                             aGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAg                                                         >  1:70760/1‑150 (MQ=255)

                                                                                                  tGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCaa                                                    >  1:1257009/1‑150 (MQ=255)

                                                                                                   gCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCaaa                                                   <  1:675423/150‑1 (MQ=255)

                                                                                                   gCCAGGCAGACAAAGGGCCATGACGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGGCaaa                                                   <  2:1176858/150‑1 (MQ=255)

                                                                                                    ccAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAat                                                  >  2:1060962/1‑150 (MQ=255)

                                                                                                    ccAGGCAGACAAAGCTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGGGGTTGTTGGGTCAACGCGGTAGCAAGTGAACCATAACGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCAAGTAGGTCaaaa                                                  >  2:1714332/1‑149 (MQ=255)

                                                                                                       ggCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATg                                               >  1:1021236/1‑150 (MQ=255)

                                                                                                       ggCAGACAAAGGGCAATGAAGCGGCTATCCCACGAGCCGCGGCCAGTCGGGCGCGGGATTGTGGTTGTTGGGTAATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATg                                               <  2:648181/150‑1 (MQ=255)

                                                                                                          aGACATAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCCGACGGGCGCGGGTATGTGGGTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATCTTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGGCAAATATGCAt                                            <  1:103258/150‑1 (MQ=255)

                                                                                                           gACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCATGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATa                                           >  2:913853/1‑150 (MQ=255)

                                                                                                           gACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATa                                           <  2:895811/150‑1 (MQ=255)

                                                                                                           gACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATa                                           <  2:373110/150‑1 (MQ=255)

                                                                                                            aCAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATac                                          <  2:1796441/150‑1 (MQ=255)

                                                                                                             cAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATaca                                         >  2:725696/1‑150 (MQ=255)

                                                                                                                           gCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCa                           >  1:1562402/1‑150 (MQ=255)

                                                                                                                            cGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGGTGGGTCATCGCGGTAGCAAGTGAGCGATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGGATACACCACTTTGCTAccct                          >  2:1653739/1‑148 (MQ=255)

                                                                                                                             ggCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATg                         >  2:1223362/1‑150 (MQ=255)

                                                                                                                                aaTCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTTAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCt                      <  2:109960/150‑1 (MQ=255)

                                                                                                                                  tCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTcc                    <  2:771526/150‑1 (MQ=255)

                                                                                                                                   cccACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCa                   >  2:805244/1‑150 (MQ=255)

                                                                                                                                    ccACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAg                  >  2:573431/1‑150 (MQ=255)

                                                                                                                                      aCGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTc                >  2:848372/1‑150 (MQ=255)

                                                                                                                                       cGAGCCGCGCCCAGACGGGCGCGGGAGTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCg               <  1:797177/150‑1 (MQ=255)

                                                                                                                                        gAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCgg              <  2:54774/150‑1 (MQ=255)

                                                                                                                                        gAGCCGCGCCCAGACGGGCGCGGGATTGGGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCGCCAGTCgg              <  1:1644939/150‑1 (MQ=255)

                                                                                                                                         cgCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGt             <  1:1819781/149‑1 (MQ=255)

                                                                                                                                          gCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTc            <  1:1370929/150‑1 (MQ=255)

                                                                                                                                          gCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGAAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTc            <  1:1362271/150‑1 (MQ=255)

                                                                                                                                           ccGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCa           <  1:1820718/150‑1 (MQ=255)

                                                                                                                                            cgcgCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAg          <  1:591017/150‑1 (MQ=255)

                                                                                                                                            cgcgCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAg          >  2:1649116/1‑150 (MQ=255)

                                                                                                                                            cgcgCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAg          <  1:1596511/150‑1 (MQ=255)

                                                                                                                                            cgcgCCCAGACGGGCGCGGGATGGTGGTTGGTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGGCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAg          <  2:1237729/150‑1 (MQ=255)

                                                                                                                                               gcCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCAc       <  2:133764/150‑1 (MQ=255)

                                                                                                                                                   aGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAg   <  1:568501/150‑1 (MQ=255)

                                                                                                                                                   aGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAg   <  1:936541/150‑1 (MQ=255)

                                                                                                                                                    gACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGa  <  1:434857/150‑1 (MQ=255)

AGCACCATTCCGGCAATCGCCAGCACAATTGTCCAGAAATGGTATACCGACACCATCTCTTCCGGGCTGGAGTTCTTAATGCTCGGTCCTATCAGAAATGCCAGGCAGACAAAGGTCAATGAAGCGGCAATCCCACGAGCCGCGCCCAGACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGA  >  NC_000913/1693677‑1693974

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