BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,655	T→A	8.6%	intergenic (‑65/+39)	eptA ← / ← adiC	phosphoethanolamine transferase EptA/arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,335,655	0	T	A	8.6%	339.8 / 15.2	115	intergenic (‑65/+39)	eptA/adiC	phosphoethanolamine transferase EptA/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (48/57); new base A (6/4); total (54/61)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

CAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAA  >  NC_000913/4335507‑4335803

cAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTc                                                                                                                                                     >  2:537201/1‑150 (MQ=255)

cAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTc                                                                                                                                                     >  1:1481344/1‑150 (MQ=255)

      gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAgg                                                                                                                                               <  2:376032/150‑1 (MQ=255)

      gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAgg                                                                                                                                               >  2:1513744/1‑150 (MQ=255)

       aTAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGATCTTAAGGt                                                                                                                                              <  1:387688/150‑1 (MQ=255)

          gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAACATAAAAGCCAACCTTAAGAACTTAAGGtt                                                                                                                                             <  2:285682/150‑1 (MQ=39)

          gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGAACTTAAGGTTgg                                                                                                                                           <  2:109224/150‑1 (MQ=255)

           agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                                                                          <  1:659137/150‑1 (MQ=255)

           agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                                                                          <  2:875874/150‑1 (MQ=255)

             atatAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCtt                                                                                                                                        >  2:918492/1‑150 (MQ=255)

              tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCttt                                                                                                                                       <  1:697023/150‑1 (MQ=255)

              tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGAACTTAAGGTTGGCtta                                                                                                                                       <  1:459430/150‑2 (MQ=255)

               ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCtttt                                                                                                                                      <  2:1078842/150‑1 (MQ=255)

                taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTa                                                                                                                                     <  1:948311/150‑1 (MQ=255)

                taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTa                                                                                                                                     >  1:787535/1‑150 (MQ=255)

                taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGATAGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTa                                                                                                                                     >  1:1340657/1‑150 (MQ=255)

                 aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTAt                                                                                                                                    >  1:386949/1‑150 (MQ=255)

                  aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATg                                                                                                                                   >  1:1447020/1‑150 (MQ=255)

                    aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGATCAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGtt                                                                                                                                 <  1:925641/150‑1 (MQ=255)

                        ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCt                                                                                                                             <  2:910242/150‑1 (MQ=255)

                        ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGGGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCt                                                                                                                             <  1:972391/150‑1 (MQ=255)

                         gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTg                                                                                                                            >  1:593312/1‑150 (MQ=255)

                                aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                     >  1:833981/1‑150 (MQ=255)

                                aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                     >  1:748892/1‑150 (MQ=255)

                                   gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgggg                                                                                                                  >  1:855762/1‑150 (MQ=255)

                                   gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgggg                                                                                                                  >  2:387324/1‑150 (MQ=255)

                                       ggCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAAtt                                                                                                              <  1:161403/150‑1 (MQ=255)

                                       ggCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAAtt                                                                                                              >  1:573896/1‑150 (MQ=255)

                                        gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTa                                                                                                             >  1:1490453/1‑150 (MQ=255)

                                           aGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATc                                                                                                          >  2:48247/1‑150 (MQ=255)

                                            gTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCt                                                                                                         >  1:490718/1‑150 (MQ=255)

                                             tAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCtt                                                                                                        >  2:807496/1‑150 (MQ=255)

                                              aaaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCttt                                                                                                       <  2:1537638/150‑1 (MQ=255)

                                               aaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTg                                                                                                      >  1:1114722/1‑150 (MQ=255)

                                                  tCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCtt                                                                                                   <  2:1453810/150‑1 (MQ=255)

                                                  tCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCtt                                                                                                   >  1:759665/1‑150 (MQ=255)

                                                   cAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTa                                                                                                  <  2:952930/150‑1 (MQ=255)

                                                     aaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                                >  2:364826/1‑150 (MQ=255)

                                                      agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                               >  2:279052/1‑150 (MQ=255)

                                                       gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTgg                                                                                              <  2:1385306/150‑1 (MQ=255)

                                                        agGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                             <  2:1264928/150‑1 (MQ=255)

                                                         gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                            <  1:90478/150‑1 (MQ=255)

                                                         gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                            <  2:810713/150‑1 (MQ=255)

                                                         gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                            <  1:387324/150‑1 (MQ=255)

                                                         gggTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                            <  2:759665/150‑1 (MQ=255)

                                                            tCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                     <  1:508584/138‑1 (MQ=255)

                                                                 ttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                                                                    <  1:537201/150‑1 (MQ=255)

                                                                  tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                                                                   <  2:151894/150‑1 (MQ=255)

                                                                  tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                                                                   <  2:496465/150‑1 (MQ=255)

                                                                  tagtagGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAgg                                                                                   <  1:441444/150‑1 (MQ=255)

                                                                       ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATAc                                                                              >  2:1098097/1‑150 (MQ=255)

                                                                          gcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACggg                                                                           <  2:1447020/150‑1 (MQ=255)

                                                                          gcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACggg                                                                           >  1:390928/1‑150 (MQ=255)

                                                                            ttCAACGTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtt                                                                         >  1:216377/1‑150 (MQ=255)

                                                                              cAACATCAAGTTGTTTCAATCTAACAAACTCCGCATATGCTCGAAAAGAAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtttt                                                                       >  2:600844/1‑150 (MQ=255)

                                                                                 cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGTAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                    <  2:593312/150‑1 (MQ=255)

                                                                                 cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                    >  1:1175187/1‑150 (MQ=255)

                                                                                  aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATg                                                                   <  1:44514/150‑1 (MQ=255)

                                                                                   tCACGGTGTTTCCATCGAACACACTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGt                                                                  <  2:1376716/150‑1 (MQ=255)

                                                                                     aCGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAg                                                                >  1:1024649/1‑150 (MQ=255)

                                                                                          gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGt                                                           >  1:89561/1‑150 (MQ=255)

                                                                                          gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGt                                                           <  1:26228/150‑1 (MQ=255)

                                                                                            ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TTAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg                                                         <  1:450083/150‑1 (MQ=255)

                                                                                            ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTGAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg                                                         <  1:1435388/150‑1 (MQ=255)

                                                                                            ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg                                                         <  2:1445521/150‑1 (MQ=255)

                                                                                            ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg                                                         <  1:1513744/150‑1 (MQ=255)

                                                                                             tCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGt                                                        <  1:997374/150‑1 (MQ=255)

                                                                                             tCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGt                                                        >  1:369385/1‑150 (MQ=255)

                                                                                               cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                      <  2:490718/150‑1 (MQ=255)

                                                                                               cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                      >  1:261235/1‑150 (MQ=255)

                                                                                                aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAt                                                     <  2:993526/150‑1 (MQ=255)

                                                                                                aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAt                                                     >  1:291966/1‑150 (MQ=255)

                                                                                                 tCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAtt                                                    <  2:1481344/150‑1 (MQ=255)

                                                                                                  cGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTc                                                   <  1:793053/150‑1 (MQ=255)

                                                                                                   gAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCa                                                  >  2:793694/1‑150 (MQ=255)

                                                                                                        aaGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCa                                             >  2:1162069/1‑150 (MQ=255)

                                                                                                          gTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATa                                           >  1:833743/1‑150 (MQ=255)

                                                                                                          gTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATa                                           >  1:659829/1‑150 (MQ=255)

                                                                                                                  atGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                                   >  1:1084857/1‑150 (MQ=255)

                                                                                                                    gCTCGCAAAGCAAA‑AT‑TACGCCAACCTTAAGTTCTTACGGTTGGATTTTATGTTTGCTGGATTCAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGGTGTTATGTAGCCGGTGGTAATTCAGCGCATACATGGCGGtg                                 <  2:527442/150‑1 (MQ=255)

                                                                                                                    gCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtg                                 >  2:367270/1‑150 (MQ=255)

                                                                                                                     cTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                                <  1:48247/150‑1 (MQ=255)

                                                                                                                          aaagcaaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGAcc                           <  1:929320/150‑1 (MQ=255)

                                                                                                                              caaa‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCa                       <  2:89561/150‑1 (MQ=255)

                                                                                                                               aaa‑aT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTTAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAg                      >  1:1059101/1‑150 (MQ=255)

                                                                                                                               aaa‑aT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAg                      <  2:573896/150‑1 (MQ=255)

                                                                                                                               aaa‑aT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAg                      >  2:543351/1‑150 (MQ=255)

                                                                                                                                 a‑aT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAggg                    <  2:1175187/150‑1 (MQ=255)

                                                                                                                                  cataaaaGCCAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAgg                     >  2:766074/6‑150 (MQ=255)

                                                                                                                                         gCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaa              <  2:291692/150‑1 (MQ=255)

                                                                                                                                          ccAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa             >  1:74306/1‑150 (MQ=255)

                                                                                                                                          ccAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa             >  1:174034/1‑150 (MQ=255)

                                                                                                                                          ccAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa             >  1:1163155/1‑150 (MQ=255)

                                                                                                                                            aaCCTTAAGAACTTAAGGTTGGCTTTTATGTTTGATGTATTAAGGGGAATTAATCTTTGCTTATTAGTTCATCTAAGGGATACGCATTTTTATTTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATg           >  2:1259737/1‑150 (MQ=255)

                                                                                                                                               cTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAcc        <  2:787535/150‑1 (MQ=255)

                                                                                                                                                ttAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACca       <  1:16573/150‑1 (MQ=255)

                                                                                                                                                  aaGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACcaca     <  1:600844/150‑1 (MQ=255)

                                                                                                                                                    gTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATa   >  1:907204/1‑150 (MQ=255)

                                                                                                                                                    gTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATa   >  1:1106753/1‑150 (MQ=255)

                                                                                                                                                     ttCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATaa  <  2:1049794/150‑1 (MQ=255)

CAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAA  >  NC_000913/4335507‑4335803

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