BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	734,377	T→C	5.3%	pseudogene (245/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,377	0	T	C	5.3%	465.5 / 11.8	151	pseudogene (245/1521 nt)	rhsO	pseudogene, Rhs family
Reads supporting (aligned to +/- strand): ref base T (72/71); new base C (5/3); total (77/74)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 152 total aligned reads displayed.

AACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTG  >  NC_000913/734232‑734525

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

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

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

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

  ccGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCAGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCGCCGGTACCATTACGACAGTCAGCACCGGCTGGTGc                                                                                                                                                 <  2:279865/149‑1 (MQ=17)

          cccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcg                                                                                                                                         >  2:941714/1‑149 (MQ=32)

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

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

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

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

                   cgcACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACaa                                                                                                                                <  2:697034/149‑1 (MQ=32)

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

                       cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATg                                                                                                                            <  3:158698/149‑1 (MQ=32)

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

                        tATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGc                                                                                                                           <  4:74730/149‑1 (MQ=32)

                        tATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCAGCCGGTGctga                                                                                                                                              <  3:216647/130‑4 (MQ=255)

                          tCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGCGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCag                                                                                                                         <  1:820929/149‑1 (MQ=17)

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

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

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

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

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

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

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

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

                                           gTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGc                                                                                                        <  1:531300/149‑1 (MQ=32)

                                                   aGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCt                                                                                                <  3:564798/149‑1 (MQ=34)

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

                                                     gCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCttt                                                                                              <  2:87735/149‑1 (MQ=35)

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

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

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

                                                       tGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTAc                                                                                            >  1:396709/1‑149 (MQ=34)

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

                                                          cGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGAc                                                                                         <  4:501337/149‑1 (MQ=33)

                                                           ggAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGAcc                                                                                        >  3:866611/1‑149 (MQ=32)

                                                            gagaAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGAccc                                                                                       <  4:596999/149‑1 (MQ=32)

                                                                aaaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTg                                                                                   <  2:65643/149‑1 (MQ=32)

                                                                 aaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAATCCGCTATCTTTa                                                                                             >  4:796557/1‑138 (MQ=17)

                                                                 aaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                                                                  >  4:902783/1‑149 (MQ=32)

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

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

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

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

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

                                                                         tCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAgg                                                                          >  3:68218/1‑149 (MQ=32)

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

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

                                                                           aTCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGAGAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCa                                                                            <  3:793765/145‑1 (MQ=17)

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

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

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

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

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

                                                                                  aaGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaaaa                                                                 <  2:191189/149‑1 (MQ=32)

                                                                                   aGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAAc                                                                >  3:373791/1‑149 (MQ=32)

                                                                                    ggggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACg                                                               >  4:443923/1‑149 (MQ=32)

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

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

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

                                                                                           tCCGCACGGATGATGAGCGCACCCGCCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGCACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTgg                                                        <  2:749641/149‑1 (MQ=21)

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

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

                                                                                           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)

                                                                                                  ggATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTg                                                 >  3:987367/1‑149 (MQ=38)

                                                                                                         gaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGa                                          <  4:713626/149‑1 (MQ=38)

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

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

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

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

                                                                                                                 ccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACgg                                  >  2:388849/1‑149 (MQ=39)

                                                                                                                    ccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCt                               <  1:616923/149‑1 (MQ=39)

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

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

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

                                                                                                                              tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCt                     <  3:443923/149‑1 (MQ=38)

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

                                                                                                                                 gACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGa                                          <  1:267323/125‑1 (MQ=38)

                                                                                                                                 gACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTc                  <  1:689139/149‑1 (MQ=38)

                                                                                                                                  aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCa                 <  4:163186/149‑1 (MQ=37)

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

                                                                                                                                   cAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCAc                >  3:85040/1‑149 (MQ=38)

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

                                                                                                                                      agtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGa             >  4:681317/4‑149 (MQ=255)

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

                                                                                                                                       gtgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGaa            >  3:527522/3‑149 (MQ=255)

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

                                                                                                                                        aGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGaaa           <  1:46967/149‑1 (MQ=38)

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

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

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

                                                                                                                                           aCCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCg        >  4:206127/1‑149 (MQ=38)

                                                                                                                                            ccGGCTGGTGCACTACACGCGGAAACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGc       <  3:613909/149‑1 (MQ=37)

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

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

                                                                                                                                               gCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTAc                                                                                            <  1:184612/61‑1 (MQ=35)

                                                                                                                                                 tGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTg  <  1:254959/149‑1 (MQ=38)

AACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTG  >  NC_000913/734232‑734525

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