BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	158,868	C→T	A87T (GCC→ACC)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	158,868	0	C	T	75.9%	386.8 / 142.3	196	A87T (GCC→ACC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base C (25/22); new base T (73/75); total (99/97)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-02

Only 100 of 216 total aligned reads displayed.

ACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCT  >  NC_000913/158768‑158967

aCTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAggcgtggc                                                                                                     <  1:210983/101‑1 (MQ=255)

  ttCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAggcgtggtgt                                                                                                   <  2:210777/101‑1 (MQ=255)

     gATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAg                                                                                                >  2:211355/1‑101 (MQ=255)

      aTAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGt                                                                                               >  1:211356/1‑101 (MQ=255)

       tAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTg                                                                                              <  1:210932/101‑1 (MQ=255)

       tAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTg                                                                                              <  1:211274/101‑1 (MQ=255)

        aaTCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTgg                                                                                             >  2:211359/1‑101 (MQ=255)

        aaTCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTgg                                                                                             >  1:211360/1‑101 (MQ=255)

         aTCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGt                                                                                            >  1:211363/1‑101 (MQ=255)

           ctcCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTa                                                                                          >  2:211366/1‑101 (MQ=255)

           ctcCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTa                                                                                          <  1:211071/101‑1 (MQ=255)

            tcCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTAc                                                                                         >  1:211368/1‑101 (MQ=255)

              cGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCAACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGt                                                                                       <  1:211213/101‑1 (MQ=255)

              cGGGCCAAACATTACATGAGCCAGACGGAAACGGAGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTgg                                                                                             >  2:211369/1‑95 (MQ=255)

                   cAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCaaaa                                                                                  >  2:211372/1‑101 (MQ=255)

                     aaCATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATc                                                                                >  2:211375/1‑101 (MQ=255)

                       cATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCtt                                                                              <  2:210923/101‑1 (MQ=255)

                       cATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCtt                                                                              <  1:211035/101‑1 (MQ=255)

                       cATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGAGGCGTTAGTGGTTACGTCAAAATCtt                                                                              <  2:210969/101‑1 (MQ=255)

                          tACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCg                                                                           >  2:211377/1‑101 (MQ=255)

                          tACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCg                                                                           >  1:211378/1‑101 (MQ=255)

                            cATGAGCCAGACGGAAACGGAGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGc                                                                         >  2:211381/1‑101 (MQ=255)

                               gAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCttt                                                                      <  1:211024/101‑1 (MQ=255)

                               gAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCttt                                                                      <  1:211240/101‑1 (MQ=255)

                                 gCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCttttt                                                                    >  2:211384/1‑101 (MQ=255)

                                 gCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCttttt                                                                    <  2:211025/101‑1 (MQ=255)

                                  ccAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTg                                                                   <  2:210943/101‑1 (MQ=255)

                                  ccAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTg                                                                   >  1:211385/1‑101 (MQ=255)

                                  ccAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTg                                                                   <  2:211168/101‑1 (MQ=255)

                                   cAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGtt                                                                                           <  1:210945/76‑1 (MQ=255)

                                   cAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGGGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGc                                                                  >  2:211387/1‑101 (MQ=255)

                                     gACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCa                                                                >  1:211390/1‑101 (MQ=255)

                                     gACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCa                                                                <  2:211133/101‑1 (MQ=255)

                                       cGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAg                                                              <  2:210966/101‑1 (MQ=255)

                                        ggAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGt                                                             >  1:211392/1‑101 (MQ=255)

                                         gAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTa                                                            <  2:211122/101‑1 (MQ=255)

                                          aaaCGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTaa                                                           <  1:210821/101‑1 (MQ=255)

                                          aaaCGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTaa                                                           >  1:211395/1‑101 (MQ=255)

                                            aCGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACa                                                         <  1:210897/101‑1 (MQ=255)

                                             cGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAg                                                        >  1:211398/1‑101 (MQ=255)

                                             cGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAg                                                        <  1:210959/101‑1 (MQ=255)

                                             cGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAg                                                        <  1:211069/101‑1 (MQ=255)

                                              ggCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAgg                                                       >  2:211399/1‑101 (MQ=255)

                                              ggCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTC‑AAATCTTTCGGCttttt                                                                    >  1:211401/1‑87 (MQ=255)

                                               gCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGt                                                      <  2:211087/101‑1 (MQ=255)

                                                cGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTc                                                     <  1:210902/101‑1 (MQ=255)

                                                   cccACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTcgcg                                                  >  2:211404/1‑101 (MQ=255)

                                                     caccaGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCa                                                <  1:210990/101‑1 (MQ=255)

                                                      accaGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAc                                               >  1:211407/1‑101 (MQ=255)

                                                      accaGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAc                                               >  2:211408/1‑101 (MQ=255)

                                                         aGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAcgcc                                            >  1:211411/1‑101 (MQ=255)

                                                         aGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAcgcc                                            >  2:211412/1‑101 (MQ=255)

                                                           gcggcAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAcgccgc                                          <  2:211181/101‑1 (MQ=255)

                                                            cggcAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAcgccgcc                                         >  1:211415/1‑101 (MQ=255)

                                                             ggcAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAcgccgccg                                        >  2:211417/1‑101 (MQ=255)

                                                             ggcAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCAcgccgccg                                        <  2:211067/101‑1 (MQ=255)

                                                               cAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACgccgccgcc                                      >  1:211420/1‑101 (MQ=255)

                                                               cAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACgccgccgcc                                      >  2:211419/1‑101 (MQ=255)

                                                               cAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAAAAGGTCGCGCACgccgccgcc                                      <  2:211215/101‑1 (MQ=255)

                                                                aGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCa                                     <  1:211010/101‑1 (MQ=255)

                                                                  ttACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAc                                   <  2:211155/101‑1 (MQ=255)

                                                                   tACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACgccgccgcc                                      <  2:210835/97‑1 (MQ=255)

                                                                   tACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAcc                                  >  1:211423/1‑101 (MQ=255)

                                                                   tACGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAcc                                  <  1:211148/101‑1 (MQ=255)

                                                                    aCGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAcca                                 <  1:211003/101‑1 (MQ=255)

                                                                    aCGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAcca                                 <  1:211051/101‑1 (MQ=255)

                                                                    aCGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAcca                                 >  2:211425/1‑101 (MQ=255)

                                                                     cGGAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTcg                                                    >  2:211430/1‑81 (MQ=255)

                                                                     cGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAccag                                >  1:211428/1‑101 (MQ=255)

                                                                       gAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAccagcc                              >  2:211434/1‑101 (MQ=255)

                                                                       gAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAccagcc                              <  2:210991/101‑1 (MQ=255)

                                                                       gAACAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAccagcc                              >  1:211432/1‑101 (MQ=255)

                                                                        aaCAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAccagcca                             <  2:211144/101‑1 (MQ=255)

                                                                        aaCAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAccagcca                             <  1:210993/101‑1 (MQ=255)

                                                                        aaCAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAAccagcca                             <  2:211392/101‑1 (MQ=255)

                                                                          cAGTTTGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAgg                           >  1:211437/1‑101 (MQ=255)

                                                                          cAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAgg                           <  1:211085/101‑1 (MQ=255)

                                                                           aGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGc                          <  2:211054/101‑1 (MQ=255)

                                                                             tttGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCtt                        >  2:211438/1‑101 (MQ=255)

                                                                             tttGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCtt                        <  1:211107/101‑1 (MQ=255)

                                                                              ttGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTc                       <  1:211119/101‑1 (MQ=255)

                                                                               tGCGCACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCg                      <  2:211050/101‑1 (MQ=255)

                                                                                gcgcACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGt                     >  1:211441/1‑101 (MQ=255)

                                                                                  gcACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTAt                   <  2:211000/101‑1 (MQ=255)

                                                                                   cACCTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATc                  >  1:211443/1‑101 (MQ=255)

                                                                                     ccTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATccc                <  2:211038/101‑1 (MQ=255)

                                                                                     ccTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATccc                >  1:211445/1‑101 (MQ=255)

                                                                                      cTGCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCg               <  1:211302/101‑1 (MQ=255)

                                                                                        gCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCt             >  1:211447/1‑101 (MQ=255)

                                                                                        gCTCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCt             >  2:211448/1‑101 (MQ=255)

                                                                                          tCAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCttt           <  1:211152/101‑1 (MQ=255)

                                                                                           cAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCttt           >  1:211452/1‑100 (MQ=255)

                                                                                           cAGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCt             <  2:211279/98‑1 (MQ=255)

                                                                                            aGGCGTGGTGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTAt         >  2:211454/1‑101 (MQ=255)

                                                                                             ggcgtggcgtTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTAtt        >  1:211455/1‑101 (MQ=255)

                                                                                              gcgtggcgtTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTg       <  1:211130/101‑1 (MQ=255)

                                                                                                gtggtgtTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAg     <  2:211077/101‑1 (MQ=255)

                                                                                                gtggtgtTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAg     <  2:210977/101‑1 (MQ=255)

                                                                                                 tggtgtTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGc    <  1:210967/101‑1 (MQ=255)

                                                                                                  ggcgtTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGcc   >  1:211459/1‑101 (MQ=255)

ACTTCGATAATCTCCGGGCCAAACATTACATGAGCCAGACGGAAACGGCGACCCACCAGGCGGCAGTTACGGAACAGTTTGCGCACCTGCTCAGGCGTGGCGTTAGTGGTTACGTCAAAATCTTTCGGCTTTTTGCCAAGTAACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCT  >  NC_000913/158768‑158967

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`