BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,560,317	A→C	5.6%	N21H (AAT→CAT)	yffL →	CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,560,317	0	A	C	5.6%	513.1 / 18.7	143	N21H (AAT→CAT)	yffL	CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (77/58); new base C (5/3); total (82/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

Only 100 of 146 total aligned reads displayed.

TTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCCAGATAATCACGAGGAAACAGATCATCACGAGGCAAAACAGCAGGAGCAAGCCTGGCGTGGCAACC‑TGCGCTACCTGGACACATTAAACCGAC  >  NC_000913/2560117‑2560587

tttgttttgtGTGTTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTGTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTGTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCACGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                                                                                                                                                                                   <  1:144193/279‑1 (MQ=255)

tttgttttgtGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                                                                                                                                                                                                       >  4:156735/1‑275 (MQ=255)

        tntgtgTTTAACTATCAACTTATAAAAAAACTTNNTTCAATTTGCCGGATTCGGCCGTCTTTTCTNTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCATGAACTTAAATCCAAGGTTGCNACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca                                                                                                                                                                                           <  3:90290/274‑1 (MQ=255)

        gtgtgGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCGGCAACATGCCGGGGGCGGTATGAGAAACAAGCAAAAGGAGAATTTCCCAGGTCGGATATGCTAGCAGTAGAAAAAGa                                                                                                                                                                                               >  2:646650/1‑275 (MQ=255)

                            gATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGGAGAACAAGAAACAAACAATGATGTTCGTc                                                                                                                                                                           >  2:452411/1‑275 (MQ=255)

                                                     cGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGaa                                                                                                                                                                                                                                                                                >  5:478238/1‑149 (MQ=255)

                                                     cGGATTCGTCCGGCTTTTCTTTTNCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                                                                                                                                                                                                   <  3:158083/162‑1 (MQ=255)

                                                       gATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAt                                                                                                                                                                                                                                                                               >  7:241901/1‑148 (MQ=255)

                                                              ccGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                                                                                                                                                                                                       <  5:541138/149‑1 (MQ=255)

                                                              ccGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                                                                                                                                                                                                       <  5:339884/149‑1 (MQ=255)

                                                                ggCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                                                                                                                                                                                                      <  7:682668/148‑1 (MQ=255)

                                                                      tcttttTCGCTCATTATTCCATTATTCCAATTAAGTAGAATAAATTATTAAATATCAAGAACTTAAATCCATGTTTGCCACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgg                                                                                                                                                                                                                                                                <  7:342155/148‑1 (MQ=255)

                                                                      tcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                                                                                                                                                                                               <  5:724313/149‑1 (MQ=255)

                                                                          tttCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                                                                                                                                                                                           <  6:824339/149‑1 (MQ=255)

                                                                                tAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                                                                                                                                                                                                     >  5:240008/1‑149 (MQ=255)

                                                                                    tattccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                                                                                                                                                                                                  <  7:401892/148‑1 (MQ=255)

                                                                                      ttccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGAGGTATGAGAAACa                                                                                                                                                                                                                                               >  6:758908/1‑149 (MQ=255)

                                                                                           ttattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                                                                                                                                                                                           >  7:24333/1‑148 (MQ=255)

                                                                                             attccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                                                                                                                                                                                        <  5:248603/149‑1 (MQ=255)

                                                                                             attccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaa                                                                                                                                                                                                                                         <  8:212466/148‑1 (MQ=255)

                                                                                              ttccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAAAGACAAAACACgaaa                                                                                                         >  4:399291/1‑274 (MQ=255)

                                                                                               tccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGctttta                                                                                                                                                                                                                                                       <  8:731783/132‑6 (MQ=38)

                                                                                               tccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                                                                                                                                                                                       >  7:312464/1‑148 (MQ=255)

                                                                                                 caATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                                                                                                                                                                                                    >  6:133288/1‑149 (MQ=255)

                                                                                                  aaTTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                                                                                                                                                                                            <  8:651756/124‑1 (MQ=255)

                                                                                                   aTTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                                                                                                                                                                                                   <  7:209187/148‑1 (MQ=255)

                                                                                                      aaGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACg                                                                                                 <  2:156735/275‑1 (MQ=255)

                                                                                                         tGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                                                                                                                                                                                           <  6:924100/134‑1 (MQ=255)

                                                                                                          ggAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                                                                                                                                                                                           <  6:997103/149‑1 (MQ=255)

                                                                                                           gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCCAGa                                                                                            >  4:439315/1‑275 (MQ=255)

                                                                                                            aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                                                                                                                                                                                          >  7:712682/1‑148 (MQ=255)

                                                                                                             aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCCAGATa                                                                                          >  4:557708/1‑275 (MQ=255)

                                                                                                                 attattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGa                                                                                                                                                                                                                    <  6:53138/149‑1 (MQ=255)

                                                                                                                 attattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGa                                                                                                                                                                                                                    >  5:499683/1‑149 (MQ=255)

                                                                                                                     ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGANNATTTCTCan                                                                                                                                                                                                                          >  3:40306/1‑138 (MQ=255)

                                                                                                                     ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAGTTCTCAt                                                                                                                                                                                                                          <  1:40306/139‑1 (MQ=255)

                                                                                                                      tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCGCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                                                                                                                                                                                                <  7:417719/148‑1 (MQ=255)

                                                                                                                      tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                                                                                                                                                                                                  >  3:8409/1‑130 (MQ=255)

                                                                                                                      tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                                                                                                                                                                                                <  8:312464/148‑1 (MQ=255)

                                                                                                                      tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                                                                                                                                                                                               <  6:125596/149‑1 (MQ=255)

                                                                                                                         atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                                                                                                                                                                                             >  7:783472/1‑148 (MQ=255)

                                                                                                                          tatCAAGAACTTAAATCCATGTTTGCAACGACAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTACAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAAAAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGCGAACAAACACGAAGAACACCCGCAAGAAAATACAGAGGaaaa                                                                             >  4:310171/1‑274 (MQ=255)

                                                                                                                             cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                                                                                                                                                                                                         <  7:968596/148‑1 (MQ=255)

                                                                                                                              aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                                                                                                                                                                                        <  8:542732/148‑1 (MQ=255)

                                                                                                                                    ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                                                                                                                                                                                                 >  6:37986/1‑149 (MQ=255)

                                                                                                                                    ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                                                                                                                                                                                                  <  7:461859/148‑1 (MQ=255)

                                                                                                                                    ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGACTAAAGCACCCAAAAACGAAGAGACTAAACACAAAGAACACACGCCAGATAATCACGAGGAAACAGATCATCACGAgg                                                               >  1:57688/1‑279 (MQ=255)

                                                                                                                                      aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                                                                                                                                                                                                <  7:86175/148‑1 (MQ=255)

                                                                                                                                      aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCAAGATAATCACGAGGAAACAGATCATCACGa                                                                 >  2:349320/1‑275 (MQ=255)

                                                                                                                                        aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                                                                                                                                                                                                           >  8:343763/1‑71 (MQ=255)

                                                                                                                                         tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                                                                                                                                                                                              <  8:409696/99‑1 (MQ=255)

                                                                                                                                         tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                                                                                                                                                                                            >  5:300800/1‑149 (MQ=255)

                                                                                                                                          ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                                                                                                                                                                                            >  7:366303/1‑148 (MQ=255)

                                                                                                                                            aTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaa                                                                                                                                                                                         <  5:600095/149‑1 (MQ=255)

                                                                                                                                              gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTACAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAACTTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTACACACGAAGAACACACGCCAGATAATCACGAGGAAACAGATCATCACGAGGCaaaa                                                          >  2:608508/1‑274 (MQ=255)

                                                                                                                                              gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACa                                                                                                                                                               >  1:105261/1‑173 (MQ=255)

                                                                                                                                              gTTTGCAACGAAAGAGCCAGAGTTCGAAAATAGGNACAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACa                                                                                                                                                               <  3:105261/173‑1 (MQ=255)

                                                                                                                                                 tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTNGCAGTAGAACAGGAAACAAACAANGATGNTCGTCAATTTTTAAACAAAATTAACGa                                                                                                                                                     >  4:729396/1‑180 (MQ=255)

                                                                                                                                                     aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg                                                                                                                                                                                 >  7:578459/1‑148 (MQ=255)

                                                                                                                                                      cGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGtt                                                                                                                                                                               <  6:1054412/149‑1 (MQ=255)

                                                                                                                                                      cGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGt                                                                                                                                                                                >  8:127030/1‑148 (MQ=255)

                                                                                                                                                       gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGtt                                                                                                                                                                               >  8:819044/1‑148 (MQ=255)

                                                                                                                                                       gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGAGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTc                                                                                                                                                                              <  5:758908/149‑1 (MQ=255)

                                                                                                                                                        aaaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTc                                                                                                                                                                              >  8:90523/1‑148 (MQ=255)

                                                                                                                                                             tCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                                                                                                                                                                                             <  5:863593/112‑1 (MQ=255)

                                                                                                                                                             tCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                                                                                                                                                                                             >  6:863593/1‑112 (MQ=255)

                                                                                                                                                                agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtttt                                                                                                                                                                     >  5:191629/1‑149 (MQ=255)

                                                                                                                                                                 gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtttt                                                                                                                                                                     >  8:743777/1‑148 (MQ=255)

                                                                                                                                                                    ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATTTCTGATATGCTATCAGTAGAACAAGANACAAACANTGANGTTCGTNCATTTTTaaa                                                                                                                                                                 >  3:32989/1‑149 (MQ=255)

                                                                                                                                                                      cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAACGCACCCAAAAACAAAGAGACTAAACACGAAGAACACACACCAGATAATCACGAGGAAACAGATCAACACGAGACAAAACAGCAGAAGCAAGACTGGCgcgaca                                 >  4:303423/1‑270 (MQ=255)

                                                                                                                                                                       gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACaa                                                                                                                                                              >  5:39908/1‑149 (MQ=255)

                                                                                                                                                                        aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtc                                                                                                                                                                                                                              >  8:63393/1‑84 (MQ=255)

                                                                                                                                                                        aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtc                                                                                                                                                                                                                              <  7:63393/84‑1 (MQ=255)

                                                                                                                                                                          aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                                                                                                                                                                                             <  2:780774/67‑1 (MQ=255)

                                                                                                                                                                               gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                                                                                                                                                                                               >  6:433906/1‑92 (MQ=255)

                                                                                                                                                                                aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGa                                                                                                                                                     >  5:149841/1‑149 (MQ=255)

                                                                                                                                                                                 tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                                                                                                                                                                                                         >  5:1078039/1‑96 (MQ=255)

                                                                                                                                                                                 tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                                                                                                                                                                                                         <  6:1078039/96‑1 (MQ=255)

                                                                                                                                                                                  cAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAg                                                                                                                                                    >  8:35358/1‑148 (MQ=255)

                                                                                                                                                                                   aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGt                                                                                                                                                   >  7:742116/1‑148 (MQ=255)

                                                                                                                                                                                   aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGt                                                                                                                                                   <  8:578459/148‑1 (MQ=255)

                                                                                                                                                                                     cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                                                                                                                                                                                        >  7:603264/1‑61 (MQ=255)

                                                                                                                                                                                      acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                                                                                                                                                                                             <  6:343732/71‑1 (MQ=255)

                                                                                                                                                                                      acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACg                                                                                                                                               <  5:37986/149‑1 (MQ=255)

                                                                                                                                                                                       caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACgg                                                                                                                                              >  5:586495/1‑149 (MQ=255)

                                                                                                                                                                                       caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACg                                                                                                                                               >  7:896847/1‑148 (MQ=255)

                                                                                                                                                                                        aaaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGa                                                                                                                                             >  6:2613/1‑149 (MQ=255)

                                                                                                                                                                                         aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtc                                                                                                                                                                                                                              >  6:108198/1‑67 (MQ=255)

                                                                                                                                                                                         aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtc                                                                                                                                                                                                                              <  5:108198/67‑1 (MQ=255)

                                                                                                                                                                                          aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCCAGATAATCACGAGGAAACAGACAATCACGGGGCAAAACACAAAGAGCAATCCTGGCGTGGCAACCATGC‑CTACCTGGAaaca            >  4:668602/1‑275 (MQ=255)

                                                                                                                                                                                              aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCCAGATAATCACGAGGAAACAGATCATCACGAGGCAAAACAGCAGGAGCAAGCCTGGCGTGGCAACC‑TGCGCTACCTGGACACATTAAACCg    >  1:95094/1‑279 (MQ=255)

                                                                                                                                                                                                 cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCa                                                                                                                                     >  8:364384/1‑148 (MQ=255)

                                                                                                                                                                                                  ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCAcc                                                                                                                                   >  6:1041920/1‑149 (MQ=255)

                                                                                                                                                                                                   cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACaaaa                                                                                                                                                            >  6:743736/1‑123 (MQ=255)

                                                                                                                                                                                                   cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCAccc                                                                                                                                  >  6:172373/1‑149 (MQ=255)

                                                                                                                                                                                                   cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCAccc                                                                                                                                  >  5:58833/1‑149 (MQ=255)

                                                                                                                                                                                                    ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGaa                                                                                                                         <  1:208597/157‑1 (MQ=255)

                                                                                                                                                                                                     cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCCAGATAATCACGAGGAAACAGATCATCACGAGGCAAAACAGCAGAGACAAGCCTGGCGTGGCAACC‑TGCCCAACCTGGACACATTAAAccagc  >  4:252014/1‑271 (MQ=255)

                                                                                                                                                                                                        aaTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCaaaaa                                                                                                                             <  6:24952/149‑1 (MQ=255)

                                                                                                                                                                                                        aaTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCaaaa                                                                                                                              <  8:241901/148‑1 (MQ=255)

TTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGCACCCAAAAACGAAGAGACTAAACACGAAGAACACACGCCAGATAATCACGAGGAAACAGATCATCACGAGGCAAAACAGCAGGAGCAAGCCTGGCGTGGCAACC‑TGCGCTACCTGGACACATTAAACCGAC  >  NC_000913/2560117‑2560587

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