BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,933,723	T→C	6.5%	C125R (TGT→CGT) ‡	rbsD →	D‑ribose pyranase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,933,723	0	T	C	6.5%	362.4 / 11.7	123	C125R (TGT→CGT) ‡	rbsD	D‑ribose pyranase
Reads supporting (aligned to +/- strand): ref base T (62/53); new base C (4/4); total (66/57)
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.21e-01

Only 100 of 126 total aligned reads displayed.

ACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCC  >  NC_000913/3933574‑3933868

actaTAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAAt                                                                                                                                                   <  2:1174564/147‑1 (MQ=255)

acCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGATAGTCAGGCGGTAATTCGCA‑GCGGAGAAt                                                                                                                                                   <  1:1136993/150‑1 (MQ=255)

acCATAATCCGCAACTCAACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAAt                                                                                                                                                   >  1:172985/1‑150 (MQ=255)

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

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

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

          gCAACTCCACGAAACGGTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCATACGGAGAACATTCTCCgc                                                                                                                                          <  1:320911/150‑2 (MQ=255)

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

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

              cTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCg                                                                                                                                     >  1:224170/1‑150 (MQ=255)

              cTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGCTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCg                                                                                                                                     >  2:571640/1‑150 (MQ=255)

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

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

                 cACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAat                                                                                                                                  <  1:485860/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)

                     aaaCGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCa                                                                                                                              <  1:1233352/150‑1 (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:1117645/1‑150 (MQ=255)

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

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

                                        tGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGAc                                                                                                           <  1:1500541/150‑1 (MQ=255)

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

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

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

                                              gcTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTg                                                                                                     <  2:1625418/150‑1 (MQ=255)

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

                                                 gCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAgg                                                                                                  >  1:1394022/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:497328/150‑1 (MQ=255)

                                                    aaaaCACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCg                                                                                               >  2:1302500/1‑150 (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)

                                                            aGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGaa                                                                                       >  2:837068/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)

                                                                     ccATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACtt                                                                              <  2:728311/150‑1 (MQ=255)

                                                                      cATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTc                                                                             <  2:321180/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)

                                                                             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)

                                                                                             aTGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGcc                                                      <  1:965521/150‑1 (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:518359/1‑150 (MQ=255)

                                                                                                       cAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCg                                            >  1:442929/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)

                                                                                                           caacaaACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTACTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTaaa                                        >  2:237984/1‑150 (MQ=255)

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

                                                                                                                aaCCGCAGATAGTCAGGCGGTAATTCGCA‑GCTGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGGGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGccc                                   <  1:863891/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                             cAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGc                      >  1:685092/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:581913/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                                           cA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATc        >  1:784699/1‑150 (MQ=255)

                                                                                                                                            atacGGAGAACATTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATc        >  1:1648340/4‑150 (MQ=255)

                                                                                                                                            a‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATcc       <  2:467303/150‑1 (MQ=255)

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

ACCATAATCCGCAACTCCACGAAACGTTGCTCACTCACCTTGAGCAGCTGCAAAAACACCAGGGAAATACCATTGAAATTCGTTACACCACGCATGAACAATTCAAACAACAAACCGCAGAAAGTCAGGCGGTAATTCGCA‑GCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCC  >  NC_000913/3933574‑3933868

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