BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,244,368	A→G	11.4%	G39G (GGT→GGC) ‡	rnt ←	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,244,368	0	A	G	11.4%	219.8 / 17.7	88	G39G (GGT→GGC) ‡	rnt	ribonuclease T
Reads supporting (aligned to +/- strand): ref base A (42/36); new base G (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 = 5.19e-01

CTCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCG  >  NC_002947/1244221‑1244510

cTCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGACAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAaccac                                                                                                                                               >  8:283221/1‑149 (MQ=255)

cTCAAGAGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGTGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAaccac                                                                                                                                               <  6:106798/149‑1 (MQ=255)

 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                                                                                                                                  <  5:60667/149‑1 (MQ=255)

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

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

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

                 cGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCGATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCa                                                                                                                                       <  6:314474/140‑1 (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                                                                                    <  8:371060/149‑1 (MQ=255)

                                                           tATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTgg                                                                                    >  2:405003/1‑149 (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                                                 <  5:203591/148‑1 (MQ=255)

                                                                                               cAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGgtcgtc                                                 >  6:203591/1‑148 (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                <  3:204520/135‑1 (MQ=255)

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

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

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

                                                                                                                                               aaCCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAgcg  >  5:185629/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)

CTCAAGTGCGGCCGGCTCGATGTTGGCGCCCTCGAACGGCTCCACTCGGTGGAAGTAGGTATGTTCGGGGAACAGGAAACCTTTCTCGTCCATGCCAATGGTGACGGCGGCAATTTCCAGCAATGCATCGGTAGCGCTGTTGAAACCACCGGTTTCAACATCCACCACGACCGGCAGATAGCCGCGGAAGCGCTCGGCCATCGGGTGGCGAGGACCGCTGCTGCCTTGGGTGTCCAGGTCGTCTTCGTAAAGGTCTTCGCTCACGCGTTTTTCTCCAGCAGGCGCCAGCG  >  NC_002947/1244221‑1244510

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: `‑`