BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,284,655	T→A	30.9%	Y83F (TAT→TTT)	actP ←	acetate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,284,655	0	T	A	30.9%	139.1 / 83.6	100	Y83F (TAT→TTT)	actP	acetate transporter
Reads supporting (aligned to +/- strand): ref base T (31/38); new base A (15/16); total (46/54)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 107 total aligned reads displayed.

GTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCT  >  NC_000913/4284508‑4284802

gtgtAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACatat                                                                                                                                                     <  6:259916/148‑1 (MQ=255)

 tgtAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACAtata                                                                                                                                                    <  3:259137/148‑1 (MQ=255)

      cgcCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTccc                                                                                                                                               >  6:18098/1‑148 (MQ=255)

       gcCCCAGGTTACGCAGACGTTCGGCAATCCGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTcccc                                                                                                                                              <  6:275745/148‑1 (MQ=255)

          ccAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGGCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCgg                                                                                                                                            >  5:274421/1‑147 (MQ=255)

           cAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGc                                                                                                                                           >  8:334019/1‑147 (MQ=255)

            aGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGTACATATAGTCCCCGGCaa                                                                                                                                         >  5:18507/1‑148 (MQ=255)

               ttACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTccc                                                                                                                                               >  5:177081/1‑139 (MQ=255)

               ttACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTccc                                                                                                                                               <  6:177081/139‑1 (MQ=255)

               ttACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCg                                                                                                                                      >  6:172591/1‑148 (MQ=255)

               ttACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCg                                                                                                                                      >  7:101194/1‑148 (MQ=255)

                tACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGc                                                                                                                                     >  5:285503/1‑148 (MQ=255)

                 aCGCAGCAGTTCGGCAATCAGGAACAAAATGATCGGCAAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGAAATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCTAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGc                                                                                                                                     <  2:88592/147‑1 (MQ=255)

                  cGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGcc                                                                                                                                    <  8:6043/147‑1 (MQ=255)

                    cAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGTGAGTAAATTAAGCCGTCATAGCCGGAGGTAACCACCAGCGCGGAGATCCGCAAGAATGAGGCGGCGGACATATAGTCCCCGGAAATCGCCAGc                                                                                                                                 >  3:251158/1‑148 (MQ=255)

                        cGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGt                                                                                                                             >  2:232832/1‑148 (MQ=255)

                          ttCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCCATCGCCAGCCCGTTc                                                                                                                           >  5:2675/1‑148 (MQ=255)

                              gCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGa                                                                                                                       <  8:259179/148‑1 (MQ=255)

                                  tCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGtt                                                                                                                            >  2:199078/1‑139 (MQ=255)

                                  tCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGtt                                                                                                                            <  1:199078/139‑1 (MQ=255)

                                  tCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGcc                                                                                                                   <  1:172774/148‑1 (MQ=255)

                                  tCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGcc                                                                                                                   <  3:232346/148‑1 (MQ=255)

                                   cAGGAACAAAATGATCGGCCAGCCCACCAGTAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTGAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGTCAATCGCCAGCCCGTTCTGGAAGc                                                                                                                    >  1:255917/1‑146 (MQ=255)

                                    aGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAg                                                                                                                 <  6:285503/148‑1 (MQ=255)

                                        aCAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGat                                                                                                             <  1:123540/148‑1 (MQ=255)

                                        aCAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGat                                                                                                             <  1:156693/148‑1 (MQ=255)

                                        aCAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGat                                                                                                             <  8:262629/148‑1 (MQ=255)

                                         ctaaaTGATCGGCCAGCCCACCAGGTAGCCCAGCGAGTAAATTAAGCCGTGATAGCCGGTGGTAAACACCAGCGCGGAGATCCCCAAGTATGAGGCGGGGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGata                                                                                                            <  3:292181/146‑1 (MQ=255)

                                            aaTGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTg                                                                                                         <  7:334019/148‑1 (MQ=255)

                                            aaTGATCGGCAAGCCCACCAGGAAGCCCAGCGAGAAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTg                                                                                                         >  2:210836/1‑148 (MQ=255)

                                             aTGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTgc                                                                                                        >  1:228084/1‑148 (MQ=255)

                                               gATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTgccg                                                                                                      >  2:130333/1‑148 (MQ=255)

                                                aTCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGcc                                                                                                                   >  6:133967/1‑134 (MQ=255)

                                                aTCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGcc                                                                                                                   <  5:133967/134‑1 (MQ=255)

                                                    gCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGc                                                                                                 <  8:143428/148‑1 (MQ=255)

                                                      cAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATATCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGAGGACATAAAGTCCACGGCAATCGCCAGCCCGTTGTGGAAGCCAGTGATATTTCCTCCTGCgt                                                                                               >  5:87930/1‑147 (MQ=255)

                                                       aGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGgt                                                                                              <  1:232832/148‑1 (MQ=255)

                                                        gCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGcc                                                                                                                   <  3:29133/126‑1 (MQ=255)

                                                        gCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGcc                                                                                                                   >  4:29133/1‑126 (MQ=255)

                                                         cccACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGgtgt                                                                                            <  5:172591/148‑1 (MQ=255)

                                                           caccaGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTgtag                                                                                          >  5:168784/1‑148 (MQ=255)

                                                            accaGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTgtagt                                                                                         <  1:130333/148‑1 (MQ=255)

                                                            accaGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTgtagt                                                                                         <  8:101194/148‑1 (MQ=255)

                                                             ccaGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTgtag                                                                                          >  4:127693/1‑146 (MQ=255)

                                                              caGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGtagtag                                                                                       <  6:274421/148‑1 (MQ=255)

                                                                ggAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAg                                                                                                                     >  2:44893/1‑116 (MQ=255)

                                                                ggAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAg                                                                                                                     <  1:44893/116‑1 (MQ=255)

                                                                 gAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCg                                                                                    <  3:250694/148‑1 (MQ=255)

                                                                    gCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTa                                                                                 <  7:244631/148‑1 (MQ=255)

                                                                    gCCCAGCGACTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGCCGGCGGACATTTAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATTTTGCCGCCTGCGGTGTAGTAGTCGCTa                                                                                 <  3:147572/148‑1 (MQ=255)

                                                                     cccAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGAGCGGAGATAACCAAGAATGAGTCGGCGGACATATAGTCCCCGGCAATAGCCATCACGTTATGGAAGCCAGTGATATTGAAGACTGAGGTGTATTATTCGCt                                                                                  >  8:75253/1‑146 (MQ=255)

                                                                      ccAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACg                                                                               <  5:78507/148‑1 (MQ=255)

                                                                      ccAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTAc                                                                                >  1:294726/1‑147 (MQ=255)

                                                                      ccAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACg                                                                               <  2:228084/148‑1 (MQ=255)

                                                                        aGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTAc                                                                                >  6:281786/1‑145 (MQ=255)

                                                                        aGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACgag                                                                             <  2:294726/148‑1 (MQ=255)

                                                                          cGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTgtagt                                                                                         <  4:69884/134‑1 (MQ=255)

                                                                          cGAGTAAAGTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTgtagt                                                                                         >  3:69884/1‑134 (MQ=255)

                                                                             gTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGt                                                                        <  5:184493/148‑1 (MQ=255)

                                                                               aaaTTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTgc                                                                                                        >  6:86670/1‑114 (MQ=255)

                                                                               aaaTTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTgc                                                                                                        <  5:86670/114‑1 (MQ=255)

                                                                                  ttAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCAAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTAcgcg                                                                   >  1:97721/1‑148 (MQ=255)

                                                                                    atcccGTCATACCCGTAGGTAACCCCCAGCGCCGCGA‑GCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGGGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGtt                                                                 <  4:128870/144‑1 (MQ=255)

                                                                                    aatccGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGt                                                                  <  5:18098/144‑1 (MQ=255)

                                                                                       ccGTCATAGCCGGAGGTGAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGTTGTAGTAGTCGCTACGAGAACGTACGCGTTTTg                                                              <  2:255917/148‑1 (MQ=255)

                                                                                        cGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGa                                                             <  6:168784/148‑1 (MQ=255)

                                                                                              aGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGcc                                                         >  4:70515/1‑146 (MQ=255)

                                                                                                ccGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAg                                                      >  6:73475/1‑147 (MQ=255)

                                                                                                ccGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGGTTTGATGCCCAGt                                                     <  2:62170/148‑1 (MQ=255)

                                                                                                   taGGTAAACACCAGCTCGGAGATCCCCAAGAATGAGGCGGCGGACATATATTCCCCGGCAATCGCCAGCCCGTTCTGGAAGAAATTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAgg                                                  <  7:75253/147‑1 (MQ=255)

                                                                                                       tAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGt                                                                        >  6:261538/1‑122 (MQ=255)

                                                                                                       tAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGt                                                                        <  5:261538/122‑1 (MQ=255)

                                                                                                       tAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGt                                                     >  6:245192/1‑141 (MQ=255)

                                                                                                       tAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGt                                                     <  5:245192/141‑1 (MQ=255)

                                                                                                          acacCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGt                                                 >  1:78090/1‑142 (MQ=255)

                                                                                                          acacCAGCCCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGt                                                 <  2:78090/142‑1 (MQ=255)

                                                                                                           cacCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCg                                          <  6:103369/148‑1 (MQ=255)

                                                                                                            acCAGCGCGGAGATCCCAAAGAATGAGGCGGCGGACATATATTCAACGGCAATCGCAAGCAAGTTCTGGAAGCCAGTGATATTGCCGACTGCGGTGTATTAGTCGCTACTAGAACGTACGCGTTTTTATGCCCAGTAGGTAATGACGa                                         >  8:309033/1‑148 (MQ=255)

                                                                                                             ccAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAg                                        <  3:224444/148‑1 (MQ=255)

                                                                                                              cAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTa                                                    >  5:50856/1‑135 (MQ=255)

                                                                                                              cAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTa                                                    <  6:50856/135‑1 (MQ=255)

                                                                                                               aGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCg                                      >  4:268720/1‑148 (MQ=255)

                                                                                                               aGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCg                                      >  6:6973/1‑148 (MQ=255)

                                                                                                                  gcgGAGATCGCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTaa                                   <  4:251158/148‑1 (MQ=255)

                                                                                                                    ggAGATCCCCAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTaaa                                  >  8:261522/1‑147 (MQ=255)

                                                                                                                       gaTCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACacg                              <  3:83883/148‑1 (MQ=255)

                                                                                                                         tCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCgg                                                                                               <  2:315531/81‑1 (MQ=255)

                                                                                                                         tCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCgg                                                                                               <  2:319994/81‑1 (MQ=255)

                                                                                                                         tCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCgg                                                                                               >  1:319994/1‑81 (MQ=255)

                                                                                                                         tCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCgg                                                                                               >  1:315531/1‑81 (MQ=255)

                                                                                                                            ccAAGAATGAGGCGGCGGACATAAAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACAcgccga                          >  5:212577/1‑147 (MQ=255)

                                                                                                                                   tGAGGCGGCGGACATATAGTCCCCTGCAATCTCCAGCCAGTTCTGGAAGCCAGTTATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGAGTTTTGATGCCCAGTAGGTAATGCCGAGCTTAAACACTACGAAAATCAgg                  <  7:309033/148‑1 (MQ=255)

                                                                                                                                      ggcggcGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTa                                                <  7:303756/115‑1 (MQ=255)

                                                                                                                                      ggcggcGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTa                                                >  8:303756/1‑115 (MQ=255)

                                                                                                                                              aCATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATaataat       <  3:127693/148‑1 (MQ=255)

                                                                                                                                                atatAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCg     <  3:70515/148‑1 (MQ=255)

                                                                                                                                                 tataGTCACCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGAAATTGCCGCCTGCGTTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCAGAGCGTAAACACGACGAAAAACAGGAACATAAAAATCGc    >  6:8826/1‑148 (MQ=255)

                                                                                                                                                   taGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGaaa                        <  3:290716/126‑1 (MQ=255)

                                                                                                                                                   taGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGaaa                        >  4:290716/1‑126 (MQ=255)

                                                                                                                                                   taGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCt  >  3:211487/1‑148 (MQ=255)

GTGTAGCGCCCCAGGTTACGCAGACGTTCGGCAATCAGGAACAAAATGATCGGCCAGCCCACCAGGAAGCCCAGCGAGTAAATTAAGCCGTCATAGCCGGAGGTAAACACCAGCGCGGAGATCCCCAAGAATGAGGCGGCGGACATATAGTCCCCGGCAATCGCCAGCCCGTTCTGGAAGCCAGTGATATTGCCGCCTGCGGTGTAGTAGTCGCTACGAGAACGTACGCGTTTTGATGCCCAGTAGGTAATGCCGAGCGTAAACACGACGAAAATCAGGAACATAATAATCGCCT  >  NC_000913/4284508‑4284802

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