BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	987,560	G→C	5.6%	intergenic (‑578/+25)	ompF ← / ← asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,560	0	G	C	5.6%	489.2 / 15.3	142	intergenic (‑578/+25)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand): ref base G (71/63); new base C (3/5); total (74/68)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

Only 100 of 145 total aligned reads displayed.

CTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTT  >  NC_000913/987332‑987796

cTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTTACAAAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGGTCTTGAAGACGCTAAATTAGAAGCTTGCGTTACGCGGa                                                                                                                                                                                                 >  2:273686/1‑275 (MQ=255)

 tGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGGCAGCGGGCGCTGGCTCTTAAAGACGCCAAATTATAAGCTGGCGTTAcgccgga                                                                                                                                                                                                >  2:224134/1‑271 (MQ=255)

                                                                           aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAATTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                                                                                                                 >  4:300896/1‑216 (MQ=255)

                                                                               aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcac                                                                                                                                                                             <  2:300896/216‑1 (MQ=255)

                                                                                tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                                                                                                                                                                               >  5:187850/1‑149 (MQ=255)

                                                                                tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                                                                                                                                                                               <  6:727595/149‑1 (MQ=255)

                                                                                   gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                                                                                                                                                                                            >  7:1452734/1‑149 (MQ=255)

                                                                                      aTAAGAATGCGGATTAGTTTTTGCGCGGAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTg                                                                                                                                                                                                                                         >  7:1293769/1‑149 (MQ=255)

                                                                                        aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGc                                                                                                                                                                                                                                       <  8:1022217/149‑1 (MQ=255)

                                                                                         aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGct                                                                                                                                                                                                                                      <  6:1168027/149‑1 (MQ=255)

                                                                                             tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTCAGTGGGAACGGAATCACATCACGTAAGCTTTTCACGCCAGGTACGTAAGCAATCAGTCGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGa                                                                                                   >  2:216208/1‑275 (MQ=255)

                                                                                             tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttt                                                                                                                                                                                                                                  <  8:195145/149‑2 (MQ=255)

                                                                                              gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                                                                                                                                                                 <  5:170740/149‑1 (MQ=255)

                                                                                                     aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACACCACGTACGTTTTTCACGCCCGTTACGTAACCACTCAGACGTTCAAAACCCAG‑ACCGAAAACTGAATGCGGACCAGTAc                                                                                            >  2:50370/1‑274 (MQ=255)

                                                                                                       tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaa                                                                                                                                                                                                                        <  6:907170/149‑1 (MQ=255)

                                                                                                               cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAg                                                                                                                                                                                                                <  5:396786/149‑1 (MQ=255)

                                                                                                                gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                                                                                                                                                                                                               >  8:354674/1‑149 (MQ=255)

                                                                                                                gTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAAAAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGCCGCCGGCTCTTGAAGACGCTTAAATAGAAGc                                                                                                                                                                                                               >  6:1093306/1‑149 (MQ=255)

                                                                                                                   aaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTgg                                                                                                                                                                                                            <  6:599201/149‑1 (MQ=255)

                                                                                                                         gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                                                                                                                                                                   <  5:600476/120‑1 (MQ=255)

                                                                                                                         gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                                                                                                                                                                   >  6:600476/1‑120 (MQ=255)

                                                                                                                         gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGctc                                                                                                                                                                                                                                     >  7:802367/1‑118 (MQ=255)

                                                                                                                               tAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                                                                                                                                                                >  8:1046438/1‑149 (MQ=255)

                                                                                                                                 atatTCCGTAACAGGATGATCGTTATGGCGGGGTTTCAAATAAATTTCTTGACTTAACAAACAGATAACTTGACACAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCCAAATTAGAAGCTGGCGTTACGCGGAGtg                                                                                                                                                                                              >  8:1108793/1‑148 (MQ=255)

                                                                                                                                 atatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTa                                                                                                                                                                                              >  5:221129/1‑149 (MQ=255)

                                                                                                                                     tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                                                                                                                                                          <  8:499445/149‑1 (MQ=255)

                                                                                                                                     tCCGTAAAAGGATGATCGTTAAGGCGGGGTTGAAAATAGATTGCTTGACTTAACAAACAGGTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCGTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAAt                                                           <  4:224134/275‑1 (MQ=255)

                                                                                                                                       cGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGATGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGAGGGAAGGGAATAACATCACGTACGTTTTGAACGCCAGTTACGTAAGCAATCAGACGTTCAAAACAACG‑ACCGAAACCTGAATGGGGAAAAGTACAGTAGCGAAGCAAATCGCGATACAACCAgacatc                                                          >  1:73056/1‑269 (MQ=255)

                                                                                                                                        gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTg                                                                                                                                                                                                             <  8:363508/127‑1 (MQ=255)

                                                                                                                                        gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGa                                                                                                                                                                                       >  8:428192/1‑149 (MQ=255)

                                                                                                                                         tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTTACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAACTAGCGTTACGCGGAGTACCTGGGAACGGAATCCCAACACATAAGTTTTGCCCCCCAGTTACGTAAGCCATCAGACGTTCAAAAACAAGCACCAAAACATAAATGCGAAAC‑CAACCCTAGCGCCACAGATCGCAATACCAACAGAAAacttc                                                       >  4:144788/1‑275 (MQ=255)

                                                                                                                                              ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                                                                                                                 >  6:986168/1‑149 (MQ=255)

                                                                                                                                              ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                                                                                                                 <  8:1371111/149‑1 (MQ=255)

                                                                                                                                                atgatCGTTATGGCGGAGATGCAAATAGATTGCCTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                                                                                                                                                               <  5:1144810/149‑1 (MQ=255)

                                                                                                                                                  gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcac                                                                                                                                                                             <  6:694074/149‑1 (MQ=255)

                                                                                                                                                   atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                                                                                                                                                            <  8:782133/149‑1 (MQ=255)

                                                                                                                                                    tCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                                                                                                                                                                           >  6:492149/1‑149 (MQ=255)

                                                                                                                                                      gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                                                                                                                                                                         >  7:1043633/1‑149 (MQ=255)

                                                                                                                                                         aTGGCGGGGTTGCAAATAGAATGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGACAGCGGGCGCTGGCGCTTGAAGACGCTAAATTAAAAGCTGGCGTTa                                                                                                                                                                                                       >  4:99463/1‑116 (MQ=255)

                                                                                                                                                           ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                                                                                                                     <  7:1362930/132‑1 (MQ=255)

                                                                                                                                                             cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCGAAATTGGAAGCTGGCGTTAcgcg                                                                                                                                                                                                   <  2:99463/116‑1 (MQ=255)

                                                                                                                                                               gggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                                                                                                                                                                    <  7:263249/145‑1 (MQ=255)

                                                                                                                                                               gggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                                                                                                                                                                <  5:986168/149‑1 (MQ=255)

                                                                                                                                                                ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt                                                                                                                                                               <  7:428192/149‑1 (MQ=255)

                                                                                                                                                                ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtttt                                                                                                                                                               >  7:936428/1‑149 (MQ=255)

                                                                                                                                                                  ttGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                                                                                                             >  7:1071709/1‑149 (MQ=255)

                                                                                                                                                                  ttGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                                                                                                             >  8:726201/1‑149 (MQ=255)

                                                                                                                                                                      aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTGAATTa                                                                                                                                                                                                                    <  5:960250/90‑1 (MQ=255)

                                                                                                                                                                      aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTa                                                                                                                                                                                                                    >  6:960250/1‑90 (MQ=255)

                                                                                                                                                                      aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                                                                                                                         >  6:237437/1‑149 (MQ=255)

                                                                                                                                                                      aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                                                                                                                         >  5:1235245/1‑149 (MQ=255)

                                                                                                                                                                      aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                                                                                                                         <  7:354674/149‑1 (MQ=255)

                                                                                                                                                                           gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt                                                                                                                                                    <  6:1235245/149‑1 (MQ=255)

                                                                                                                                                                                cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                                                                                                                                                                      >  8:1370116/1‑110 (MQ=255)

                                                                                                                                                                                cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                                                                                                                                                                      <  7:1370116/110‑1 (MQ=255)

                                                                                                                                                                                cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                                                                                                                                                          <  2:377412/138‑1 (MQ=255)

                                                                                                                                                                                cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa                                                                                                                                               >  8:943866/1‑149 (MQ=255)

                                                                                                                                                                                cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa                                                                                                                                               >  7:505188/1‑149 (MQ=255)

                                                                                                                                                                                   gACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGc                                                                                                                                            >  7:699973/1‑149 (MQ=255)

                                                                                                                                                                                    aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                                                                                                                 >  5:595698/1‑111 (MQ=255)

                                                                                                                                                                                    aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCa                                                                                                                                           >  6:25296/1‑149 (MQ=255)

                                                                                                                                                                                      ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                                                                                                                         >  5:1151999/1‑149 (MQ=255)

                                                                                                                                                                                      ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                                                                                                                         >  5:41941/1‑149 (MQ=255)

                                                                                                                                                                                        aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAGTACCGTAGCGa                                                                                    >  1:455/1‑199 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatcac                                                                                                                                                                        <  8:810102/113‑1 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                                                                                                                >  7:1102289/1‑105 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                                                                                                                <  8:1102289/105‑1 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                                                                                                                    >  7:37869/1‑149 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                                                                                                                    <  8:936428/149‑1 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                                                                                                                    <  8:572339/149‑1 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                                                                                                                    <  7:547058/149‑1 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                                                                                                                    >  7:634393/1‑149 (MQ=255)

                                                                                                                                                                                           aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                                                                                                                    >  6:1416979/1‑149 (MQ=255)

                                                                                                                                                                                            aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCa                                                                                <  3:455/199‑1 (MQ=255)

                                                                                                                                                                                             acaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGt                                                                                                                                  <  8:634393/149‑1 (MQ=255)

                                                                                                                                                                                               aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTc                                                                                                                                >  6:1144573/1‑149 (MQ=255)

                                                                                                                                                                                               aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTc                                                                                                                                >  5:790703/1‑149 (MQ=255)

                                                                                                                                                                                               aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGttctt  <  2:144788/274‑1 (MQ=255)

                                                                                                                                                                                                gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCa                                                                                                                               <  8:1452734/149‑1 (MQ=255)

                                                                                                                                                                                                 aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                                                                                                                                                                                 >  6:1430038/1‑66 (MQ=255)

                                                                                                                                                                                                 aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                                                                                                                                                                                 <  5:1430038/66‑1 (MQ=255)

                                                                                                                                                                                                 aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCaa                                                                                                                              <  5:1093306/149‑1 (MQ=255)

                                                                                                                                                                                                      ttGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCa                                                                                                                         >  7:880379/1‑149 (MQ=255)

                                                                                                                                                                                                       tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCaa                                                                                                                        <  7:356365/149‑1 (MQ=255)

                                                                                                                                                                                                         aCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑a                                                                                                                     <  7:1046438/149‑1 (MQ=255)

                                                                                                                                                                                                          cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑Ac                                                                                                                    <  7:943866/149‑1 (MQ=255)

                                                                                                                                                                                                           aGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑a                                                                                                                     >  7:838628/1‑147 (MQ=255)

                                                                                                                                                                                                           aGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑a                                                                                                                     <  8:838628/147‑1 (MQ=255)

                                                                                                                                                                                                            gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc                                                                                                                                                  <  6:668331/118‑1 (MQ=255)

                                                                                                                                                                                                            gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc                                                                                                                                                  >  5:668331/1‑118 (MQ=255)

                                                                                                                                                                                                             aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                                                                                                                                                                                    <  8:1234693/67‑1 (MQ=255)

                                                                                                                                                                                                             aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                                                                                                                                                                                    >  7:1234693/1‑67 (MQ=255)

                                                                                                                                                                                                                   tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTg                                                                                                           >  5:1353954/1‑149 (MQ=255)

                                                                                                                                                                                                                      aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAAt                                                                                                        <  5:223222/149‑1 (MQ=255)

                                                                                                                                                                                                                      aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAAt                                                                                                        >  7:461002/1‑149 (MQ=255)

                                                                                                                                                                                                                         aaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCg                                                                                                     <  8:437458/149‑1 (MQ=255)

                                                                                                                                                                                                                           gagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGa                                                                                                   >  7:144523/1‑149 (MQ=255)

                                                                                                                                                                                                                            agCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGaa                                                                                                  <  6:790703/149‑1 (MQ=255)

                                                                                                                                                                                                                               cAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAg                                                                                               >  5:798344/1‑149 (MQ=255)

                                                                                                                                                                                                                                    ggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAGTACCg                                                                                          >  7:646748/1‑149 (MQ=255)

CTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAG‑ACCGAAACCTGAATGCGGAACAGTACCGTAGCGACGCAGATCGCGATACCACCAGTAATCTTCTTTATTCAGGCCCATTTCCAGCATACGCTCGTCCAGCACGTCCAGACGTTCTT  >  NC_000913/987332‑987796

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