BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,707,320	A→T	5.9%	intergenic (‑196/+2)	lepA ← / ← rseC	back‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,707,320	0	A	T	5.9%	725.7 / 25.5	203	intergenic (‑196/+2)	lepA/rseC	back‑translocating elongation factor EF4, GTPase/SoxR iron‑sulfur cluster reduction factor component
Reads supporting (aligned to +/- strand): ref base A (109/82); new base T (3/9); total (112/91)
Fisher's exact test for biased strand distribution p-value = 3.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 208 total aligned reads displayed.

CGGTCAGACAGCGTCGATTTACCGTGGTCAATGTGAGCTATGATCGAAAAGTTACGTATATTCTTCATAAAGTATGATTATTGTGCCTTAATGCCCGGTTAACCAGGCTTTTAGAAGTCGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGG‑CGACATATAAACCAGTAATGCGGAGCTAAGCAGGC  >  NC_000913/2707057‑2707584

             tccaTTTGCCGTGGGCAAGTTATGCTATTATTGAAAAGTTACGTTAATTTTTCATAAAGTATGATTATTTTGCCTTAATGCCCGGTTAACCAGGCTTTTAGAAGTCGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAg                                                                                                                                                                                                                                                   <  4:322464/272‑1 (MQ=255)

                                                                                                                        cTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCAc                                                                                                                                                                                                                                                                      >  5:853875/1‑149 (MQ=255)

                                                                                                                              tGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtt                                                                                                                                                                                                                                                                <  6:862916/149‑2 (MQ=255)

                                                                                                                              tGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtt                                                                                                                                                                                                                                                                <  7:287011/149‑2 (MQ=255)

                                                                                                                              tGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTc                                                                                                                                                                                                                                                                >  6:1095370/1‑149 (MQ=255)

                                                                                                                               gAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCttt                                                                                                                                                                                                                                                               <  5:1314819/149‑3 (MQ=255)

                                                                                                                               gAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCttt                                                                                                                                                                                                                                                               <  8:469061/149‑3 (MQ=255)

                                                                                                                                       cGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACATAGCCAGTGAAATCACTGGCTCGCGTCTTcc                                                                                                                                                                                                                                                       <  5:719759/149‑1 (MQ=255)

                                                                                                                                        gTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCg                                                                                                                                                                                                                                                      >  7:577246/1‑149 (MQ=255)

                                                                                                                                        gTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCg                                                                                                                                                                                                                                                      <  7:23025/149‑1 (MQ=255)

                                                                                                                                          cTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGaa                                                                                                                                                                                                                                                    <  6:279154/149‑1 (MQ=255)

                                                                                                                                           tGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgt                                                                                                                                                                                                                                                             <  8:511380/139‑5 (MQ=255)

                                                                                                                                           tGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGATTTCACTGGCtttgt                                                                                                                                                                                                                                                             >  7:511380/1‑135 (MQ=255)

                                                                                                                                               ttAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGt                                                                                                                                                                                                                                               <  5:37798/149‑1 (MQ=255)

                                                                                                                                                tAATAGAAACGCCCATTTAAAAACTACAACGTTTGGCCGAGGAAATATTCATAACTTAGGGTTTGGGGTATCTGGAAACTTCCTCTCATTTGTTATGCAAAATGCAACAAAGCCAGTGAACTCCCGGGCTCCCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGca                                                                                                                <  4:238906/275‑1 (MQ=255)

                                                                                                                                                  aTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTc                                                                                                                                                                                                                                            >  8:699019/1‑149 (MQ=255)

                                                                                                                                                   tAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCa                                                                                                                                                                                                                                           >  8:1347546/1‑149 (MQ=255)

                                                                                                                                                   tAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCa                                                                                                                                                                                                                                           >  6:1196370/1‑149 (MQ=255)

                                                                                                                                                    aGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaa                                                                                                                                                                                                                                          >  7:1078968/1‑149 (MQ=255)

                                                                                                                                                    aGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaa                                                                                                                                                                                                                                          >  5:195390/1‑149 (MQ=255)

                                                                                                                                                    aGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaa                                                                                                                                                                                                                                          >  6:122457/1‑149 (MQ=255)

                                                                                                                                                     gAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCaaa                                                                                                                                                                                                                                         <  8:1224912/149‑1 (MQ=255)

                                                                                                                                                          gccgcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGca                                                                                                                                                                                                                                    <  8:93411/149‑1 (MQ=255)

                                                                                                                                                            cgcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCAAAAGGCCAGGCGGCAAGGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGAAGGAAACCTACGCAAGTAGCCGGCAGACATCAGGAACCCGCCAATACCACCGAGGATAGcacccaaatatgct                                                                                                    >  4:53234/1‑265 (MQ=255)

                                                                                                                                                             gcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGa                                                                                                                                                                                                                                                  <  7:475331/132‑1 (MQ=255)

                                                                                                                                                              cATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAg                                                                                                                                                                                                                                <  7:1436186/149‑1 (MQ=255)

                                                                                                                                                                   tacactacaACGTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAg                                                                                                                                                                                                                           >  5:771622/1‑149 (MQ=255)

                                                                                                                                                                   tacactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAg                                                                                                                                                                                                                           >  5:915294/1‑149 (MQ=255)

                                                                                                                                                                       ctacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggcgg                                                                                                                                                                                                                       >  7:704594/1‑149 (MQ=255)

                                                                                                                                                                           aacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAAt                                                                                                                                                                                                                   <  5:121539/149‑1 (MQ=255)

                                                                                                                                                                              aTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGcc                                                                                                                                                                                                                >  6:127471/1‑149 (MQ=255)

                                                                                                                                                                              aTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGcc                                                                                                                                                                                                                >  7:156830/1‑149 (MQ=255)

                                                                                                                                                                               ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCa                                                                                                                                                                                                               >  8:590652/1‑149 (MQ=255)

                                                                                                                                                                                tGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCAc                                                                                                                                                                                                              <  6:915294/149‑1 (MQ=255)

                                                                                                                                                                                  aGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGc                                                                                                                                                                                                            <  6:1356775/149‑1 (MQ=255)

                                                                                                                                                                                      gAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTaa                                                                                                                                                                                                        <  7:965662/149‑1 (MQ=255)

                                                                                                                                                                                       aGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGACCGCAGCACATAATGCTGCGACGTCAGAAGAAAAGAGTAACTGa                                                                         >  4:35092/1‑275 (MQ=255)

                                                                                                                                                                                       aGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAAGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATCACACCGCGGATCGCACCACATAAGGCTGCGACGTCAGAAGCAGAGAGAAGCGGa                                                                         >  1:8165/1‑275 (MQ=255)

                                                                                                                                                                                        ggAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAga                                                                                                                                                                                                      >  7:763061/1‑149 (MQ=255)

                                                                                                                                                                                             tGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATc                                                                                                                                                                                                 <  7:276040/149‑1 (MQ=255)

                                                                                                                                                                                              gTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCg                                                                                                                                                                                                >  6:1179062/1‑149 (MQ=255)

                                                                                                                                                                                                tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGt                                                                                                                                                                                              >  8:1080022/1‑149 (MQ=255)

                                                                                                                                                                                                      cGTATGGATTGTGGTATCTGGAAACTTCCTCCCATTTGTTAGGCAAAATGCAACCAAGCCATTGAAATCACTGGCTGGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGAt                                                           <  2:53234/274‑1 (MQ=255)

                                                                                                                                                                                                       gTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGt                                                                                                                                                                                                                                               >  5:444028/1‑93 (MQ=255)

                                                                                                                                                                                                        tATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCcgcgcg                                                                                                                                                   <  4:145683/184‑1 (MQ=255)

                                                                                                                                                                                                        tATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAg                                                        >  2:174407/1‑275 (MQ=255)

                                                                                                                                                                                                        tATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGACCGCACCACCTAATGCGTCGACGGCAGAAACAAAGAGTCGCTAAAATagccaacggaaaa                                                         >  4:278963/1‑267 (MQ=255)

                                                                                                                                                                                                         aTGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGCTGATCGGTTGCCATTCTGCCCGCGCAGCAAACTTACGCGAGTAGCCGCGCTCAATCAGGAACCCGCCAATACCACCGACGAACGCACCACATAATGCTGCGACGTCCGAAGCAAAGAGTAGCTGCAATAGCGAAGCGGTAAg                                                        >  1:64898/1‑274 (MQ=255)

                                                                                                                                                                                                         aTGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACg                                                                                              >  4:88982/1‑236 (MQ=255)

                                                                                                                                                                                                           ggATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGc                                                                                                                                                                                   <  8:531733/149‑1 (MQ=255)

                                                                                                                                                                                                            gagtgtgGTAGCTGGAAACTTCCTCGCATTTGTTATGCAAAATGCAACAAGGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGAt                                                           <  2:278963/265‑1 (MQ=255)

                                                                                                                                                                                                             aTTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCcg                                                                                                                                                       >  2:39957/1‑175 (MQ=255)

                                                                                                                                                                                                             aTTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAg                                                                                          <  2:88982/236‑1 (MQ=255)

                                                                                                                                                                                                              ttGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCg                                                                                                                                                                                <  5:1191815/149‑1 (MQ=255)

                                                                                                                                                                                                                 tgGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGAGCGCACCACATAATGCTGCGACGTCAGAAGCAAAGCGTAGCTGAAATagcgaagc                                                              >  1:44087/1‑260 (MQ=255)

                                                                                                                                                                                                                       cTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCaaa                                                                                                                                                                       <  8:128805/149‑1 (MQ=255)

                                                                                                                                                                                                                        tGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAAc                                                                                                                                                                      >  6:185132/1‑149 (MQ=255)

                                                                                                                                                                                                                         ggAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACt                                                                                                                                                                     >  6:299882/1‑149 (MQ=255)

                                                                                                                                                                                                                            aaCGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAAGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAc                                                                                                                                                                  <  7:1208979/149‑1 (MQ=255)

                                                                                                                                                                                                                            aaCGCCCGCCCATTTTTTATGCAAAATGCAACAAAGCGAGTGAAACCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTGGCCATTCTGCCCGGGCAGCAAACTTACGCAGGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAgcgg‑cg                                   <  2:35092/275‑1 (MQ=255)

                                                                                                                                                                                                                              cGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcgc                                                                                                                                                                <  8:419434/149‑1 (MQ=255)

                                                                                                                                                                                                                               gTCCTCGCATTTTTTATGCAAAATGCAACAAAGCGAGTGAAATCACTGGCACGCGTCTTCGGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGGTGATCGGTGGCAATTCTGCCCGGGCAGCAAACTTACGCGAGTAGAcgcg                                                                                                                                                     <  2:133537/159‑1 (MQ=255)

                                                                                                                                                                                                                                tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCCCAAGGCCAGGCGGCAATGCCACGCTGAAGATGATCGGTGGCCATTCTGCCCGGGCAGCAAACTTACGCGa                                                                                                                                                              <  8:358269/149‑1 (MQ=255)

                                                                                                                                                                                                                                tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGa                                                                                                                                                              <  7:1083322/149‑1 (MQ=255)

                                                                                                                                                                                                                                tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGa                                                                                                                                                              <  5:1415209/149‑1 (MQ=255)

                                                                                                                                                                                                                                tCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGACAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTAcgcaa                                                                                                                                                              <  8:760034/149‑3 (MQ=255)

                                                                                                                                                                                                                                 ccTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAg                                                                                                                                                             >  8:494439/1‑149 (MQ=255)

                                                                                                                                                                                                                                  cTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGt                                                                                                                                                            >  7:1001168/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGTGGCAATGCCACGCTTAAga                                                                                                                                                                                                      <  8:1007277/106‑1 (MQ=255)

                                                                                                                                                                                                                                   tCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTa                                                                                                                                                           >  5:1444646/1‑149 (MQ=255)

                                                                                                                                                                                                                                   tCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTa                                                                                                                                                           >  7:1306346/1‑149 (MQ=255)

                                                                                                                                                                                                                                    cGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAg                                                                                                                                                          <  8:778233/149‑1 (MQ=255)

                                                                                                                                                                                                                                      cATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGcc                                                                                                                                                        <  6:554514/149‑1 (MQ=255)

                                                                                                                                                                                                                                        tttGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAgatga                                                                                                                                                                                                   <  7:36831/104‑1 (MQ=255)

                                                                                                                                                                                                                                        tttGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCcgc                                                                                                                                                      >  5:1167215/1‑149 (MQ=255)

                                                                                                                                                                                                                                          tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCcgcgc                                                                                                                                                    <  6:781050/149‑1 (MQ=255)

                                                                                                                                                                                                                                          tGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATGGCACCACATAATGCTGCGACGTAGAAAGCAAAGAGTAGCCGAAATGAAGAAGGGAAAAGGAATAAAACCACCAGCGGCCGACataaagacagg                      >  4:72603/1‑267 (MQ=255)

                                                                                                                                                                                                                                          ggttagccAAAATCCACAAAAGACAGTGAAGCCCCTGGCCCGCGTCTTCTGAAGATTTTTCAAATCGCACAAGGCCGNGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATGCTGCCCGGGCAGCAAGCTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGG‑CGACATATAAACCAGt                     <  3:65476/268‑1 (MQ=255)

                                                                                                                                                                                                                                               tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGtt                                                                                                                                                                                                                                              >  8:888213/1‑54 (MQ=255)

                                                                                                                                                                                                                                               tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGtt                                                                                                                                                                                                                                              <  7:888213/54‑1 (MQ=255)

                                                                                                                                                                                                                                               tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAAt                                                                                                                                               <  5:1294217/149‑1 (MQ=255)

                                                                                                                                                                                                                                               tGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAAt                                                                                                                                               >  7:43483/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaaTGCAACAAATCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTACAAATCGCACAAGGCCAGGCGCCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAg                                                                                                                                            >  8:206797/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAg                                                                                                                                            >  7:94261/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAg                                                                                                                                            >  6:632627/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAg                                                                                                                                            >  7:421461/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAg                                                                                                                                            >  5:1245637/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAg                                                                                                                                            <  8:156830/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  aaaaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGAGCAATCAg                                                                                                                                            <  5:103061/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aaTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGa                                                                                                                                          >  8:284017/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGaa                                                                                                                                         >  5:1333299/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aTGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGaa                                                                                                                                         >  5:32058/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       gCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggc                                                                                                                                                                                                                         >  7:666137/1‑67 (MQ=255)

                                                                                                                                                                                                                                                       gCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAggc                                                                                                                                                                                                                         <  8:666137/67‑1 (MQ=255)

                                                                                                                                                                                                                                                              aGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAAt                                                                                                                                >  5:354088/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              aGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGG‑CGACATATAAACCAGTAATGCGGAGATAAGCAGGc  >  2:302941/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                ccAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATac                                                                                                                              <  5:28300/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      aaaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGa                                                                                                                        <  7:101478/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       aaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGtt                                                                                                                                                                                             <  8:1132120/79‑1 (MQ=255)

                                                                                                                                                                                                                                                                       aaTCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGtt                                                                                                                                                                                             >  7:1132120/1‑79 (MQ=255)

CGGTCAGACAGCGTCGATTTACCGTGGTCAATGTGAGCTATGATCGAAAAGTTACGTATATTCTTCATAAAGTATGATTATTGTGCCTTAATGCCCGGTTAACCAGGCTTTTAGAAGTCGCTGTTCTGAGCTTAACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCACTGGCTCGCGTCTTCCGAAGATGTTTCAAATCGCACAAGGCCAGGCGGCAATGCCACGCTTAAGATGATCGGTTGCCATTCTGCCCGGGCAGCAAACTTACGCGAGTAGCCGCGCGCAATCAGGAACCCGCCAATACCACCGAGGATCGCACCACATAATGCTGCGACGTCAGAAGCAAAGAGTAGCTGAAATAGCGAAGCGATAAGGAATAATCCCACCAGCGG‑CGACATATAAACCAGTAATGCGGAGCTAAGCAGGC  >  NC_000913/2707057‑2707584

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