BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,244,370	C→T	11.3%	G39S (GGT→AGT) ‡	rnt ←	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,244,370	0	C	T	11.3%	231.1 / 15.5	88	G39S (GGT→AGT) ‡	rnt	ribonuclease T
Reads supporting (aligned to +/- strand): ref base C (42/36); new base T (7/3); total (49/39)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTT  >  NC_002947/1244222‑1244518

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAaccacc                                                                                                                                                      >  3:93953/1‑149 (MQ=255)

tCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGT‑‑CGGCGGCAATATACACCTATCGAAAGGTGGAGCTCTTTAAACCACCgg                                                                                                                                                    >  3:308237/1‑149 (MQ=255)

  aaGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCgg                                                                                                                                                    >  8:288709/1‑149 (MQ=255)

        ggccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCGGTGGTTTCaa                                                                                                                                              <  5:332794/149‑1 (MQ=255)

        ggccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCaa                                                                                                                                              >  3:402525/1‑149 (MQ=255)

        ggccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCaa                                                                                                                                              <  8:23556/149‑1 (MQ=255)

          ccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCGGTGGTTTCAACa                                                                                                                                            <  1:92630/149‑1 (MQ=255)

          ccggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACa                                                                                                                                            >  6:158443/1‑149 (MQ=255)

            ggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATc                                                                                                                                          >  1:362650/1‑149 (MQ=255)

            ggcTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATc                                                                                                                                          <  5:60667/149‑1 (MQ=255)

               tCGATCTTGGAGCCCTCTCAGTGCTCCACTCGTTGGAAGTATGTATGTTCTGGGAACAGTATACCATTCTCGTCAATGGCAATGGTGTGGGCTTCATTTTCAAGCAATGGATCTGTAGCGCTTTTTAAACCACCGGTTTCAACATcctc                                                                                                                                       <  4:147885/149‑3 (MQ=255)

                cGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCGATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCa                                                                                                                                               <  6:314474/140‑1 (MQ=255)

                cGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCGATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCa                                                                                                                                               >  5:314474/1‑140 (MQ=255)

                       ggCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCCCGACCg                                                                                                                               >  7:251582/1‑149 (MQ=255)

                       ggCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGAcc                                                                                                                                <  3:296429/148‑1 (MQ=255)

                       ggCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGGGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGAcc                                                                                                                                <  4:330527/148‑1 (MQ=255)

                          gcCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCgg                                                                                                                              >  1:260537/1‑147 (MQ=255)

                            ccTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCAACCACGACCGGCAg                                                                                                                           >  3:235446/1‑148 (MQ=255)

                             cTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGAt                                                                                                                         <  4:93953/149‑1 (MQ=255)

                              tCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGa                                                                                                                          >  6:244493/1‑147 (MQ=255)

                                gAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAaccacc                                                                                                                                                      <  6:215334/117‑1 (MQ=255)

                                gAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAaccacc                                                                                                                                                      >  5:215334/1‑117 (MQ=255)

                                gAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATa                                                                                                                        <  7:283221/147‑1 (MQ=255)

                                gAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCACCAATGCATCGGTAGCGCTGTTGAAACCACCGGTATCAACATCCACCACGAGCGGCAGATAGc                                                                                                                      >  6:313034/1‑149 (MQ=255)

                                 aaCGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGcc                                                                                                                     <  2:260537/149‑1 (MQ=255)

                                      cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCcg                                                                                                                    <  1:269935/145‑1 (MQ=255)

                                      cTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCcg                                                                                                                    >  2:269935/1‑145 (MQ=255)

                                         cACTCGGTGGAAGTAGGTATGTTCGGGGAACATGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGaa                                                                                                               <  7:363780/147‑1 (MQ=255)

                                         cACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAgc                                                                                                             >  6:258428/1‑149 (MQ=255)

                                            tCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCt                                                                                                          >  5:255354/1‑149 (MQ=255)

                                                   aaGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCAt                                                                                                   >  2:157054/1‑149 (MQ=255)

                                                    aGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATc                                                                                                  <  5:399961/149‑1 (MQ=255)

                                                    aGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATc                                                                                                  <  1:399811/149‑1 (MQ=255)

                                                     gTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCg                                                                                                 <  6:255354/149‑1 (MQ=255)

                                                        ggTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATTGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGt                                                                                              <  8:140948/149‑1 (MQ=255)

                                                        ggTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGt                                                                                              <  7:288709/149‑1 (MQ=255)

                                                          tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTgg                                                                                            >  2:405003/1‑149 (MQ=255)

                                                          tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTgg                                                                                            <  8:371060/149‑1 (MQ=255)

                                                              ttCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAg                                                                                        <  5:147400/149‑1 (MQ=255)

                                                                 ggggAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGGCGCGGAAGCGCTCGGCGATCGGGTGGCGAGGAc                                                                                     <  4:235446/149‑1 (MQ=255)

                                                                   ggAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCg                                                                                   >  1:340626/1‑149 (MQ=255)

                                                                     aaCAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgct                                                                                 <  1:157054/149‑1 (MQ=255)

                                                                        aGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgctgct                                                                              <  2:340626/149‑1 (MQ=255)

                                                                          gAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCtgctgc                                                                            <  4:48611/149‑1 (MQ=255)

                                                                          gAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCtgctgc                                                                            >  3:48611/1‑149 (MQ=255)

                                                                           aaaCCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGcc                                                                           >  6:330411/1‑149 (MQ=255)

                                                                               cTTTCTCCTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTgg                                                                       <  4:402525/149‑1 (MQ=255)

                                                                                tttCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTggg                                                                      >  5:98120/1‑149 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTc                                                                  <  5:330411/147‑1 (MQ=255)

                                                                                      gTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTAGGGTGTc                                                                  >  7:93571/1‑147 (MQ=255)

                                                                                       tCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAg                                                               <  5:258428/149‑1 (MQ=255)

                                                                                          aTGCCAATGGTGATGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtc                                                            <  4:388159/149‑1 (MQ=255)

                                                                                          aTGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtc                                                            >  2:244854/1‑149 (MQ=255)

                                                                                              cAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtcgtc                                                         >  6:203591/1‑148 (MQ=255)

                                                                                              cAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtcgtc                                                         <  5:203591/148‑1 (MQ=255)

                                                                                                 tGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCg                                                     >  4:58021/1‑149 (MQ=255)

                                                                                                   gTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTc                                                                  <  3:32732/134‑1 (MQ=255)

                                                                                                   gTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTc                                                                  >  4:32732/1‑134 (MQ=255)

                                                                                                      aCGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                <  5:244493/149‑1 (MQ=255)

                                                                                                      aCGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                <  8:93571/149‑1 (MQ=255)

                                                                                                           ggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgctg                                                                                >  5:19472/1‑112 (MQ=255)

                                                                                                           ggcAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCgctg                                                                                <  6:19472/112‑1 (MQ=255)

                                                                                                             cAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCg                                         >  3:239155/1‑149 (MQ=255)

                                                                                                              aaTTTCCAGCAATGCATCGGTAGCGCTGATGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGc                                        >  1:227851/1‑149 (MQ=255)

                                                                                                                tttCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTc                                      <  5:158443/149‑1 (MQ=255)

                                                                                                                     aGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCAcgcg                                 >  3:49001/1‑149 (MQ=255)

                                                                                                                       cAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGtt                               >  6:153028/1‑149 (MQ=255)

                                                                                                                        aaTGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGttt                              >  8:221062/1‑149 (MQ=255)

                                                                                                                             aTCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGa                                                                                         <  6:395150/85‑1 (MQ=255)

                                                                                                                                 gTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCcag                      >  2:119546/1‑148 (MQ=255)

                                                                                                                                  tAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCcag                      >  7:98479/1‑147 (MQ=255)

                                                                                                                                     cgcTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCca                       >  8:185274/1‑143 (MQ=255)

                                                                                                                                     cgcTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCca                       <  7:185274/143‑1 (MQ=255)

                                                                                                                                        tGTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCTCTCACGCGTTTTTCTCCAGCAGgcgc               >  4:90308/1‑148 (MQ=255)

                                                                                                                                        tGTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgc               >  3:27390/1‑148 (MQ=255)

                                                                                                                                        tGTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGcc              >  8:134702/1‑149 (MQ=255)

                                                                                                                                        tGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGcc              >  5:177368/1‑149 (MQ=255)

                                                                                                                                         gTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgc               >  4:106008/1‑147 (MQ=255)

                                                                                                                                         gTTGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCCGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgc               >  8:303771/1‑147 (MQ=255)

                                                                                                                                         gTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGgcgc               >  3:269612/1‑147 (MQ=255)

                                                                                                                                          ttGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                <  3:306245/113‑1 (MQ=255)

                                                                                                                                          ttGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAg                                                >  4:306245/1‑113 (MQ=255)

                                                                                                                                          ttGAAACCGGTGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAg            >  7:23854/1‑149 (MQ=255)

                                                                                                                                          ttGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCTGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAg            <  2:227851/149‑1 (MQ=255)

                                                                                                                                          ttGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAg            >  6:368823/1‑149 (MQ=255)

                                                                                                                                          ttGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGCAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAg            >  1:290578/1‑149 (MQ=255)

                                                                                                                                            gAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTcc                        >  4:204520/1‑135 (MQ=255)

                                                                                                                                            gAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTcc                        <  3:204520/135‑1 (MQ=255)

                                                                                                                                            gAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgcg          >  7:320639/1‑149 (MQ=255)

                                                                                                                                              aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgcg          >  5:185629/1‑147 (MQ=255)

                                                                                                                                              aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgcg          >  7:206210/1‑147 (MQ=255)

                                                                                                                                                ccaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTGCAGGTCGTCTTCGTAAAGGTCt                                            <  4:95810/111‑1 (MQ=255)

                                                                                                                                                ccaccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCt                                            >  3:95810/1‑111 (MQ=255)

                                                                                                                                                  accGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGttt                              <  5:154842/123‑1 (MQ=255)

                                                                                                                                                  accGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGttt                              >  6:154842/1‑123 (MQ=255)

                                                                                                                                                   ccGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGtttt   >  4:126815/1‑149 (MQ=255)

                                                                                                                                                    cGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGttttt  >  8:11243/1‑149 (MQ=255)

TCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCGCAGTTTTT  >  NC_002947/1244222‑1244518

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