BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,015,912	T→A	7.8%	intergenic (+30/+40)	rmf → / ← fabA	ribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,015,912	0	T	A	7.8%	1051.4 / 56.6	305	intergenic (+30/+40)	rmf/fabA	ribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase
Reads supporting (aligned to +/- strand): ref base T (148/132); new base A (15/9); total (163/142)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-02

Only 100 of 306 total aligned reads displayed.

ATTGCACATTTAGTAATCACTGTTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTCACTTTTTTCGCTGTCGGCAGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCC  >  NC_000913/1015653‑1016181

aTTGCACATTTAGTAATCACTGTTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCCCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGca                                                                                                                                                                                                                                                                >  2:5344/1‑274 (MQ=255)

                                     gAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCCAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCCTGGCGGACAGGGTGGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGCCTGATTATAACCAGATGGCAAAc                                                                                                                                                                                                                           >  4:693155/1‑275 (MQ=255)

                                                                                    gCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGt                                                                                                                                                                                                                         >  4:311572/1‑230 (MQ=255)

                                                                                                                 cAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAct                                                                                                                                                                                                                                                      <  4:266350/172‑1 (MQ=255)

                                                                                                                    gCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                                                                                                                                                                                          >  7:85708/1‑149 (MQ=255)

                                                                                                                      cGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAg                                                                                                                                                                                                                                                                        <  6:749313/149‑1 (MQ=255)

                                                                                                                        gCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGt                                                                                                                                                                                                                                                                      <  8:249901/149‑1 (MQ=255)

                                                                                                                         cATCGCCGGACGCCCAAAAGAACTGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                                                                                                                                                >  1:129636/1‑202 (MQ=255)

                                                                                                                          aTCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGttt                                                                                                                                                                                                                                                                    >  5:1289378/1‑149 (MQ=255)

                                                                                                                                ggACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCtt                                                                                                                                                                                                                                                              >  6:1442516/1‑149 (MQ=255)

                                                                                                                                    gCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTCTCGCCTTTTGa                                                                                                                                                                                                                                                          <  7:1092007/149‑1 (MQ=255)

                                                                                                                                      tCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCtttt                                                                                                                                                                                                                                                            <  6:706560/145‑1 (MQ=255)

                                                                                                                                      tCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATa                                                                                                                                                                                                                                                        >  6:958083/1‑149 (MQ=255)

                                                                                                                                      tCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGt                                                                                                                                                                                                                                                                      <  7:96486/135‑1 (MQ=255)

                                                                                                                                       aaaaaGAAATGTGTCCCTATCAGACGCTGAATCAAATGTCACAATGGCGGGGAGGCTGGCGAGAAGCCATGGCGNACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCg                                                                                                                                                                                       <  3:237555/212‑1 (MQ=255)

                                                                                                                                        aaaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGa                                                                                                                                                                                                                                                          <  7:1542725/145‑1 (MQ=255)

                                                                                                                                        aaaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGa                                                                                                                                                                                                                                                          >  8:1542725/1‑145 (MQ=255)

                                                                                                                                        aaaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAct                                                                                                                                                                                                                                                      >  8:1125669/1‑149 (MQ=255)

                                                                                                                                         aaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATActc                                                                                                                                                                                                                                                     >  5:1337498/1‑149 (MQ=255)

                                                                                                                                              aaTGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtct                                                                                                                                                                                                                                                <  5:54903/149‑1 (MQ=255)

                                                                                                                                                tgtgTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGa                                                                                                                                                                                                                                              <  6:8552/149‑1 (MQ=255)

                                                                                                                                                tgtgTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCg                                                                                                                                                                                             >  4:360721/1‑198 (MQ=255)

                                                                                                                                                  tgtCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATa                                                                                                                                                                                                                                                        >  1:52635/1‑137 (MQ=255)

                                                                                                                                                    tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCg                                                                                                                                                                                                                                                                           <  3:51479/116‑1 (MQ=255)

                                                                                                                                                    tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATtat                                                                                                                                                                                                                                          >  6:606809/1‑149 (MQ=255)

                                                                                                                                                    tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTgg                                                                                                                                                                                         <  2:360725/198‑1 (MQ=255)

                                                                                                                                                    tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTaa                                                                                                        >  1:264441/1‑279 (MQ=255)

                                                                                                                                                      ccTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATaa                                                                                                                                                                                                                                        <  8:748679/149‑1 (MQ=255)

                                                                                                                                                       cTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAAt                                                                                                                                                                                                                                       >  5:1622553/1‑149 (MQ=255)

                                                                                                                                                        tATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGcc                                                                                                                                                                                                                                                                <  2:765694/123‑1 (MQ=255)

                                                                                                                                                         aTCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTaaa                                                                                                       <  2:693160/275‑1 (MQ=255)

                                                                                                                                                                  tGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGACTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGa                                                                                                                                                                                                                            >  7:1426935/1‑149 (MQ=255)

                                                                                                                                                                   gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGAc                                                                                                                                                                                                                           <  5:958083/149‑1 (MQ=255)

                                                                                                                                                                   gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGAc                                                                                                                                                                                                                           <  5:183262/149‑1 (MQ=255)

                                                                                                                                                                   gAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGAc                                                                                                                                                                                                                           >  8:919415/1‑149 (MQ=255)

                                                                                                                                                                     aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgt                                                                                                                                                                                                                                                  <  7:648079/124‑1 (MQ=255)

                                                                                                                                                                        aaaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtc                                                                                                                                                                                                                                                 <  7:1281684/122‑1 (MQ=255)

                                                                                                                                                                        aaaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtc                                                                                                                                                                                                                                                 >  8:1281684/1‑122 (MQ=255)

                                                                                                                                                                        aaaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATc                                                                                                                                                                                                                      <  7:1611927/149‑1 (MQ=255)

                                                                                                                                                                         aaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCttt                                                                                                                                                                                                                                                             >  8:50338/1‑109 (MQ=255)

                                                                                                                                                                         aaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATcc                                                                                                                                                                                                                     <  5:466103/149‑1 (MQ=255)

                                                                                                                                                                         aaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATcc                                                                                                                                                                                                                     >  8:551123/1‑149 (MQ=255)

                                                                                                                                                                          aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGAt                                                                                                                                                                                                                                             <  6:314189/124‑1 (MQ=255)

                                                                                                                                                                          aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGAt                                                                                                                                                                                                                                             >  5:314189/1‑124 (MQ=255)

                                                                                                                                                                            gTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTgg                                                                                                                                                                                                                  >  5:1035139/1‑149 (MQ=255)

                                                                                                                                                                              cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACg                                                                                                                                                                                                                          <  4:739572/139‑1 (MQ=255)

                                                                                                                                                                              cacaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCa                                                                                                                                >  2:315020/1‑229 (MQ=255)

                                                                                                                                                                                 aaTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGaa                                                                                                                                                                                                                                  >  8:799479/1‑128 (MQ=255)

                                                                                                                                                                                  aTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                                                                                                                                            <  7:135246/149‑1 (MQ=255)

                                                                                                                                                                                  aTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                                                                                                                                            >  7:1109340/1‑149 (MQ=255)

                                                                                                                                                                                   tGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGAc                                                                                                                                                                                                           <  5:1206914/149‑1 (MQ=255)

                                                                                                                                                                                              cTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCttt                                                                                                                                                                                                                                                             >  6:178899/1‑88 (MQ=255)

                                                                                                                                                                                               tGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCCGGTCGCTGGc                                                                                                                                                                                        >  4:132323/1‑156 (MQ=255)

                                                                                                                                                                                                ggCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTc                                                                                                                                                                                              >  6:1160302/1‑149 (MQ=255)

                                                                                                                                                                                                  cGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGc                                                                                                                                                                                            <  6:455320/149‑1 (MQ=255)

                                                                                                                                                                                                   gagaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCa                                                                                                                                            <  4:49936/196‑1 (MQ=255)

                                                                                                                                                                                                   gagaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGACACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGta                                                                                                                                                                                    <  2:132325/156‑1 (MQ=255)

                                                                                                                                                                                                              gCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGa                                                                                                                                                                                                                                                                         >  8:603913/1‑60 (MQ=255)

                                                                                                                                                                                                              gCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGa                                                                                                                                                                                                                                                                         <  7:603913/60‑1 (MQ=255)

                                                                                                                                                                                                                       gtagtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAgacgac                                                                                                                                                                       <  5:1087324/149‑1 (MQ=255)

                                                                                                                                                                                                                        tagtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGAcc                                                                                                                                                                      <  5:70617/149‑1 (MQ=255)

                                                                                                                                                                                                                              tGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAAc                                                                                                                                                                <  5:638876/149‑1 (MQ=255)

                                                                                                                                                                                                                                gCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGttt                                                                                                                                                                                                                                                                    >  4:234139/1‑47 (MQ=255)

                                                                                                                                                                                                                                gCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTGTCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAGTCAGACGACGGTTAACAATACGTTTAAAATGAATACGGGTGGTCACTTTTTTCGCTGTCGGCCGTCCCTGACCAGTGACTTTCCCTTCGCCAACGccc                                 >  4:648483/1‑274 (MQ=255)

                                                                                                                                                                                                                                gCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTCACTTTTTTCGCTGTCGGCAGTACCTGACCAGTGAATTTCACTTCGCCAACGccc                                 >  2:434518/1‑274 (MQ=255)

                                                                                                                                                                                                                                    gATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGcc                                                                                                                                                                                                                                                                <  2:234142/47‑1 (MQ=255)

                                                                                                                                                                                                                                    gATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTCACTTTTTTCGCt                                                                           <  2:281607/228‑1 (MQ=255)

                                                                                                                                                                                                                                      ttCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACtt                                                                                                                                                        >  7:795472/1‑149 (MQ=255)

                                                                                                                                                                                                                                       tCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTc                                                                                                                                                       <  5:1426088/149‑1 (MQ=255)

                                                                                                                                                                                                                                       tCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTc                                                                                                                                                       <  8:1627286/149‑1 (MQ=255)

                                                                                                                                                                                                                                       tCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTc                                                                                                                                                       <  5:81047/149‑1 (MQ=255)

                                                                                                                                                                                                                                         tGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACaa                                                                                                                 <  2:673844/185‑1 (MQ=255)

                                                                                                                                                                                                                                          gTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGcc                                                                                                                                                    >  7:1232966/1‑149 (MQ=255)

                                                                                                                                                                                                                                           tctcTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCa                                                                                                                                                   >  6:1615173/1‑149 (MQ=255)

                                                                                                                                                                                                                                           tctcTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAAAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCa                                                                                                                                                   >  8:115791/1‑149 (MQ=255)

                                                                                                                                                                                                                                             tctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATc                                                                                                                                                 >  5:90224/1‑149 (MQ=255)

                                                                                                                                                                                                                                              ctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATcc                                                                                                                                                <  6:1289378/149‑1 (MQ=255)

                                                                                                                                                                                                                                               tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGGTTATAATCAGAAGGCAGACGTATCCTGGAACAGAc                                                                                                                                                                                                           <  8:666359/89‑1 (MQ=255)

                                                                                                                                                                                                                                               tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCg                                                                                                                                               >  6:1410540/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTACCTTTTTTCGCTTTCGGAACTACCTGACCAGTGAATTTCACTTCCCCAACGCCCACCCCGCCGCCTTTa                  >  4:418851/1‑274 (MQ=255)

                                                                                                                                                                                                                                                 tAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGcc                                                                                                                                             >  6:910287/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGcc                                                                                                                                             >  6:1677843/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCa                                                                                                                                            >  5:1360540/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                                                                                                                                          >  7:448959/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAgg                                                                                                                                          >  5:1641770/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaaGAAACCTCCGCATTGCGGAGGTTTCGCCGTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTgg                                                                                                                                                                                         <  2:627434/102‑1 (MQ=255)

                                                                                                                                                                                                                                                     aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAgac                                                                                                                                                                          >  7:40020/1‑116 (MQ=255)

                                                                                                                                                                                                                                                     aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAgac                                                                                                                                                                          <  8:40020/116‑1 (MQ=255)

                                                                                                                                                                                                                                                     aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGc                                                                                                                                         >  5:1001633/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGc                                                                                                                                         >  8:141525/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGc                                                                                                                                         >  5:381549/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       gAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGccc                                                                                                                                       >  7:351155/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        aaaCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGc                                                                                                                                         <  6:1066717/146‑1 (MQ=255)

                                                                                                                                                                                                                                                         aaCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCAt                                                                                                                                     >  8:1373708/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           ccTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATaa                                                                                                                                   >  7:304409/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            cTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAAt                                                                                                                                  >  7:726797/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCa                                                                                                                                      <  7:526828/144‑1 (MQ=255)

                                                                                                                                                                                                                                                             tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATc                                                                                                                                 >  6:1175034/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              ccGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCa                                                                                                                                <  6:1360540/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              ccGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCa                                                                                                                                <  5:1706600/149‑1 (MQ=255)

ATTGCACATTTAGTAATCACTGTTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGACGGTTAACAATGCGTTTAAAGTGAATACGGTAGGTCACTTTTTTCGCTGTCGGCAGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCC  >  NC_000913/1015653‑1016181

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