BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,654	T→A	9.6%	intergenic (‑64/+40)	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,654	0	T	A	9.6%	331.5 / 20.4	115	intergenic (‑64/+40)	eptA/adiC	phosphoethanolamine transferase EptA/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (48/56); new base A (6/5); total (54/61)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 117 total aligned reads displayed.

TTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAA  >  NC_000913/4335505‑4335803

ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGt                                                                                                                                                       >  1:910242/1‑150 (MQ=255)

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

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

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

        gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAgg                                                                                                                                               <  2:376032/150‑1 (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)

             agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCGAA‑AT‑TAAGCCAACCTTAAGAACTTAAGGTTGGc                                                                                                                                          <  1:1526044/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)

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

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

                   aaaaaGCGGCCAACAATAGCCAGGCGAGTAAAATCAAAGAGCGTCTTTTTAGTAGGCGCTTCAACATCACGGGGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCCAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTAt                                                                                                                                    <  2:1089775/150‑1 (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)

                            ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTgg                                                                                                                           >  1:1264928/1‑150 (MQ=255)

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

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

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

                                     gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgggg                                                                                                                  >  1:855762/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)

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

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

                                                aaaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCttt                                                                                                       >  1:810713/1‑150 (MQ=255)

                                                 aaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTg                                                                                                      <  2:1351242/150‑1 (MQ=255)

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

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

                                                    tCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCtt                                                                                                   <  2:1453810/150‑1 (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)

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

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

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

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

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

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

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

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

                                                                        agGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATa                                                                               >  2:352237/1‑150 (MQ=255)

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

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

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

                                                                                 aaCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttttt                                                                      <  1:279052/150‑1 (MQ=255)

                                                                                   cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGTAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                                    <  2:593312/150‑1 (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)

                                                                                              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)

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

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

                                                                                                ccATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTa                                                       <  1:925150/150‑1 (MQ=255)

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

                                                                                                 cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                      <  2:490718/150‑1 (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)

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

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

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

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

                                                                                                                  atatGCTCGCAAAGCAAT‑AT‑TAAGCCAACCTTAAGTTCTTTAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGc                                     <  1:537144/150‑1 (MQ=255)

                                                                                                                    atGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                                   >  1:1084857/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)

                                                                                                                                             cAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAAt            >  2:1491813/1‑150 (MQ=255)

                                                                                                                                              aaCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATg           <  1:685024/150‑1 (MQ=255)

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

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

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

                                                                                                                                                 cTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAcc        <  2:833981/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)

TTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑AT‑TAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCACATAA  >  NC_000913/4335505‑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: `‑`