BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,135,357	G→A	17.3%	Q339* (CAG→TAG) ‡	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,135,357	0	G	A	17.3%	459.1 / 78.9	196	Q339* (CAG→TAG)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand): ref base G (85/77); new base A (18/16); total (103/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 208 total aligned reads displayed.

GAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACA  >  NC_000913/3135258‑3135456

  aGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                   <  1:4084674/100‑1 (MQ=255)

  aGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                   <  2:2455593/100‑1 (MQ=255)

      ttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTgcg                                                                                               >  2:2432900/1‑100 (MQ=255)

      ttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTgcg                                                                                               <  1:2578058/100‑1 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCt                                                                                             <  1:3482005/100‑1 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCt                                                                                             <  2:4291940/100‑1 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCt                                                                                             <  1:57275/100‑1 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCt                                                                                             <  1:2842039/100‑1 (MQ=255)

         tCAGCAGGTTCTTGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                            <  2:3954165/100‑3 (MQ=255)

         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                            <  2:3399265/100‑3 (MQ=255)

         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                            <  2:1637540/100‑3 (MQ=255)

         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                            <  2:2121748/100‑3 (MQ=255)

         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                            <  1:883849/100‑3 (MQ=255)

         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgcc                                                                                             <  1:2180942/99‑2 (MQ=255)

          cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccgc                                                                                           <  2:2094390/100‑4 (MQ=255)

          cagcagGCTCTTGTCTTCTTTACTGAAGCCTGGCAGTTGGGCTAGCTTCGCGGAGGTAGCAGAGATACAAAACATAATGCGGCGCAgctggctgtg                                                                                               >  1:1622034/1‑94 (MQ=255)

           agcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAc                                                                                          <  1:292529/100‑1 (MQ=255)

            gcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACt                                                                                         >  2:556650/1‑100 (MQ=255)

            gcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccgcat                                                                                         <  2:872011/100‑6 (MQ=255)

               ggTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCAc                                                                                      >  2:2274228/1‑100 (MQ=255)

               ggTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCAc                                                                                      >  2:1994848/1‑100 (MQ=255)

                 ttCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGc                                                                                    <  1:962181/100‑1 (MQ=255)

                 ttCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGc                                                                                    <  2:1509347/100‑1 (MQ=255)

                    tGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTa                                                                                 <  1:4129050/100‑1 (MQ=255)

                     ggTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTaa                                                                                >  1:2678750/1‑100 (MQ=255)

                     ggTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTaa                                                                                >  2:120562/1‑100 (MQ=255)

                     ggTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTaa                                                                                >  2:1406514/1‑100 (MQ=255)

                     ggTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTa                                                                                 >  1:1680329/1‑99 (MQ=255)

                      gTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAAc                                                                               <  2:1444345/100‑1 (MQ=255)

                      gTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAAc                                                                               <  2:4218084/100‑1 (MQ=255)

                      gTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCCTAAc                                                                               <  1:234766/100‑1 (MQ=255)

                       tcttctTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACg                                                                              <  2:4332600/100‑1 (MQ=255)

                          tcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCt                                                                           >  1:3071470/1‑100 (MQ=255)

                           cttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGGGACGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACCCACGCTTAACGGCTc                                                                          >  1:176995/1‑100 (MQ=255)

                              tACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTa                                                                       >  2:1319937/1‑100 (MQ=255)

                              tACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTa                                                                       >  2:3819055/1‑100 (MQ=255)

                               aCTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTaa                                                                      <  1:3621612/100‑1 (MQ=255)

                               aCTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTaa                                                                      <  1:670108/100‑1 (MQ=255)

                                cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                     >  1:4272129/1‑100 (MQ=255)

                                cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                     >  1:3183/1‑100 (MQ=255)

                                cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                     <  1:3202586/100‑1 (MQ=255)

                                cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                     >  2:93595/1‑100 (MQ=255)

                                cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                     >  2:624971/1‑100 (MQ=255)

                                cTGACGCCTGGCAGTTTCGCTACCTTCGCGGTGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                     <  1:2965487/100‑1 (MQ=255)

                                  gACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACtt                                                                   <  1:1135981/100‑1 (MQ=255)

                                    cGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTc                                                                 >  1:2754977/1‑100 (MQ=255)

                                      ccTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTccg                                                               >  2:563141/1‑99 (MQ=255)

                                      ccTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCa                                                               <  1:3150076/100‑1 (MQ=255)

                                      ccTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCa                                                               >  1:299641/1‑100 (MQ=255)

                                       cTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCat                                                              <  1:2744944/100‑1 (MQ=255)

                                       cTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCa                                                               >  1:2731964/1‑99 (MQ=255)

                                         ggCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCatat                                                            <  1:1024665/100‑1 (MQ=255)

                                           cAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATAtt                                                           >  2:2915598/1‑99 (MQ=255)

                                           cAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTg                                                          <  2:3627563/100‑1 (MQ=255)

                                           cAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTg                                                          <  1:2400821/100‑1 (MQ=255)

                                                 cGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgc                                                    >  2:905436/1‑100 (MQ=255)

                                                   cTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgcag                                                  <  2:181619/100‑1 (MQ=255)

                                                    tAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgcagc                                                 <  2:34363/100‑1 (MQ=255)

                                                       cTTCGCGGAGGTATCAGAGATGCACCGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATc                                              >  1:2864281/1‑100 (MQ=255)

                                                        ttCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCg                                             <  1:1873628/100‑1 (MQ=255)

                                                         tCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGt                                            >  1:3928865/1‑100 (MQ=255)

                                                          cgcgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGtt                                           <  2:1308162/100‑1 (MQ=255)

                                                            cgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGtttt                                         <  2:3342198/100‑1 (MQ=255)

                                                             ggAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTa                                        <  2:4233324/100‑1 (MQ=255)

                                                              gAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTaa                                       <  1:3466507/100‑1 (MQ=255)

                                                              gAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTaa                                       <  2:2747232/100‑1 (MQ=255)

                                                                  tATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcg                                   >  2:4300031/1‑100 (MQ=255)

                                                                     cAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACgcgcgc                                <  2:1707077/100‑1 (MQ=255)

                                                                     cAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACgcgcgc                                >  1:4190953/1‑100 (MQ=255)

                                                                     cAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACgcgcgc                                <  2:1359383/100‑1 (MQ=255)

                                                                      agagATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCa                               >  2:711572/1‑100 (MQ=255)

                                                                      agagATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCa                               >  2:1766988/1‑100 (MQ=255)

                                                                         gaTGCACAGCATAATGCGGCGCCGCTGGCTGCGCTGACTCACGCTTAACGGCTCGGAACTTTCCATATTTCCTAGCAGCATCGTTTTAACGCGCGCaa                              >  1:1971209/1‑98 (MQ=255)

                                                                             cacaGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCAt                        >  1:1303918/1‑100 (MQ=255)

                                                                             cacaGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCAt                        >  2:3067780/1‑100 (MQ=255)

                                                                               caGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCa                      <  1:3407149/100‑1 (MQ=255)

                                                                                 gCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaaa                    >  1:1116034/1‑100 (MQ=255)

                                                                                 gCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaaa                    >  2:2493501/1‑100 (MQ=255)

                                                                                 gCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaaa                    >  1:2281504/1‑100 (MQ=255)

                                                                                 gCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaaa                    >  1:2912234/1‑100 (MQ=255)

                                                                                   aTAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAgg                  >  2:260903/1‑100 (MQ=255)

                                                                                   aTAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAgg                  <  2:1555965/100‑1 (MQ=255)

                                                                                   aTAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAgg                  <  2:4298399/100‑1 (MQ=255)

                                                                                   aTAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCCCAATAGCATCAAAgg                  >  1:992694/1‑100 (MQ=255)

                                                                                    tAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGt                 <  2:2723059/100‑1 (MQ=255)

                                                                                      agtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTAt               >  1:2176378/6‑100 (MQ=255)

                                                                                       gtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCACAGGTAtt              >  1:4239554/5‑100 (MQ=255)

                                                                                       ggcagcgcAGCCAGCTGCGCTGACTCAAGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAACGTAtt              >  1:1382563/5‑100 (MQ=255)

                                                                                        tcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTAttt             >  2:3799493/4‑100 (MQ=255)

                                                                                         cagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTg            >  2:1263571/3‑100 (MQ=255)

                                                                                          agcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATCTTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGc           >  1:3799939/2‑100 (MQ=255)

                                                                                           gcgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGcc          >  2:2203293/1‑100 (MQ=255)

                                                                                             gcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCgg        <  2:1116038/100‑1 (MQ=255)

                                                                                             gcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCgg        <  2:547309/100‑1 (MQ=255)

                                                                                             gcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCgg        <  1:736376/100‑1 (MQ=255)

                                                                                             gcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCgg        <  2:335532/100‑1 (MQ=255)

                                                                                                gctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATg     <  1:3349090/100‑1 (MQ=255)

                                                                                                gctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATg     <  2:330950/100‑1 (MQ=255)

                                                                                                 ctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGa    >  2:4123310/1‑100 (MQ=255)

                                                                                                   ggctgCGCTGACTCACGCTTAACGGCTCGTACATTTCCATATTGCCTGGCAGCATCGTTTTAACGCCCGCAATAGCATCAAAGGTATTTGCCGGAGGACa  >  2:4178348/1‑100 (MQ=255)

GAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACA  >  NC_000913/3135258‑3135456

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