BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,135,356	T→C	15.3%	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	15.3%	226.2 / 31.0	98	Q339R (CAG→CGG)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand): ref base T (42/40); new base C (7/8); total (49/49)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 106 total aligned reads displayed.

GAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGA  >  NC_000913/3135258‑3135454

gAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCCGAGATGCACAGCATAATGCGGCGCAcctg                                                                                                   >  2:331144/1‑100 (MQ=255)

  aGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                 >  1:330121/1‑100 (MQ=255)

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

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

     tttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCGCAGCATAATGCGGCGCAGCTGGCTgc                                                                                              >  2:2284197/1‑100 (MQ=255)

     tttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTgc                                                                                              >  1:2548337/1‑100 (MQ=255)

     tttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTgc                                                                                              >  1:54651/1‑100 (MQ=255)

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

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

      ttttCAGCAGGCTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTGCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCGgctggctg                                                                                               <  2:1229336/98‑1 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTgcgc                                                                                            <  2:595409/99‑1 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCt                                                                                           >  1:322116/1‑100 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCt                                                                                           >  1:1072187/1‑100 (MQ=255)

        ttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgcc                                                                                           <  1:2011675/100‑2 (MQ=255)

         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTg                                                                                          >  1:1716415/1‑100 (MQ=255)

         tCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTg                                                                                          <  1:1867587/100‑1 (MQ=255)

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

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

          cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccgc                                                                                         <  1:1568310/100‑4 (MQ=255)

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

          cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                          <  2:2340754/99‑3 (MQ=255)

           agcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccgca                                                                                        <  1:1691050/100‑5 (MQ=255)

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

               ggTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCAc                                                                                    <  2:1612375/100‑1 (MQ=255)

                gTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACg                                                                                   <  1:2125710/100‑1 (MQ=255)

                gTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACg                                                                                   >  2:1215656/1‑100 (MQ=255)

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

                  tCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCt                                                                                 <  2:1318779/100‑1 (MQ=255)

                      gTCTTCTCTACTGACGCCTGGCAGTTTCGCTGGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAAc                                                                             <  2:806173/100‑1 (MQ=255)

                       tcttctTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACg                                                                            >  2:679249/1‑100 (MQ=255)

                          tcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGc                                                                          >  2:1004253/1‑99 (MQ=255)

                           cttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTc                                                                        >  2:1842142/1‑100 (MQ=255)

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

                             ttACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGt                                                                      <  2:1818741/100‑1 (MQ=255)

                             ttACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGt                                                                      <  2:1818161/100‑1 (MQ=255)

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

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

                              tACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTa                                                                     <  1:931716/100‑1 (MQ=255)

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

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

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

                                    cGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAAGGGCTCGTAACTTTc                                                               >  2:1745586/1‑100 (MQ=255)

                                     gCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTcc                                                              <  1:2496033/100‑1 (MQ=255)

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

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

                                        tGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCat                                                            >  2:506524/1‑99 (MQ=255)

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

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

                                            aGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATAtt                                                         >  2:1646424/1‑98 (MQ=255)

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

                                             gTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGcc                                                      >  2:2272663/1‑100 (MQ=255)

                                             gTTTCGCTAGCTTCGCGGAGGTATCAGAGAGGCACAGAATAATGCGGCGCAGCGGGCGCCGCTGCCTCCCGCTAAACGGCTCGTAACTTTCCATATGGcc                                                      <  2:1397321/100‑1 (MQ=255)

                                               ttCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTg                                                    <  2:1515074/100‑1 (MQ=255)

                                               ttCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTg                                                    <  2:1474596/100‑1 (MQ=255)

                                                  gCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgca                                                 <  1:490071/100‑1 (MQ=255)

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

                                                     aGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgcagca                                              >  2:180319/1‑100 (MQ=255)

                                                      gCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGATTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgcagca                                              >  2:2152501/1‑99 (MQ=255)

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

                                                       cTTCGCGGAGGTATTAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATc                                            <  1:1966967/100‑1 (MQ=255)

                                                       cTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATc                                            <  2:657865/100‑1 (MQ=255)

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

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

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

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

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

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

                                                                ggTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcg                                   <  2:1707110/100‑1 (MQ=255)

                                                                ggTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcg                                   <  2:1402582/100‑1 (MQ=255)

                                                                 gTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgc                                  <  2:51697/100‑1 (MQ=255)

                                                                 gTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgc                                  >  1:2441700/1‑100 (MQ=255)

                                                                  tATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATCTAGCCTGGCAGCATCGtt                                         >  2:1587536/1‑92 (MQ=255)

                                                                  tATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcg                                 <  2:31472/100‑1 (MQ=255)

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

                                                                   ntCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCATAACTTTCCATATTGCCTGGCATCATCGTTGTaa                                     >  2:129/2‑95 (MQ=255)

                                                                       gagaTGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGCATCGTTTTAACGCGCGCaa                            >  1:1748205/1‑100 (MQ=255)

                                                                          aTGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAg                         <  2:2324551/100‑1 (MQ=255)

                                                                           tGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGc                        >  1:2135666/1‑100 (MQ=255)

                                                                             nacaGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCa                       >  2:1191291/2‑99 (MQ=255)

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

                                                                              acaGCATAATGCGGCGCAGCTGGCTGCGCTGTCTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATc                     <  1:1128962/100‑1 (MQ=255)

                                                                               caGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTCGCAGCATCGTTTTAACGCGCGCAATAGCAt                      >  2:1210882/1‑98 (MQ=255)

                                                                                aGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaa                   >  2:559498/1‑100 (MQ=255)

                                                                                aGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaa                   <  1:2246084/100‑1 (MQ=255)

                                                                                 gCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCaaa                  <  1:481741/100‑1 (MQ=255)

                                                                                 gCACAATTCTGCGCAGCTGGCTTCTCTGCCTCACTCTCAATGGCTCGTAACGTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCa                    <  1:2288525/98‑1 (MQ=255)

                                                                                  cgtgagtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGGGCAATAGCATCAAAg                 >  1:2520848/10‑100 (MQ=255)

                                                                                  cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAg                 <  2:1025974/100‑1 (MQ=255)

                                                                                  cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAg                 >  2:2300071/1‑100 (MQ=255)

                                                                                     gcgtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGGTTTAACGCGCGCAATAGCATCAAAGGt               >  1:1731614/7‑99 (MQ=255)

                                                                                      agtcagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTCTTAACGCGCGCAATAGCATCa                    >  2:1763524/6‑93 (MQ=255)

                                                                                       tGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTACCTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTAtt            >  2:304224/1‑100 (MQ=255)

                                                                                         cggcgCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTAttt           >  2:1054744/1‑99 (MQ=255)

                                                                                         cagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTg          >  1:157089/3‑100 (MQ=255)

                                                                                         cagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTACTTg          >  1:2214603/3‑100 (MQ=255)

                                                                                          agcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGTATCAAAGGTATTTGc         >  1:2525098/2‑100 (MQ=255)

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

                                                                                           gcgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGcc        <  1:1823451/100‑1 (MQ=255)

                                                                                                gctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGa     <  2:1043039/98‑1 (MQ=255)

                                                                                                 ctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGa  <  1:180321/100‑1 (MQ=255)

GAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGA  >  NC_000913/3135258‑3135454

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