BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,079,859	C→G	5.3%	D593H (GAT→CAT)	tktA ←	transketolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,079,859	0	C	G	5.3%	387.8 / 8.5	114	D593H (GAT→CAT)	tktA	transketolase 1
Reads supporting (aligned to +/- strand): ref base C (48/59); new base G (4/2); total (53/61)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 115 total aligned reads displayed.

TCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTC  >  NC_000913/3079711‑3080007

 ctGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCt                                                                                                                                                    <  1:225706/150‑1 (MQ=255)

      ggAGCTGATTCTCCGAAGTTGGTCATACCGACGATATCACCTTTCAGGCCAAGATACTTGTAGCAGTAGTCATCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGTTAAGCAGAATCCTGCttt                                                                                                                                               <  2:1336229/150‑2 (MQ=255)

      ggAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTg                                                                                                                                               >  2:1352697/1‑150 (MQ=255)

      ggAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTg                                                                                                                                               <  1:1554970/150‑1 (MQ=255)

      ggAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTg                                                                                                                                               <  2:431363/150‑1 (MQ=255)

        aGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTc                                                                                                                                             >  2:1046837/1‑150 (MQ=255)

         gCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGGATGCTGCttaccg                                                                                                                                            <  2:1435677/150‑5 (MQ=255)

           aGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCaaa                                                                                                                                          >  1:293631/1‑150 (MQ=255)

            gATTGACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAAt                                                                                                                                         <  2:1641416/150‑1 (MQ=255)

               tCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCg                                                                                                                                      <  2:1328143/150‑1 (MQ=255)

                 aCCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTc                                                                                                                                    <  1:1063728/150‑1 (MQ=255)

                 aCCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTc                                                                                                                                    >  1:293779/1‑150 (MQ=255)

                  ccGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCg                                                                                                                                   >  2:518694/1‑150 (MQ=255)

                               tACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATgg                                                                                                                      <  1:1454961/150‑1 (MQ=255)

                                aCCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGa                                                                                                                     >  2:1005773/1‑150 (MQ=255)

                                 ccGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGac                                                                                                                    >  1:1366014/1‑150 (MQ=255)

                                 ccGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGac                                                                                                                    >  2:416857/1‑150 (MQ=255)

                                   gacgaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTGTCGGCAGTACGGATTCCCGGTAAGCAGCATCCTGCTTGTCAAATGCGGCGGTAGACGGCATGGaccc                                                                                                                  >  2:1444927/1‑148 (MQ=255)

                                    acgaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcacc                                                                                                                 <  2:7565/150‑1 (MQ=255)

                                    acgaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcacc                                                                                                                 <  1:1072964/150‑1 (MQ=255)

                                      gaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcaccac                                                                                                               <  1:1352697/150‑1 (MQ=255)

                                          gCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                                                                                           >  2:1657081/1‑150 (MQ=255)

                                          gCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                                                                                           >  2:1290980/1‑150 (MQ=255)

                                          gCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                                                                                           >  1:880364/1‑150 (MQ=255)

                                          gCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                                                                                           >  2:1107271/1‑150 (MQ=255)

                                           cACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgcg                                                                                                          <  2:880364/150‑1 (MQ=255)

                                            aCCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACgcgcgc                                                                                                         >  2:336231/1‑150 (MQ=255)

                                                ttCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                                                     <  2:822095/150‑1 (MQ=255)

                                                ttCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                                                     <  1:744444/150‑1 (MQ=255)

                                                   aGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACg                                                                                                  <  2:1038445/150‑1 (MQ=255)

                                                   aGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACg                                                                                                  <  2:1468483/150‑1 (MQ=255)

                                                   aGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACg                                                                                                  <  2:293779/150‑1 (MQ=255)

                                                     gCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGcc                                                                                                <  2:9060/150‑1 (MQ=255)

                                                      ccAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATgg                                                                                                                      >  1:1501612/1‑127 (MQ=255)

                                                      ccAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATgg                                                                                                                      <  2:1501612/127‑1 (MQ=255)

                                                         aCATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCg                                                                                            <  2:1208616/150‑1 (MQ=255)

                                                          cATACTTGTACAAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATACTGCTTGTCAAATGCGTCGGTAGCCGGCATGGACACCCCGCGCGCTTTCACGCCTTCgg                                                                                           >  2:1550770/1‑150 (MQ=255)

                                                           aTACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                                                                                          <  2:487141/150‑1 (MQ=255)

                                                             aCTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcag                                                                                        <  1:597992/150‑1 (MQ=255)

                                                              cTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagt                                                                                       <  1:290703/150‑1 (MQ=255)

                                                                tGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtca                                                                                     >  1:1371436/1‑150 (MQ=255)

                                                                tGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                                                                                          >  1:1022775/1‑145 (MQ=255)

                                                                tGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                                                                                          <  2:1022775/145‑1 (MQ=255)

                                                                     cAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTTCGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtttt                                                                                >  2:1115669/1‑150 (MQ=255)

                                                                      agtagtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttttt                                                                               <  1:1534619/150‑1 (MQ=255)

                                                                      agtaatCATCAATAACCGAAACAACAGAAAAGCGTGCAGTAACCGCTTTCGGAAGTACCAATTAACAATAAACAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCCCGAGCGCTTTCACGCCTTCGCCAGTCAGttttt                                                                               >  2:209725/1‑150 (MQ=255)

                                                                        tagtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCg                                                                             <  1:252449/150‑1 (MQ=255)

                                                                        tagtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCg                                                                             <  2:1005493/150‑1 (MQ=255)

                                                                        tagtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCg                                                                             <  2:78925/150‑1 (MQ=255)

                                                                            cAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAgg                                                                         >  2:170139/1‑150 (MQ=255)

                                                                             aGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGAAGTCAGTTTTTCGTAGgc                                                                        >  1:1614995/1‑150 (MQ=255)

                                                                                   aCCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAAc                                                                  >  2:1120309/1‑150 (MQ=255)

                                                                                    cccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACa                                                                 >  1:299669/1‑150 (MQ=255)

                                                                                      cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                                                               >  1:433240/1‑150 (MQ=255)

                                                                                       gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGcc                                                              <  1:882112/150‑1 (MQ=255)

                                                                                       gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGcc                                                              <  2:694569/150‑1 (MQ=255)

                                                                                         ttCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAg                                                            <  1:806166/150‑1 (MQ=255)

                                                                                           cTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGtt                                                          >  1:794313/1‑150 (MQ=255)

                                                                                            tACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAAGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTc                                                         <  1:1108231/150‑1 (MQ=255)

                                                                                              cAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgc                                                                        <  2:276888/133‑1 (MQ=255)

                                                                                              cAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgc                                                                        >  1:276888/1‑133 (MQ=255)

                                                                                              cAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCaa                                                       >  1:1122196/1‑150 (MQ=255)

                                                                                              cAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCaa                                                       >  1:1615516/1‑150 (MQ=255)

                                                                                                gCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACt                                                     >  2:1612967/1‑150 (MQ=255)

                                                                                                gCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACt                                                     >  2:1249771/1‑150 (MQ=255)

                                                                                                 cAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACtt                                                    <  1:559255/150‑1 (MQ=255)

                                                                                                 cAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACtt                                                    <  1:518694/150‑1 (MQ=255)

                                                                                                  aaCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTc                                                   <  1:1005773/150‑1 (MQ=255)

                                                                                                    cgcgTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTg                                                 >  2:929021/1‑150 (MQ=255)

                                                                                                      cgTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGaa                                               >  2:467150/1‑150 (MQ=255)

                                                                                                      cgTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGaa                                               >  2:947991/1‑150 (MQ=255)

                                                                                                          cAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAACATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCgg                                           <  1:1063019/150‑1 (MQ=255)

                                                                                                            gTAACCGCTTTCGGCAGTACGCATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGGCGGTAGAACGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACCGACAGATCAACTACTGAACCGGTa                                         >  2:1008021/1‑150 (MQ=255)

                                                                                                              aaCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGc                                       >  2:949972/1‑150 (MQ=255)

                                                                                                              aaCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGc                                       <  1:548757/150‑1 (MQ=255)

                                                                                                               aCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCg                                      >  2:56912/1‑150 (MQ=255)

                                                                                                                        ggCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATc                             <  1:618017/150‑1 (MQ=255)

                                                                                                                        ggCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATc                             >  2:1001196/1‑150 (MQ=255)

                                                                                                                        ggCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATc                             <  1:1046837/150‑1 (MQ=255)

                                                                                                                        ggCAGTACGGATACAAGGTAAGCAGCATCCTGCTTTTCAAATGCGTCGGTAGACGGCATCGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGCTAGCGATGAAAATc                             >  2:782502/1‑150 (MQ=255)

                                                                                                                         gCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCa                            <  1:1107271/150‑1 (MQ=255)

                                                                                                                          cAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTGGTAGGCAGCAACAGCCAGTTCTACTTCTGAACCGGTAGCGATGaaattgtg                           <  2:1614995/150‑6 (MQ=255)

                                                                                                                            gTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGtt                         >  1:1658074/1‑150 (MQ=255)

                                                                                                                               cGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACCCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCg                      >  1:685797/1‑150 (MQ=255)

                                                                                                                                 gATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACTCCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGc                    <  1:949972/150‑1 (MQ=255)

                                                                                                                                  aTTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCt                   <  1:1120309/150‑1 (MQ=255)

                                                                                                                                     cACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGAc                >  2:1165609/1‑150 (MQ=255)

                                                                                                                                     cACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGAc                <  2:1371436/150‑1 (MQ=255)

                                                                                                                                       tgacaagaagGATGCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGAGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCg              >  2:1451710/5‑150 (MQ=255)

                                                                                                                                       tgacaaGCAGGATGCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCg              >  2:1126645/5‑150 (MQ=255)

                                                                                                                                        gacaaGCAGGATGCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCgg             >  2:588092/4‑150 (MQ=255)

                                                                                                                                         acaaGCAGGATGCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGgc            >  1:579396/3‑150 (MQ=255)

                                                                                                                                           aaGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGgcgc          <  1:1001196/150‑1 (MQ=255)

                                                                                                                                              cagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCGGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGt       <  1:882183/150‑1 (MQ=255)

                                                                                                                                               agcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTc      >  1:695498/1‑150 (MQ=255)

                                                                                                                                                gcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCt     >  1:406145/1‑150 (MQ=255)

                                                                                                                                                 caTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCACCTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGCCCGGCGCAGTCtt    >  2:1627779/1‑150 (MQ=255)

                                                                                                                                                 caTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCtt    >  1:525962/1‑150 (MQ=255)

                                                                                                                                                  aTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCttt   <  2:794313/150‑1 (MQ=255)

                                                                                                                                                   tCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTc  <  1:572193/150‑1 (MQ=255)

TCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTC  >  NC_000913/3079711‑3080007

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