BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,106,959	A→T	5.1%	Y90N (TAT→AAT)	wbbH ←	O‑antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,106,959	0	A	T	5.1%	505.0 / 13.1	176	Y90N (TAT→AAT)	wbbH	O‑antigen polymerase
Reads supporting (aligned to +/- strand): ref base A (91/76); new base T (6/3); total (97/79)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

Only 100 of 193 total aligned reads displayed.

TTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCATAT  >  NC_000913/2106813‑2107107

ttCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTACCCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTTTGCACTTTCTTCGATGGTATGCTATaaa                                                                                                                                                    >  1:1063974/1‑149 (MQ=255)

ttCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATaaa                                                                                                                                                    >  1:536381/1‑149 (MQ=255)

ttCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATaaa                                                                                                                                                    >  1:925957/1‑149 (MQ=255)

  cAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATa                                                                                                                                                  >  2:1418531/1‑149 (MQ=255)

   aaCATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAg                                                                                                                                                 >  2:157133/1‑149 (MQ=255)

        cAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATTTATAGCatac                                                                                                                                            <  1:1100014/149‑3 (MQ=255)

         aGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCattat                                                                                                                                           >  1:1593304/1‑149 (MQ=255)

          gCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCattatt                                                                                                                                          <  1:315062/149‑1 (MQ=255)

          gCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCattatt                                                                                                                                          <  2:1094128/149‑1 (MQ=255)

            aTCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGa                                                                                                                                        <  1:168215/149‑1 (MQ=255)

                 ttATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACtttt                                                                                                                                   >  2:1385292/1‑149 (MQ=255)

                 ttATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGACTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACtttt                                                                                                                                   <  2:1057709/149‑1 (MQ=255)

                   aTCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCt                                                                                                                                 >  2:462376/1‑149 (MQ=255)

                   aTCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCt                                                                                                                                 >  1:568970/1‑149 (MQ=255)

                   aTCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCt                                                                                                                                 <  1:1561637/149‑1 (MQ=255)

                    tCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATTCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCtt                                                                                                                                >  1:220854/1‑148 (MQ=255)

                      aaaTTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGat                                                                                                                              >  1:96732/1‑149 (MQ=255)

                       aaTTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGata                                                                                                                             <  1:1260955/149‑1 (MQ=255)

                        aTTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGatat                                                                                                                            >  1:595999/1‑149 (MQ=255)

                         ttCATAAAGCTAAGTAATCTAGTCCCTAACTGGTAGTTACTTAACCTCATGAATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATTCACTTTCTTCGATTTTATGCTATAAATAGCATTATTGACTTTTCTGATAtt                                                                                                                           <  1:24184/149‑1 (MQ=255)

                          tCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATAttt                                                                                                                          >  1:337860/1‑149 (MQ=255)

                           cATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTa                                                                                                                         <  2:212087/149‑1 (MQ=255)

                             tataGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGa                                                                                                                       >  2:724239/1‑149 (MQ=255)

                             tataGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGa                                                                                                                       <  2:658540/149‑1 (MQ=255)

                             tataGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAAATAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGa                                                                                                                       >  1:582479/1‑149 (MQ=255)

                              ataGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGAt                                                                                                                      <  1:1340708/149‑1 (MQ=255)

                                  cTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTaa                                                                                                                  <  1:1509092/149‑1 (MQ=255)

                                   taagtaagCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAAt                                                                                                                 >  1:512893/1‑149 (MQ=255)

                                   taagtaagCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAAt                                                                                                                 <  1:788854/149‑1 (MQ=255)

                                     agtaagCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATaa                                                                                                               >  1:188169/1‑148 (MQ=255)

                                       taagCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATacac                                                                                                             <  2:601820/149‑1 (MQ=255)

                                              gTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGt                                                                                                      <  2:911554/149‑1 (MQ=255)

                                               tCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTc                                                                                                     <  2:764228/149‑1 (MQ=255)

                                                 ccGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACTAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCaa                                                                                                   >  1:1485220/1‑149 (MQ=255)

                                                 ccGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAATTAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCaa                                                                                                   >  2:287549/1‑149 (MQ=255)

                                                  cGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAAt                                                                                                  >  1:710171/1‑149 (MQ=255)

                                                        ggTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACAt                                                                                            <  2:986319/149‑1 (MQ=255)

                                                         gTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATg                                                                                           >  1:488953/1‑149 (MQ=255)

                                                               acctAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAggg                                                                                     <  2:685637/146‑1 (MQ=255)

                                                               aCTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAggg                                                                                     <  2:268267/149‑1 (MQ=255)

                                                                cTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTTCATTATGCACTTTCTTCGATGTTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGt                                                                                    <  1:891221/149‑1 (MQ=255)

                                                                cTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAg                                                                                                                                                 <  1:1512280/88‑1 (MQ=255)

                                                                cTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAg                                                                                                                                                 >  2:1512280/1‑88 (MQ=255)

                                                                   aaCCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTaaa                                                                                 <  2:377289/149‑1 (MQ=255)

                                                                   aaCCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTaaa                                                                                 >  1:1361103/1‑149 (MQ=255)

                                                                    aCCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAAt                                                                                <  2:826376/149‑1 (MQ=255)

                                                                     ccTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATg                                                                               >  2:1320483/1‑149 (MQ=255)

                                                                     ccTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAACAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATg                                                                               >  1:648416/1‑149 (MQ=255)

                                                                       tCATGTCTATATATATCATCTAACAAGAAATAAATAACAAGCCTACATTCTGCCCTTTCTTCGATGTGCTGCTATAAATAGCATTATTGCCTTTTCTGAGATTTAGATCTAATACACTTTACGTCAATAAACATGACAGTGTAAATGTc                                                                             <  2:168593/149‑1 (MQ=255)

                                                                           gCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCaaaa                                                                         <  2:1176191/149‑1 (MQ=255)

                                                                           gCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGATATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATAAACTTTCCGTCAATAAACATGACAGTGTAAATGTCaaaa                                                                         >  2:776480/1‑149 (MQ=255)

                                                                             atatatatatCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACa                                                                       >  1:117019/1‑149 (MQ=255)

                                                                                tatatatCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTTTAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTg                                                                    <  2:1245592/149‑1 (MQ=255)

                                                                                tatatatCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTg                                                                    <  2:1108093/149‑1 (MQ=255)

                                                                                    tatCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCaaa                                                                >  2:118970/1‑149 (MQ=255)

                                                                                     atcatcGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAg                                                               >  2:488433/1‑149 (MQ=255)

                                                                                     atcatcGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAg                                                               >  2:909785/1‑149 (MQ=255)

                                                                                       catcGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTa                                                             <  2:1287279/149‑1 (MQ=255)

                                                                                         tcGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATATCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGa                                                           <  2:199270/149‑1 (MQ=255)

                                                                                          cGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGaa                                                          >  2:408324/1‑149 (MQ=255)

                                                                                          cGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGaa                                                          >  1:151776/1‑149 (MQ=255)

                                                                                            aaaaaGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGaaaa                                                        >  2:1149530/1‑149 (MQ=255)

                                                                                                  aaaTAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAAc                                                  <  2:568970/149‑1 (MQ=255)

                                                                                                     taactaacAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTa                                               <  2:1005386/149‑1 (MQ=255)

                                                                                                     taactaacAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTa                                               <  1:1301734/149‑1 (MQ=255)

                                                                                                     taactaacAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATATAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATAAACGTa                                               >  1:1223953/1‑149 (MQ=255)

                                                                                                     aaactaacAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTa                                               <  1:141412/148‑1 (MQ=255)

                                                                                                        ctaacAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCg                                            >  2:231477/1‑149 (MQ=255)

                                                                                                         taacAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGt                                           >  1:767325/1‑149 (MQ=255)

                                                                                                            aaaGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCAt                                        <  1:1116804/148‑1 (MQ=255)

                                                                                                              aGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATAGCAAAGTAGAAAAATCAACGTAGCGTCAttt                                      <  2:618312/149‑1 (MQ=255)

                                                                                                               gCCTACATTATTCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTATATATAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAAAGTAGCGTCATTTa                                     >  1:259880/1‑149 (MQ=255)

                                                                                                               gCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATTTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTa                                     >  2:78325/1‑149 (MQ=255)

                                                                                                                ccTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTaa                                    <  1:724239/149‑1 (MQ=255)

                                                                                                                 cTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAAc                                   <  1:1385292/149‑1 (MQ=255)

                                                                                                                 cTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAAc                                   <  2:512969/149‑1 (MQ=255)

                                                                                                                  tACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACt                                  >  2:1336496/1‑149 (MQ=255)

                                                                                                                  tACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGAAAAGTAGAAAAATCAACGTATCGTCATTTAACt                                  >  1:146313/1‑149 (MQ=255)

                                                                                                                        aTGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATa                                                                                                                >  1:385635/1‑65 (MQ=255)

                                                                                                                          gCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCa                          <  2:312889/149‑1 (MQ=255)

                                                                                                                          gCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCa                          <  2:160576/149‑1 (MQ=255)

                                                                                                                             cTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCatat                       <  2:316644/149‑1 (MQ=255)

                                                                                                                                 cttcGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCatatatat                   <  1:1149530/149‑1 (MQ=255)

                                                                                                                                  ttcGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATCCACTTTCCGTCAATAACCATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAAACAACGTAGCGTCATTAAACTGAAAAGCatatatata                  >  2:1160231/1‑148 (MQ=255)

                                                                                                                                  ttcGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATc                  >  2:1033872/1‑149 (MQ=255)

                                                                                                                                    gtcaataaTGCTATTTATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTg                >  1:850488/8‑149 (MQ=255)

                                                                                                                                     gATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCtca               >  1:411323/1‑147 (MQ=255)

                                                                                                                                       tGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATTTATCTGACg             >  2:555928/1‑149 (MQ=255)

                                                                                                                                       aataaTGCTATTTATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACg             >  2:615188/5‑149 (MQ=255)

                                                                                                                                           aTGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTAt         >  2:1333620/1‑149 (MQ=255)

                                                                                                                                            tGCTATTTATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAggg                                                                                     <  1:605098/72‑1 (MQ=255)

                                                                                                                                            tGCTATTTATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAggg                                                                                     >  2:605098/1‑72 (MQ=255)

                                                                                                                                              cTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTc      >  2:1482894/1‑149 (MQ=255)

                                                                                                                                               tataAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAGTCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACTTTATTTCa     >  2:881179/1‑149 (MQ=255)

                                                                                                                                               tataAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCa     <  2:710171/149‑1 (MQ=255)

                                                                                                                                               tataAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCa     <  2:1485220/149‑1 (MQ=255)

                                                                                                                                               tataAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCa     <  1:488433/149‑1 (MQ=255)

                                                                                                                                               tataAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCa     >  2:49628/1‑149 (MQ=255)

                                                                                                                                               tataAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCa     >  2:277185/1‑149 (MQ=255)

                                                                                                                                               tataAATAGCATTATTGACTTTTATGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGAAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCa     >  2:299055/1‑149 (MQ=255)

TTCAACATCAGCATCTCTTATCAAATTCATATAGCTAAGTAAGCTAGTCCCGAACTGGTAGTTACTTAACCTCATGCATATATATATCATCGAAAAAGAAATAACTAACAAGCCTACATTATGCACTTTCTTCGATGGTATGCTATAAATAGCATTATTGACTTTTCTGATATTTAGATCTAATACACTTTCCGTCAATAAACATGACAGGGTAAATGTCAAAACATTGCAAAGTAGAAAAATCAACGTAGCGTCATTTAACTGAAAAGCATATATATCTGACGTTATTTCATAT  >  NC_000913/2106813‑2107107

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