BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,338,553	A→G	H96R (CAC→CGC)	yciS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,338,553	0	A	G	98.2%	518.0 / ‑1.1	166	H96R (CAC→CGC)	yciS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (65/98); total (68/98)
Fisher's exact test for biased strand distribution p-value = 6.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

Only 100 of 175 total aligned reads displayed.

GAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCC  >  NC_000913/1338463‑1338642

gAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCt                                                                                         <  1:4472602/93‑1 (MQ=255)

 aGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTc                                                                                        <  2:456435/93‑1 (MQ=255)

  gTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCg                                                                                       >  1:4524633/1‑93 (MQ=255)

   ttCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGt                                                                                      <  1:2245589/93‑1 (MQ=255)

    tCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTc                                                                                     <  1:3430580/93‑1 (MQ=255)

       tgtTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAgc                                                                                  >  1:2856017/1‑93 (MQ=255)

        gtTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAgcg                                                                                 >  1:1092576/1‑93 (MQ=255)

        gtTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAgcg                                                                                 >  1:3213814/1‑93 (MQ=255)

         tttCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCTACTGACGTGGCTGTAGTGCCGCGCTCGTCAgcgg                                                                                >  1:1891509/1‑93 (MQ=255)

         tttCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAgcgg                                                                                >  2:3548803/1‑93 (MQ=255)

          ttCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTGGTGCCGCGCTCGTCAgcggc                                                                               <  1:1032109/93‑1 (MQ=255)

          ttCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAgcggc                                                                               >  1:1152525/1‑93 (MQ=255)

           tCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAgcggcg                                                                              <  2:3553983/93‑1 (MQ=255)

             ccTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGaa                                                                            <  1:436787/93‑1 (MQ=255)

             ccTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGaa                                                                            <  1:1689365/93‑1 (MQ=255)

             ccTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCCACTGACGTGGCTGCAGTGCCCCGCTCGCCAGCGGCCaa                                                                            >  2:3202663/1‑93 (MQ=255)

              cTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAg                                                                           <  1:500563/93‑1 (MQ=255)

                ggCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGa                                                                         >  2:4042251/1‑93 (MQ=255)

                     gcgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAAc                                                                    >  1:1188833/1‑93 (MQ=255)

                      cgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACt                                                                   >  1:3629307/1‑93 (MQ=255)

                      cgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACt                                                                   <  1:4595146/93‑1 (MQ=255)

                       gcTGGACGTAAAATAAAGCGACGGGAAAACCAGTTTTCCCCCGCGACTGGGGTGGCTGTAGTGCCGCGCTTCTCAGCGGCGAAGTAATAACtt                                                                  <  2:893008/93‑1 (MQ=255)

                       gcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACtt                                                                  <  2:526750/93‑1 (MQ=255)

                         tGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACtttt                                                                >  1:1209793/1‑92 (MQ=255)

                         tGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTc                                                                <  2:3927771/93‑1 (MQ=255)

                         tGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTc                                                                <  2:2666153/93‑1 (MQ=255)

                          gAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCt                                                               >  1:3369519/1‑93 (MQ=255)

                           aaCGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTa                                                              >  1:2582993/1‑93 (MQ=255)

                           aaCGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTa                                                              >  1:2582994/1‑93 (MQ=255)

                              gTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGc                                                           >  2:2472288/1‑93 (MQ=255)

                               tAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCt                                                          <  1:474426/93‑1 (MQ=255)

                                aaaaTAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTg                                                         <  2:4467618/93‑1 (MQ=255)

                                   aTAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAg                                                      <  1:1843836/93‑1 (MQ=255)

                                   aTAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAg                                                      <  1:388969/93‑1 (MQ=255)

                                   aTAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAg                                                      <  2:2217423/93‑1 (MQ=255)

                                   aTAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAg                                                      <  2:179424/93‑1 (MQ=255)

                                      aaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttg                                                   >  2:3477836/1‑93 (MQ=255)

                                      aaGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttg                                                   <  2:4681711/93‑1 (MQ=255)

                                       aGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgt                                                  <  1:4312089/93‑1 (MQ=255)

                                       aGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgt                                                  <  2:3261421/93‑1 (MQ=255)

                                        gCGACTGGAAAACCAGCTTTCCCCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgtt                                                 <  2:3347313/93‑1 (MQ=255)

                                        gCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgtt                                                 >  2:2397389/1‑93 (MQ=255)

                                         cGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttg                                                <  1:971485/93‑1 (MQ=255)

                                         cGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttg                                                >  1:2813743/1‑93 (MQ=255)

                                          gACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttgt                                               >  2:3126587/1‑93 (MQ=255)

                                          gACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttgt                                               >  1:739281/1‑93 (MQ=255)

                                          gACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttgt                                               >  2:4377336/1‑93 (MQ=255)

                                          gACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttgt                                               >  1:2906648/1‑93 (MQ=255)

                                           aCTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttgtt                                              >  1:54828/1‑93 (MQ=255)

                                           aCTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAggttttgtt                                              >  2:2837684/1‑93 (MQ=255)

                                            cTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGttt                                             >  1:816142/1‑93 (MQ=255)

                                            cTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGttt                                             <  2:2313025/93‑1 (MQ=255)

                                             tGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTc                                            <  1:2610814/93‑1 (MQ=255)

                                             tGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTc                                            <  1:2365778/93‑1 (MQ=255)

                                              ggAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCt                                           >  2:3400041/1‑93 (MQ=255)

                                                 aaaCCAGCTTTCGCCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCt                                        <  1:1406407/93‑1 (MQ=255)

                                                 aaaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCt                                        <  2:3676883/93‑1 (MQ=255)

                                                  aaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCtt                                       >  1:3657615/1‑93 (MQ=255)

                                                  aaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCtt                                       <  2:3852303/93‑1 (MQ=255)

                                                  aaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCtt                                       <  2:2309741/93‑1 (MQ=255)

                                                  aaCCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGCTTCTGCtt                                       >  1:4359301/1‑93 (MQ=255)

                                                     cAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTg                                    >  2:1441089/1‑93 (MQ=255)

                                                      aGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGc                                   >  1:2042803/1‑93 (MQ=255)

                                                      aGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGc                                   <  2:2002390/93‑1 (MQ=255)

                                                        cTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCt                                 <  1:4207585/93‑1 (MQ=255)

                                                          ttCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGt                               <  1:2713528/93‑1 (MQ=255)

                                                           tCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                              <  2:2020519/93‑1 (MQ=255)

                                                           tCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                              <  2:123572/93‑1 (MQ=255)

                                                           tCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTa                              <  2:4272934/93‑1 (MQ=255)

                                                           tCACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTCCTTTTGCCTGtg                              >  1:230817/1‑92 (MQ=255)

                                                            cACCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAg                             >  1:3010483/1‑93 (MQ=255)

                                                             ccccGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGc                            <  2:2743088/92‑1 (MQ=255)

                                                             aCCCGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGc                            >  1:3293060/1‑93 (MQ=255)

                                                              cccGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGcc                           >  1:1823286/1‑93 (MQ=255)

                                                              cccGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGcc                           >  1:3812350/1‑93 (MQ=255)

                                                              cccGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGcc                           <  2:4524633/93‑1 (MQ=255)

                                                              cccGCGACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGcc                           >  1:1394981/1‑93 (MQ=255)

                                                                 gcgACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCt                        >  2:1813007/1‑93 (MQ=255)

                                                                 gcgACTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCt                        >  2:995813/1‑93 (MQ=255)

                                                                    aCTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGcc                     >  1:525767/1‑93 (MQ=255)

                                                                     cTGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCt                    >  2:1219840/1‑93 (MQ=255)

                                                                      tGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTa                   <  1:1687891/93‑1 (MQ=255)

                                                                      tGACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTa                   >  1:852008/1‑93 (MQ=255)

                                                                       gACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTAt                  >  1:3966650/1‑93 (MQ=255)

                                                                       gACGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTAt                  <  2:1801699/93‑1 (MQ=255)

                                                                        atgTGGCTGTAGTGCCTCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATg                 <  1:2166617/91‑1 (MQ=255)

                                                                        aCGTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATg                 <  1:1280229/93‑1 (MQ=255)

                                                                        aCGTGGATGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATg                 <  1:3148140/93‑1 (MQ=255)

                                                                          gTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGc               <  2:2448342/93‑1 (MQ=255)

                                                                          gTGGCTGTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGc               <  1:1633958/93‑1 (MQ=255)

                                                                             gCTGTAGTGCCGCGCTCGTCAGCCGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTgg            <  2:2582183/93‑1 (MQ=255)

                                                                                gTAGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGtat         <  2:1853293/93‑1 (MQ=255)

                                                                                  aGTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTatat       <  1:4720781/93‑1 (MQ=255)

                                                                                   gTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATg      <  1:2464798/93‑1 (MQ=255)

                                                                                   gTGCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGGATATg      >  2:3772167/1‑93 (MQ=255)

                                                                                    tGCCGCGCTCGGCAGCTGCGAAGGAATAACTGGCGATGTTGGAGTTGTTGTTCCTGCTGTTGCCTGTAGCCGCTGCCTATGGCTGGTATATgg     <  2:2422544/93‑1 (MQ=255)

                                                                                     gCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATggg    <  2:2337097/93‑1 (MQ=255)

                                                                                     gCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATggg    <  2:1917104/93‑1 (MQ=255)

                                                                                     gCCGCGCTCGTCAGCGGCGAAGGAATAACTTTCGATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATggg    <  1:4582654/93‑1 (MQ=255)

                                                                                      ccGCGCTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGc   <  1:600349/93‑1 (MQ=255)

GAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCC  >  NC_000913/1338463‑1338642

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