BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	734,374	G→A	5.4%	pseudogene (242/1521 nt)	rhsO →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	734,374	0	G	A	5.4%	449.1 / 12.3	147	pseudogene (242/1521 nt)	rhsO	pseudogene, Rhs family
Reads supporting (aligned to +/- strand): ref base G (71/68); new base A (5/3); total (76/71)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

Only 100 of 150 total aligned reads displayed.

GGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAA  >  NC_000913/734228‑734522

ggATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGc                                                                                                                                                    <  4:14761/149‑1 (MQ=32)

    aaCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGt                                                                                                                                                <  3:807207/149‑1 (MQ=32)

     aCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                                                               >  1:112810/1‑149 (MQ=32)

      ccGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGc                                                                                                                                              >  2:14589/1‑149 (MQ=32)

      ccGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGc                                                                                                                                              <  4:483032/149‑1 (MQ=32)

       cGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCa                                                                                                                                             >  4:197356/1‑149 (MQ=32)

         tATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctga                                                                                                                                           <  2:719579/149‑4 (MQ=12)

               ccGTGACGTGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcgc                                                                                                                                     <  3:546345/149‑1 (MQ=17)

               ccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcgc                                                                                                                                     >  3:676209/1‑149 (MQ=32)

                 gTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                   >  4:564798/1‑149 (MQ=32)

                  tGACGCGCACTATCTTTACCGGTATGACGGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGa                                                                                                                                  <  4:1081967/149‑1 (MQ=17)

                  tGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGa                                                                                                                                  >  2:283760/1‑149 (MQ=32)

                   gACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGac                                                                                                                                 >  1:384355/1‑149 (MQ=32)

                      gcgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGAcaca                                                                                                                              >  3:501337/1‑149 (MQ=32)

                          aCTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAatat                                                                                                                          >  1:1098023/1‑149 (MQ=32)

                           cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATg                                                                                                                         >  1:455667/1‑149 (MQ=32)

                            tATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAata                                                                                                                           >  1:486923/1‑146 (MQ=32)

                            tATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctga                                                                                                                                           >  4:216647/1‑127 (MQ=12)

                                tttACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCagag                                                                                                                    >  4:861318/1‑149 (MQ=32)

                                tttACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCagag                                                                                                                    >  3:858504/1‑149 (MQ=32)

                                 ttACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAatat                                                                                                                          >  4:132120/1‑142 (MQ=32)

                                 ttACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAatat                                                                                                                          <  3:132120/142‑1 (MQ=32)

                                       gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTgg                                                                                                             <  3:947606/149‑1 (MQ=32)

                                       gTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTgg                                                                                                             >  3:596999/1‑149 (MQ=32)

                                        tATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGt                                                                                                            <  4:22014/149‑1 (MQ=32)

                                        tATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGt                                                                                                            >  1:550791/1‑149 (MQ=32)

                                         aTGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTc                                                                                                           <  2:324988/149‑1 (MQ=32)

                                           gACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGa                                                                                                         >  2:254959/1‑149 (MQ=32)

                                            aCCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGaa                                                                                                        >  1:65643/1‑149 (MQ=32)

                                             ccGTCACTGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGaaa                                                                                                       >  2:616923/1‑149 (MQ=17)

                                               gTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGc                                                                                                     >  1:659336/1‑149 (MQ=33)

                                                        ggCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCtt                                                                                            <  2:1098023/149‑1 (MQ=34)

                                                         gCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCttt                                                                                           <  2:112810/149‑1 (MQ=34)

                                                           tGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTAc                                                                                         <  2:550791/149‑1 (MQ=35)

                                                           tGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTAc                                                                                         >  2:783676/1‑149 (MQ=35)

                                                            gACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACg                                                                                        >  1:893493/1‑149 (MQ=35)

                                                                     aaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                               >  3:109313/1‑149 (MQ=32)

                                                                     aaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGACAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTa                                                                                          <  3:796557/138‑1 (MQ=18)

                                                                       ccGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGc                                                                             <  2:1077170/149‑1 (MQ=32)

                                                                        cGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                               >  1:708125/1‑146 (MQ=32)

                                                                        cGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                               <  2:708125/146‑1 (MQ=32)

                                                                        cGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGcc                                                                            >  3:898697/1‑149 (MQ=32)

                                                                         gACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCGGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCg                                                                           >  3:751342/1‑149 (MQ=17)

                                                                             tCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAgg                                                                       >  4:197637/1‑149 (MQ=32)

                                                                              cATCCCGGAAGGGGTTATCTGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAggg                                                                      >  4:723202/1‑149 (MQ=17)

                                                                              cATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAggg                                                                      <  2:659336/149‑1 (MQ=32)

                                                                               aTCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTATACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                     >  4:595575/1‑149 (MQ=17)

                                                                               aTCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                     <  4:127881/149‑1 (MQ=32)

                                                                               aTCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACAGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCa                                                                         >  4:793765/1‑145 (MQ=17)

                                                                                tCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTg                                                                    >  4:979383/1‑149 (MQ=32)

                                                                                tCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTg                                                                    <  4:858504/149‑1 (MQ=32)

                                                                                tCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGGACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTg                                                                    >  1:559353/1‑149 (MQ=17)

                                                                                  ccGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGc                                                                  >  1:191189/1‑149 (MQ=32)

                                                                                  ccGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGc                                                                  <  4:68218/149‑1 (MQ=32)

                                                                                    ggAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaa                                                                >  3:280528/1‑149 (MQ=32)

                                                                                        ggggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCACTGGGCCGCAGGGTGGCAAAACg                                                            >  3:157219/1‑149 (MQ=17)

                                                                                         gggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCGGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                           <  4:751342/149‑1 (MQ=17)

                                                                                         gggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                           >  4:243799/1‑149 (MQ=32)

                                                                                           gTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCACTGGGCCGCAGGGTGGCAAAACGGgt                                                         <  4:157219/149‑1 (MQ=17)

                                                                                               tCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                               >  1:172206/1‑123 (MQ=32)

                                                                                               tCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                               <  2:172206/123‑1 (MQ=33)

                                                                                               tCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTgg                                                     >  3:739155/1‑149 (MQ=34)

                                                                                               tCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGATGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTgg                                                     <  1:142411/149‑1 (MQ=21)

                                                                                                  gCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGcc                                                                                                    >  1:12815/1‑99 (MQ=32)

                                                                                                   cACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACTCGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGAc                                                 <  3:4416/149‑1 (MQ=37)

                                                                                                   cACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGACAAAACGGGTGTGGCGAc                                                 <  4:676209/149‑1 (MQ=25)

                                                                                                             gaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGATGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGa                                       <  4:704113/149‑1 (MQ=37)

                                                                                                               gcgcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGAcc                                     <  1:941714/149‑1 (MQ=39)

                                                                                                                   acccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGAc                                 <  4:109313/149‑1 (MQ=39)

                                                                                                                     ccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACgg                               >  3:163186/1‑149 (MQ=38)

                                                                                                                     ccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACgg                               >  2:689139/1‑149 (MQ=38)

                                                                                                                     ccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACgg                               <  3:861318/149‑1 (MQ=38)

                                                                                                                     ccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACgg                               >  4:354405/1‑149 (MQ=38)

                                                                                                                       accGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGc                             >  3:387438/1‑149 (MQ=37)

                                                                                                                        ccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTATACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCt                            <  3:595575/149‑1 (MQ=37)

                                                                                                                         cGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTg                           <  3:243799/149‑1 (MQ=37)

                                                                                                                           gTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGa                         >  2:856333/1‑149 (MQ=37)

                                                                                                                            tACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGAt                        <  2:455667/149‑1 (MQ=39)

                                                                                                                                 ttACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCttt                                                                                           <  1:1058464/77‑1 (MQ=35)

                                                                                                                                 ttACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCttt                                                                                           >  2:1058464/1‑77 (MQ=34)

                                                                                                                                    cGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGTGACCTGACGGGCTGGATGTCGCt                  <  3:1143775/147‑1 (MQ=37)

                                                                                                                                    cGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGt                >  3:563124/1‑149 (MQ=37)

                                                                                                                                     gACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTc               >  1:703971/1‑149 (MQ=37)

                                                                                                                                     gACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCACAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGa                                       >  2:267323/1‑125 (MQ=38)

                                                                                                                                      aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCa              <  3:197356/149‑1 (MQ=37)

                                                                                                                                      aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCa              <  1:283487/149‑1 (MQ=39)

                                                                                                                                        gtagtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACg            >  4:673841/6‑149 (MQ=255)

                                                                                                                                          tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCt                            >  2:469046/1‑131 (MQ=38)

                                                                                                                                          tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCt                            <  1:469046/131‑1 (MQ=38)

                                                                                                                                          tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGa          >  4:450945/1‑149 (MQ=39)

                                                                                                                                          tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGa          >  2:925403/1‑149 (MQ=37)

                                                                                                                                          tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGa          >  3:65158/1‑149 (MQ=37)

                                                                                                                                           gtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                               >  2:881967/3‑79 (MQ=11)

                                                                                                                                           gtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                               <  1:881967/77‑1 (MQ=11)

                                                                                                                                            aGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCt                                    <  2:1149700/121‑1 (MQ=40)

                                                                                                                                             gCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAAc       <  4:739155/149‑1 (MQ=38)

                                                                                                                                             gCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAAc       >  1:514029/1‑149 (MQ=38)

                                                                                                                                              cACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAAcc      <  1:1023992/149‑1 (MQ=38)

                                                                                                                                                  ggCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCaa  <  3:450945/149‑1 (MQ=38)

                                                                                                                                                  ggCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCaa  <  4:987367/149‑1 (MQ=38)

GGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAA  >  NC_000913/734228‑734522

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