BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,044,385	Δ1 bp	100%	coding (187/1431 nt)	yihF →	DUF945 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,044,385	0	T	.	100.0%	529.4 / NA	116	coding (187/1431 nt)	yihF	DUF945 family protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base . (52/64); total (52/64)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 120 total aligned reads displayed.

CCGTAATCTGGGGTGGTGGCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATGGTAATGTC  >  NC_000913/4044239‑4044528

ccGTAATCTGGGGTGGTGGCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttt                                                                                                                                               <  7:10644/149‑1 (MQ=255)

 cGTAATCTGGGGTGGTGGCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttt                                                                                                                                              <  1:103351/149‑1 (MQ=255)

 cGTAATCTGGGGTGGTGGCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttt                                                                                                                                              <  7:147570/149‑1 (MQ=255)

   tAATCTGGGGTGGTGGCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttttt                                                                                                                                             >  5:375200/1‑148 (MQ=255)

        tGGGGTGGTGGCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTcc                                                                                                                                      >  4:12029/1‑149 (MQ=255)

             tggtggCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACACGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGttt                                                                                                                                 <  2:32904/149‑1 (MQ=255)

                   cacaTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaat                                                                                                                           <  3:3418/149‑1 (MQ=255)

                   cacaTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaat                                                                                                                           <  6:220268/149‑1 (MQ=255)

                   cacaTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaat                                                                                                                           <  1:308298/149‑1 (MQ=255)

                    acaTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaatg                                                                                                                          <  1:216229/149‑1 (MQ=255)

                     caTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaatg                                                                                                                          >  7:151931/1‑148 (MQ=255)

                       tGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaatgcaa                                                                                                                       <  5:110861/149‑1 (MQ=255)

                        ggTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAATAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCACGTTTTcaaatgcaaa                                                                                                                      >  1:81780/1‑149 (MQ=255)

                        ggTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaatgcaaa                                                                                                                      >  1:303493/1‑149 (MQ=255)

                         ggacacAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaatgcaaa                                                                                                                      <  6:35130/146‑1 (MQ=255)

                          tACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTACCGATGCTAAAAAGAAAGGTGACCATGCATACGATATGACGTTATTTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaatgtaaat                                                                                                                     >  3:95523/1‑148 (MQ=255)

                          tACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaatgcaaat                                                                                                                     >  4:14516/1‑148 (MQ=255)

                          tACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTcaaa                                                                                                                            <  6:357613/141‑1 (MQ=255)

                               aGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACAtt                                                                                                               >  6:80082/1‑149 (MQ=255)

                                    cGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATa                                                                                                          <  7:307440/149‑1 (MQ=255)

                                       aaaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACg                                                                                                       >  1:138165/1‑149 (MQ=255)

                                        aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACgg                                                                                                      >  5:172319/1‑149 (MQ=255)

                                             aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa                                                                                                  <  3:12286/148‑1 (MQ=255)

                                             aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa                                                                                                  <  8:111876/148‑1 (MQ=255)

                                             aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAc                                                                                                 >  2:339586/1‑149 (MQ=255)

                                             aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAc                                                                                                 >  7:377163/1‑149 (MQ=255)

                                              gCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAcc                                                                                                <  8:377163/149‑1 (MQ=255)

                                              gCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAcc                                                                                                <  3:204952/149‑1 (MQ=255)

                                               ccTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAccc                                                                                               >  6:74566/1‑149 (MQ=255)

                                                   gtgtCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGAtc                                                                                           >  3:357737/1‑149 (MQ=255)

                                                      tCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGAtctc                                                                                         >  8:194846/1‑148 (MQ=255)

                                                        gAATAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACTTGCCTACGATTTGACGTTATGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCatt                                                                                      <  6:327905/149‑3 (MQ=255)

                                                        gAAAAATTTATTAAAGATTTTAACGATTCTAAAAAGAAAGGTGAACATGCCTACGATATGACTTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCaa                                                                                       <  4:78686/148‑1 (MQ=255)

                                                        gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCaa                                                                                       >  3:264482/1‑148 (MQ=255)

                                                        gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAat                                                                                      <  5:268574/149‑1 (MQ=255)

                                                         aaaaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAata                                                                                     >  2:212358/1‑149 (MQ=255)

                                                         aaaaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAata                                                                                     <  2:329019/149‑1 (MQ=255)

                                                             aTTTATTAAAGATTTTAACGATGCTAAAAAGACAGGTGAACCTGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑GTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCaa                                                                                 <  2:150792/149‑1 (MQ=255)

                                                                tattaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgc                                                                               <  3:12029/148‑1 (MQ=255)

                                                                  ttaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgtcag                                                                            >  8:24936/1‑149 (MQ=255)

                                                                  ttaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccag                                                                            <  1:291983/149‑1 (MQ=255)

                                                                     aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccaggc                                                                          >  7:343268/1‑148 (MQ=255)

                                                                     aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccaggc                                                                          >  7:86206/1‑148 (MQ=255)

                                                                     aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccaggc                                                                          >  1:382141/1‑148 (MQ=255)

                                                                     aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTCTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccaggc                                                                          >  5:34850/1‑148 (MQ=255)

                                                                      aGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccaggcca                                                                        <  1:356689/149‑1 (MQ=255)

                                                                       gATTTTAACGATGCTAAAAAGAAAGGGTAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATACCCGTTTTCAAATGTAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccaggccag                                                                       <  4:95523/149‑1 (MQ=255)

                                                                         ttttAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGTCAGGCCAGa                                                                      >  2:177838/1‑148 (MQ=255)

                                                                              aCGATGCTAAAAAGTAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGtt                                                                 <  2:131322/148‑1 (MQ=255)

                                                                              aCGATGCTAAAAAGATAGGGTAACATGCCTACGATATGACGTTAAGTTATCTAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGTTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGtt                                                                 <  8:371789/148‑1 (MQ=255)

                                                                              aCGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGtt                                                                 <  4:221312/148‑1 (MQ=255)

                                                                              aCGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGtt                                                                 <  8:187728/148‑1 (MQ=255)

                                                                              aCGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGtt                                                                 <  7:284017/148‑1 (MQ=255)

                                                                                 aTGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGTCAGGCCAGAAAGTTGTAt                                                             <  1:177838/149‑1 (MQ=255)

                                                                                  tGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAAATCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTAtt                                                            <  6:275346/149‑1 (MQ=255)

                                                                                     tAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGAAAAAGGC‑TTTTTAATTCCCGTTTTCAAATTCAAATGACTTTCGATAACGGTGCAACCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGa                                                         >  5:32588/1‑149 (MQ=255)

                                                                                      aaaaaGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGAt                                                        >  1:59652/1‑149 (MQ=255)

                                                                                       aaaaGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTAAATTACCGATTTCAAATGCAAATGACATTCGATAACTGTTCACACGATATCAATATCAAGACAGGACAGAAAGTTGTATTTGAtg                                                       >  1:198775/1‑149 (MQ=255)

                                                                                        aaagaaagGTGAACATGACTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGAtgt                                                      <  6:239732/149‑1 (MQ=255)

                                                                                         aagaaagGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGAtgtg                                                     >  4:281116/1‑149 (MQ=255)

                                                                                            aaagGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGAt                                                  <  2:157703/149‑1 (MQ=255)

                                                                                            aaagGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGAt                                                  <  5:49281/149‑1 (MQ=255)

                                                                                             aagGTGAACATGCATACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATg                                                 <  5:74566/149‑1 (MQ=255)

                                                                                              agGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCAAGAAAGTTGTATTTGATGTGGATGt                                                >  6:338623/1‑149 (MQ=255)

                                                                                                 tGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGa                                             <  5:4851/149‑1 (MQ=255)

                                                                                                  gAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAg                                            >  3:358910/1‑149 (MQ=255)

                                                                                                    aCATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCa                                          <  7:92336/149‑1 (MQ=255)

                                                                                                        gCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGt                                      >  1:151529/1‑149 (MQ=255)

                                                                                                             cGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGtt                                 <  7:124/149‑1 (MQ=255)

                                                                                                             cGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGt                                  <  5:80082/148‑1 (MQ=255)

                                                                                                                tattaagTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCAGTTTTCAAATGCAAATGAAATTCTATAACGGTGCACCCTATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGc                               <  2:81780/142‑1 (MQ=255)

                                                                                                                tatGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAat                                                                                      <  1:395832/93‑1 (MQ=255)

                                                                                                                tatGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAat                                                                                      >  2:395832/1‑93 (MQ=255)

                                                                                                                      gTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATca                         >  2:67655/1‑148 (MQ=255)

                                                                                                                           gTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGt                   <  4:6651/149‑1 (MQ=255)

                                                                                                                           gTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGt                   <  5:276061/149‑1 (MQ=255)

                                                                                                                            ttATCAATATTTTGACAATGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACTGTGCTCCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGtt                  <  6:375200/149‑1 (MQ=255)

                                                                                                                             tATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATACCGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGtt                  >  2:359146/1‑148 (MQ=255)

                                                                                                                              aTCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGCCATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGAACGGTCCGTTGCCCATCACAATGTaaa                >  5:178663/1‑149 (MQ=255)

                                                                                                                               tCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAAt               <  2:138165/149‑1 (MQ=255)

                                                                                                                               tCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAAt               <  6:19198/149‑1 (MQ=255)

                                                                                                                                 aaaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGaa                                                                     <  8:395861/93‑1 (MQ=255)

                                                                                                                                 aaaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGaa                                                                     >  7:395861/1‑93 (MQ=255)

                                                                                                                                 aaaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAAt               >  6:296779/1‑147 (MQ=255)

                                                                                                                                 aaaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAAt               <  5:296779/147‑1 (MQ=255)

                                                                                                                                  aaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCa            <  6:172319/149‑1 (MQ=255)

                                                                                                                                  aaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCa            <  2:151529/149‑1 (MQ=255)

                                                                                                                                  aaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCa            <  4:357737/149‑1 (MQ=255)

                                                                                                                                  aaaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACCGTGCACCCGATCTCAATATCAAGCCAGGCCTGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACATTGTTAATGCa            >  3:20571/1‑149 (MQ=255)

                                                                                                                                   aaTTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCAt           <  2:59652/149‑1 (MQ=255)

                                                                                                                                     ttttGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATgg         >  6:346196/1‑149 (MQ=255)

                                                                                                                                     ttttGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATgg         >  6:311892/1‑149 (MQ=255)

                                                                                                                                     ttttGACAAAGGC‑TTTTTAATTCCCCTTTTCCAATGCAAATGACATTCGTTAACGGTGCACCCGATCTCAAAAACCAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATgg         >  4:10779/1‑149 (MQ=255)

                                                                                                                                      tttGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATgg         >  4:268608/1‑148 (MQ=255)

                                                                                                                                         gACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATGGTaa      <  7:194846/148‑1 (MQ=255)

                                                                                                                                         gACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGAGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATGGTaa      <  8:151931/148‑1 (MQ=255)

                                                                                                                                          aCAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATGGTAAt     <  3:268608/148‑1 (MQ=255)

                                                                                                                                           ccaaGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAAAGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATTGTAATGt   <  5:122453/147‑1 (MQ=255)

                                                                                                                                           cAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATGGTAATGt   >  3:213294/1‑149 (MQ=255)

                                                                                                                                            aaaGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATGGTAATGTc  <  2:382141/149‑1 (MQ=255)

CCGTAATCTGGGGTGGTGGCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAATGCATGGTAATGTC  >  NC_000913/4044239‑4044528

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