BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,135,356	T→C	17.4%	Q339R (CAG→CGG) ‡	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,135,356	0	T	C	17.4%	432.5 / 83.8	196	Q339R (CAG→CGG)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand): ref base T (85/77); new base C (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 = 5.20e-01

Only 100 of 210 total aligned reads displayed.

CGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGA  >  NC_000913/3135257‑3135454

cGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgct                                                                                                    >  1:4233294/1‑100 (MQ=255)

cGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgct                                                                                                    >  1:4298368/1‑100 (MQ=255)

cGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgct                                                                                                    <  1:474960/100‑1 (MQ=255)

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

      tttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTgc                                                                                              <  1:76011/100‑1 (MQ=255)

      tttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCAGCGCAGCCAGCTgc                                                                                              <  2:3512419/100‑1 (MQ=255)

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

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

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

           cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCGCAGCATAAGGCGGCGCAGCTGGCTGCGCTGa                                                                                         >  1:2962403/1‑100 (MQ=255)

           cagcagGTTCTGGTCTTCCTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGCGGTATAAGAGATGCACAGCAAACTCCGGCGCAGATGGCTgagc                                                                                            <  1:647569/97‑1 (MQ=255)

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

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

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

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

               agGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCa                                                                                     >  2:2149116/1‑100 (MQ=255)

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

                  ttCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGc                                                                                  >  2:3305970/1‑100 (MQ=255)

                  ttCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCAGCGCAGCCAGCTgcgccgcattatgc                                                                                  <  2:3807155/100‑11 (MQ=37)

                     tGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATAAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTa                                                                               >  2:1060882/1‑100 (MQ=255)

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

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

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

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

                       gTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAAc                                                                             >  1:3202486/1‑100 (MQ=255)

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

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

                       gTCTTCTTTACTGACGCCTGGCAGGTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAAc                                                                             >  2:1299800/1‑100 (MQ=255)

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

                          ttcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGc                                                                          >  2:2074346/1‑100 (MQ=255)

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

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

                             tttACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCg                                                                       <  2:2391171/100‑1 (MQ=255)

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

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

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

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

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

                                 cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                   >  1:3183/1‑100 (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)

                                          ggCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCatat                                                          >  2:3632787/1‑100 (MQ=255)

                                           gCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATAtt                                                         <  2:4053494/100‑1 (MQ=255)

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

                                            cAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTg                                                        >  2:4140188/1‑100 (MQ=255)

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

                                             aGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGc                                                       <  1:4173335/100‑1 (MQ=255)

                                              gTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGcc                                                      <  2:4272160/100‑1 (MQ=255)

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

                                                   gCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgca                                                 >  2:3896595/1‑100 (MQ=255)

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

                                                       gCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAt                                             >  2:2119767/1‑100 (MQ=255)

                                                        cTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATc                                            >  1:2266332/1‑100 (MQ=255)

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

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

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

                                                            gcgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGttt                                        <  1:2923462/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)

                                                                   tATCAGAGATGCACAGCATAATGCGGCCCAGCTGGCTACGCTGACTCACGCTTCACGACCCGGAACTTTCCGTATCGCCTGGCAGCACCGTTTAAAcacg                                 >  1:412563/1‑100 (MQ=25)

                                                                    aTCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcgc                                >  2:1338023/1‑100 (MQ=255)

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

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

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

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

                                                                       agagATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCa                             <  1:2980646/100‑1 (MQ=255)

                                                                       agagATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTCTAACGCGCGCa                             >  2:1427588/1‑100 (MQ=255)

                                                                          gaTGCACATCATAATGCGGCGTAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACGTTCCATATTGGCTGGCAGCATCGTTTTAACGCGCGCAATa                          >  1:2851559/1‑100 (MQ=255)

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

                                                                             gCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCa                       >  1:4263246/1‑100 (MQ=255)

                                                                             gCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCa                       >  2:4294513/1‑100 (MQ=255)

                                                                             gCACAGCATAATGCGGCGCAGCTGGCTGCACTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCa                       <  2:1245876/100‑1 (MQ=255)

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

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

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

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

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

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

                                                                                   cgtgagtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAg                 >  2:2267654/10‑100 (MQ=255)

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

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

                                                                                     tgagtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGt               >  2:1397912/8‑100 (MQ=255)

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

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

                                                                                         gcggcgCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTAttt           >  1:719704/1‑100 (MQ=255)

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

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

                                                                                           agcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGc         >  2:1534817/2‑100 (MQ=255)

                                                                                           agcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGc         >  2:1482229/2‑100 (MQ=255)

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

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

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

CGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGA  >  NC_000913/3135257‑3135454

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