BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,654	T→A	5.7%	intergenic (‑64/+40)	eptA ← / ← adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,335,654	0	T	A	5.7%	799.2 / 24.4	245	intergenic (‑64/+40)	eptA/adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (107/124); new base A (3/11); total (110/135)
Fisher's exact test for biased strand distribution p-value = 9.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 251 total aligned reads displayed.

AGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCA  >  NC_000913/4335508‑4335798

 ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCt                                                                                                                                                  >  7:337683/1‑149 (MQ=255)

  gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt                                                                                                                                                 <  1:37199/149‑1 (MQ=255)

     gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg                                                                                                                                              <  6:112591/149‑1 (MQ=255)

      aTAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAgg                                                                                                                                             >  8:4913/1‑149 (MQ=255)

         gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                           <  8:147230/148‑1 (MQ=255)

         gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTg                                                                                                                                          <  7:292120/149‑1 (MQ=255)

         gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTg                                                                                                                                          <  2:312429/149‑1 (MQ=255)

          agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTgg                                                                                                                                         >  6:311997/1‑149 (MQ=255)

           gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGc                                                                                                                                        <  6:180677/149‑1 (MQ=255)

           gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTATTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAGGTTGGc                                                                                                                                        <  1:103564/149‑1 (MQ=255)

            atatAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                                                                       >  8:274297/1‑149 (MQ=255)

             tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTg                                                                                                                                          <  5:342390/145‑1 (MQ=255)

             tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCtt                                                                                                                                      <  2:145360/149‑1 (MQ=255)

               taAAACGCCGCCAACAATCGCCAGGCGAGTAACTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCtt‑‑tt                                                                                                                                  <  4:230374/149‑1 (MQ=255)

                  aaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATg                                                                                                                               >  7:318496/1‑149 (MQ=255)

                       ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGc                                                                                                                         <  1:66822/149‑1 (MQ=255)

                          caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTgg                                                                                                                      >  7:97555/1‑149 (MQ=255)

                           aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGATGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGa                                                                                                                     <  1:174630/149‑1 (MQ=255)

                             caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGAtt                                                                                                                   <  5:453211/149‑1 (MQ=255)

                                   ccAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAggg                                                                                                              >  4:377593/1‑148 (MQ=255)

                                     aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTgg                                                                                                                                         <  7:347152/122‑1 (MQ=255)

                                     aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAg                                                                                                                >  1:131472/1‑144 (MQ=255)

                                     aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGaa                                                                                                           >  1:297264/1‑149 (MQ=255)

                                     aGGCGAGTAAAATCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGaa                                                                                                           <  3:371310/149‑1 (MQ=255)

                                      ggCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAAt                                                                                                          <  5:15697/149‑1 (MQ=255)

                                          aGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTaa                                                                                                                 <  1:458085/138‑1 (MQ=255)

                                            tAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCt                                                                                                    >  1:150788/1‑149 (MQ=255)

                                                  cAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCt                                                                                               >  4:285487/1‑148 (MQ=255)

                                                    aaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                            >  6:83947/1‑149 (MQ=255)

                                                     agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                           >  5:242777/1‑149 (MQ=255)

                                                     agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                           >  3:250006/1‑149 (MQ=255)

                                                     agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                           >  5:321236/1‑149 (MQ=255)

                                                      gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtta                                                                                          <  2:459570/149‑2 (MQ=255)

                                                      gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                          <  4:195657/149‑1 (MQ=255)

                                                        gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                        <  5:9160/149‑1 (MQ=255)

                                                        gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                        <  7:180374/149‑1 (MQ=255)

                                                         ggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                       <  8:337683/149‑1 (MQ=255)

                                                             tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGa                                                                                                            <  8:91265/124‑1 (MQ=255)

                                                             tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGa                                                                                                            >  7:91265/1‑124 (MQ=255)

                                                             tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATc                                                                                   >  8:320317/1‑149 (MQ=255)

                                                                 tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTaa                                                                                <  6:361836/148‑1 (MQ=255)

                                                                  agtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATc                                                                                   <  1:91177/144‑1 (MQ=255)

                                                                  agtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATA‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                                                               <  6:148240/148‑1 (MQ=255)

                                                                   gtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                                                              <  5:311997/148‑1 (MQ=255)

                                                                      ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATa                                                                          >  3:137046/1‑149 (MQ=255)

                                                                        cgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACg                                                                        <  3:285487/149‑1 (MQ=255)

                                                                             cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGt                                                                     <  4:137046/147‑1 (MQ=255)

                                                                               aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                >  4:160741/1‑118 (MQ=255)

                                                                               aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                <  3:160741/118‑1 (MQ=255)

                                                                               aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                          >  3:89400/1‑124 (MQ=255)

                                                                               aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                          <  4:89400/124‑1 (MQ=255)

                                                                               aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                                                                 <  4:271908/149‑1 (MQ=255)

                                                                                 aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCttt                                                                                                  <  8:118263/114‑1 (MQ=255)

                                                                                 aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCttt                                                                                                  >  7:118263/1‑114 (MQ=255)

                                                                                 aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                               <  2:86124/149‑1 (MQ=255)

                                                                                        tgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtttt                                                                  >  8:395068/1‑139 (MQ=255)

                                                                                         gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCg                                                        >  2:350587/1‑148 (MQ=255)

                                                                                         gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCg                                                        <  1:350587/148‑1 (MQ=255)

                                                                                              cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGt                                                   >  2:453236/1‑148 (MQ=255)

                                                                                              cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGt                                                   >  3:115072/1‑148 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg                                                                                                                                              >  8:23204/1‑59 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAAt                                                                                                          <  6:338284/92‑1 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCttt                                                                                                  <  6:32256/100‑1 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGc                                                          >  1:328574/1‑140 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGc                                                          <  2:328574/140‑1 (MQ=255)

                                                                                                tCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAt                                                >  5:419577/1‑149 (MQ=255)

                                                                                                tCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAt                                                >  4:123008/1‑149 (MQ=255)

                                                                                                    aCAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTg                                                                                                                       >  2:198947/1‑74 (MQ=255)

                                                                                                    aCAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGACAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCTTCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAtt                                               >  5:334116/1‑146 (MQ=255)

                                                                                                     cAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgc                                           <  6:242777/149‑1 (MQ=255)

                                                                                                        aGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCa                                        >  5:126188/1‑149 (MQ=255)

                                                                                                         gTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCa                                        >  3:328497/1‑148 (MQ=255)

                                                                                                                 atGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCg                               >  6:274322/1‑149 (MQ=255)

                                                                                                                    cTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGt                             <  6:438208/148‑1 (MQ=255)

                                                                                                                    cTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGt                             <  1:432400/148‑1 (MQ=255)

                                                                                                                      tccagcaaacataaAAGCCAACCTTAAGAACTTAAGGTTGGCTTAATTTTGCTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGt                             <  4:97347‑M1/135‑1 (MQ=255)

                                                                                                                       gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                                                                 >  7:241890/1‑109 (MQ=255)

                                                                                                                       gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                         >  8:352422/1‑149 (MQ=255)

                                                                                                                       gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                         >  7:228376/1‑149 (MQ=255)

                                                                                                                       gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                         <  7:7649/149‑1 (MQ=255)

                                                                                                                       gcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgat                          <  3:430961/148‑1 (MQ=255)

                                                                                                                        caaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                        <  4:31956/149‑1 (MQ=255)

                                                                                                                         aaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGAGTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTGTATGTAGCCGGTTGTAATTCAGCGCCTACATGGCGGtgatgag                       >  6:406659/1‑148 (MQ=255)

                                                                                                                           agcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCa                     <  1:21143/149‑1 (MQ=255)

                                                                                                                            gcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCAt                    >  5:246058/1‑149 (MQ=255)

                                                                                                                             aacataaAAGCCAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCAt                    <  5:258118‑M1/141‑1 (MQ=255)

                                                                                                                                aaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGTGGTGATGACCATCAgg                <  5:391445/149‑1 (MQ=255)

                                                                                                                                aaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCAGCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAgg                >  8:202085/1‑149 (MQ=255)

                                                                                                                                  ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt              <  6:246058/149‑1 (MQ=255)

                                                                                                                                      GCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                                                                 <  6:454656‑M1/94‑1 (MQ=255)

                                                                                                                                        cAACCTTAAGTTCTTAAGGTTGGATT‑‑TAATA‑ATAGCTCGATTAAGGAAAAATAATCTAAGAAATTTTATAAATCTAAGGTATACGTATAATTATAAAAAAACTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTTACaac        >  2:399757/1‑148 (MQ=255)

                                                                                                                                        CAACCTTAAGAACTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa        <  7:134984‑M1/149‑1 (MQ=255)

                                                                                                                                         AACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgat                          <  2:232222‑M1/130‑1 (MQ=255)

                                                                                                                                         AACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAAt       >  5:59920‑M1/1‑149 (MQ=255)

                                                                                                                                          ACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCg                               <  4:425121‑M1/124‑1 (MQ=255)

                                                                                                                                               AAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACca  <  7:414482‑M1/148‑1 (MQ=255)

                                                                                                                                                  ttCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                              >  2:409534/1‑117 (MQ=255)

                                                                                                                                                  ttCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                              <  1:409534/117‑1 (MQ=255)

                                                                                                                                                  ttCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                         >  7:374925/1‑122 (MQ=255)

                                                                                                                                                  ttCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatg                         <  8:374925/122‑1 (MQ=255)

AGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTT‑‑TTATG‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCA  >  NC_000913/4335508‑4335798

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