BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,654	T→A	5.8%	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	5.8%	389.2 / 14.7	139	intergenic (‑64/+40)	eptA/adiC	phosphoethanolamine transferase EptA/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (57/73); new base A (3/5); total (61/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-03

Only 100 of 143 total aligned reads displayed.

GATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAC  >  NC_000913/4335513‑4335796

gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAACGCAAA‑ATT‑AAGCCCAGCTTAAGtt                                                                                                                                               >  1:1109189/1‑143 (MQ=255)

gATAGAGATATAAAAAGCGACCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGtt                                                                                                                                               <  1:438370/143‑1 (MQ=255)

  tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGGAAA‑ATT‑AAGCCAACATTAAGTTCt                                                                                                                                             >  1:199253/1‑143 (MQ=255)

  tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCACA‑ATT‑AAGCCAACCTCACGTTCt                                                                                                                                             >  2:306557/1‑143 (MQ=255)

      gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAg                                                                                                                                         <  1:1230434/143‑1 (MQ=255)

         ataACACGCGTCAAACAATAGACAGGCGAGTAAATGGAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCTCGGTGTTTCCATTGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTACGTTCttttgttt                                                                                                                                      <  2:221143/143‑7 (MQ=255)

         ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGAACTTAAGGtt                                                                                                                                      <  1:1158949/143‑1 (MQ=255)

         ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGAACTTAAGGtt                                                                                                                                      <  1:1184452/143‑1 (MQ=255)

         ataAAAAGCGGCCAACAAGAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCGTAAGGtt                                                                                                                                      <  2:1142301/143‑1 (MQ=255)

          taAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGAACTTAAGGTTg                                                                                                                                     <  1:1273148/143‑1 (MQ=255)

           aaaaaGTGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTgg                                                                                                                                    <  1:534381/143‑1 (MQ=255)

           aaaaaGCGGCCAACCATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTgg                                                                                                                                    <  1:214961/143‑1 (MQ=255)

           aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTACGGTTgg                                                                                                                                    >  1:917174/1‑143 (MQ=255)

           aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTgg                                                                                                                                    <  1:724551/143‑1 (MQ=255)

            aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAACATAAAAGCCAACCTTAAGAACTTAAGGTTg                                                                                                                                     <  2:656048/143‑1 (MQ=38)

                 cGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ACT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTa                                                                                                                              >  1:731513/1‑143 (MQ=255)

                  ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTAt                                                                                                                             <  1:976562/143‑1 (MQ=255)

                    ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTGTTATGt                                                                                                                           >  1:397888/1‑143 (MQ=255)

                       acaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTg                                                                                                                        <  1:1193952/143‑1 (MQ=255)

                         aaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCt                                                                                                                      <  2:937377/143‑1 (MQ=255)

                         aaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTCCTTAAGGTTGGCTTTTATGTTTGCt                                                                                                                      >  2:1021808/1‑143 (MQ=255)

                         aaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCGCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCt                                                                                                                      <  1:1229665/143‑1 (MQ=255)

                            aGCCAGGCGAGTAAATTCAAAGAGGGTCTTTGTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTCTATGTTTGCTGGa                                                                                                                   <  1:552800/143‑1 (MQ=255)

                            aGCCAGGCGAGTAAATTCAAAGAGGGTCCTTTTAGTAGGCGCTCCAACATCACGGTTTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGa                                                                                                                   <  1:1069332/143‑1 (MQ=255)

                                aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑CTT‑AAGCCAACCTTAAGTTCTTAAGGCTGGCTTTTATGTTTGCTGGactga                                                                                                               >  1:587260/1‑139 (MQ=255)

                                      gTAAATTCAAAGAGGGTCTTTCTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                         <  1:645510/143‑1 (MQ=255)

                                       tAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTTCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTCTTATGTTTGCTGGATTAAGGGGAAt                                                                                                        >  1:152431/1‑143 (MQ=255)

                                       tAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAAt                                                                                                        <  1:1089500/143‑1 (MQ=255)

                                         aaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTa                                                                                                      <  1:1257817/143‑1 (MQ=255)

                                           ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAAt                                                                                                    <  1:712205/143‑1 (MQ=255)

                                              aaaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTCAGGGGTATTAATCtt                                                                                                 >  1:159313/1‑143 (MQ=255)

                                              aaaGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCtt                                                                                                 >  1:828597/1‑143 (MQ=255)

                                                     gTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                          >  2:166134/1‑143 (MQ=255)

                                                     gTCTTTTTAGTAGGCGCTTCAAAATAACGGTCTCTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAt                                                                                          >  1:879526/1‑143 (MQ=255)

                                                      tCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                         <  1:1127102/143‑1 (MQ=255)

                                                       cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTg                                                                                        <  1:1224372/143‑1 (MQ=255)

                                                        tttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTgg                                                                                       >  2:784514/1‑143 (MQ=255)

                                                         ttttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                      <  2:890500/143‑1 (MQ=255)

                                                          tttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                     >  1:993881/1‑143 (MQ=255)

                                                          tttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                     <  1:500609/143‑1 (MQ=255)

                                                          tttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGATCGCAAAGCAAA‑ATT‑GCGCCAACCTTAAGTTCTTAAGGTTGGCGTTTATGTTTGCTGGTGGACGGGGACTTAATCTTTGCTTATTGGtc                                                                                     >  1:614214/1‑142 (MQ=255)

                                                          tttAGTAGGCGCTTCAAACTCAAGGTGCTTCCACCGAACAAAGTGCACATATGCTCGAACTGAGAA‑TTC‑AAGCCATCATTAAGTTCTTAAGGTTGGCTTTTCAATTTGCTTGATTAAGGGGAATTAATCTTTGCTTATTGGTg                                                                                     <  2:325753/143‑1 (MQ=255)

                                                           ttAGTAGGCGCTTCAACGTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                                                                    <  2:115113/143‑1 (MQ=255)

                                                           ttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGc                                                                                    >  2:947816/1‑143 (MQ=255)

                                                               tagGCGCTTCAACATCACGGTGTTTCCATCGAACAAGGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCt                                                                                >  1:607393/1‑143 (MQ=255)

                                                                 ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTaa                                                                              <  1:644061/143‑1 (MQ=255)

                                                                  gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑GTT‑GAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                                                             <  2:1110131/143‑1 (MQ=255)

                                                                  gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAg                                                                             <  1:222829/143‑1 (MQ=255)

                                                                     cTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTACGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGa                                                                          <  1:973362/143‑1 (MQ=255)

                                                                      ttCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTACGGTGGTAAGGTAGGCTGTTATGGTTGCTGGAGTACGGGGAATTAATCTTTGCTTATGGGTGCAACTAAGGGAt                                                                         >  1:442698/1‑143 (MQ=255)

                                                                          aCATGACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                                                                     <  2:549685/143‑1 (MQ=255)

                                                                          aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTAAGGGTTGGCTTTTATGTTTGCTGGATTAAGGGGACTTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                                                                     >  1:89283/1‑143 (MQ=255)

                                                                          aCATCACGGTGGTTCAATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACgg                                                                     <  2:481300/143‑1 (MQ=255)

                                                                             tCACGGTGTTCCCCTCGAACAAAGTGCGCATATGCTCGCAAAGCCAC‑ATT‑AAGCCAACCTCAAGTTCGTCAGGTTGGCTTCTATCTTTGCTGGATTAAGGGGAATTCCTCCTTGCTTCTTGGTGCTGCTAAGGGATACgcggc                                                                  >  2:418877/1‑139 (MQ=255)

                                                                                 ggTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                              <  1:241219/143‑1 (MQ=255)

                                                                                  gtgtTTCTATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTCAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                             <  2:1034441/143‑1 (MQ=255)

                                                                                  gtgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                             <  1:95275/143‑1 (MQ=255)

                                                                                  gtgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGCTGGCTTTTATGTTTGCTGGGTTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTAt                                                             >  2:1255773/1‑143 (MQ=255)

                                                                                    gtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCGTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGt                                                           >  1:17636/1‑143 (MQ=255)

                                                                                     ctgccATCAAACAAAGTGCGCATATGCGCGCAAAGCAAA‑ATT‑CAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa                                                          <  1:1165526/140‑1 (MQ=255)

                                                                                      ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAg                                                         >  2:182606/1‑143 (MQ=255)

                                                                                       tCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATAAGGGGATTTATGTAga                                                        >  2:307317/1‑142 (MQ=255)

                                                                                         cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCg                                                      >  1:329860/1‑143 (MQ=255)

                                                                                         cATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCg                                                      >  2:1014500/1‑143 (MQ=255)

                                                                                          aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCgg                                                     >  2:60212/1‑143 (MQ=255)

                                                                                          aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCgg                                                     <  1:8703/143‑1 (MQ=255)

                                                                                          aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCgg                                                     <  1:265874/143‑1 (MQ=255)

                                                                                             gAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTg                                                  <  1:356950/143‑1 (MQ=255)

                                                                                                    gTGCGCATATGCTCACAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGAGACGGGCGTTTATGTAGCCGGTTGTAATACa                                           >  2:1192756/1‑143 (MQ=255)

                                                                                                     tGCGCATATGTTCTCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAg                                          <  1:1143247/143‑1 (MQ=255)

                                                                                                     tGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAATGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAg                                          <  2:1117986/143‑1 (MQ=255)

                                                                                                     tGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAg                                          >  1:23573/1‑143 (MQ=255)

                                                                                                       cgtatatGGTCGCAAAGCAGA‑GTT‑AAGCTGACCTTAAGTTCTTAAGGTTGGCTTGTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgcg                                        <  2:863520/140‑1 (MQ=255)

                                                                                                       cgcATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgcg                                        <  2:225293/143‑1 (MQ=255)

                                                                                                           tatGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATa                                    >  1:1114125/1‑143 (MQ=255)

                                                                                                                tCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATgg                               >  2:285886/1‑143 (MQ=255)

                                                                                                                tCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGGTGGCGTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCGGCGCATACATgg                               >  1:256784/1‑143 (MQ=255)

                                                                                                                 cGCAAAGGAAA‑ATT‑ATGCCAACCGTATGTGCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAATGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGc                              <  2:571992/143‑1 (MQ=255)

                                                                                                                  gcaaagcaaa‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCATTGCTGATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGGTGTCATTCAGCGGATACATGGCg                             >  2:82015/1‑143 (MQ=255)

                                                                                                                   caaagcaaa‑ATT‑AAGCCAACCTTAAGTTCTTAAGTTTGGCTTTTATGTTTTCTTGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCgg                            <  2:1036404/143‑1 (MQ=255)

                                                                                                                    aaagcaaa‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGt                           >  1:885814/1‑143 (MQ=255)

                                                                                                                    aaagcaaa‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCAGCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGt                           >  2:38125/1‑143 (MQ=255)

                                                                                                                     aagcaaa‑ATT‑AAGCCCACCTTAAGTGCTTAAGGTGGGCTTCTATGTTTGCTGGATTAAGGGGAATTCATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGTCGGTTGTAATTCAGCGCATACATGGCGGtg                          <  2:285624/143‑1 (MQ=255)

                                                                                                                        caaa‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTCGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTACTTCAGCCCATACCGGGCGGtgctg                       >  2:1212001/1‑143 (MQ=255)

                                                                                                                         aaa‑aTT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                      >  1:168805/1‑143 (MQ=255)

                                                                                                                         aaa‑aTT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                      <  1:1266935/143‑1 (MQ=255)

                                                                                                                           a‑aTT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTACTTCGGCGCATACATGGCGGTGTTGAcc                    >  2:1275951/1‑143 (MQ=255)

                                                                                                                             aTT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCAGTGATGGCCa                   >  1:816804/1‑143 (MQ=255)

                                                                                                                                 aaGCCAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGccccgca                >  2:1078762/1‑139 (MQ=255)

                                                                                                                                   gCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGGGGTGATGAAGATCAgg              >  2:638011/1‑143 (MQ=255)

                                                                                                                                   gCCAACCTCAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCGTTGCTGATTGGTGCATCTAAGGGATACGGGTTTTGATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAgg              <  2:751958/143‑1 (MQ=255)

                                                                                                                                    ccAACCTTTAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCAGCTAAGGGATACGGGTTTTGATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAggg             <  2:917174/143‑1 (MQ=255)

                                                                                                                                    ccAACCTTAAGTTCTTAAGGTTGGCCTTTATGTTTGTTTGAGTAAGGGGAATTAATCTTTGCGTATTGGTGCAGCTAAGGGAGACGGGGTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGAGCATCAggg             >  2:774588/1‑143 (MQ=255)

                                                                                                                                     cAACCTTAAGTTCTTAAGGTTCGCTTCTCCCCTTGATGCGGTGGGAGGAATTAATCTTTGCTTATTTGTGCACCTGAGGGATACGGGAGCTTATGTAGCCGGTTGTAATTCAGCGCATACAGGGCGGGGATGGCCATGAGGGt            >  2:153119/1‑143 (MQ=255)

                                                                                                                                     cAACCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt            >  2:599278/1‑143 (MQ=255)

                                                                                                                                       aCCTTGAGTTCTTAAGGTTGGATTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGa          <  1:492166/143‑1 (MQ=255)

                                                                                                                                       aCCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGa          <  2:380845/143‑1 (MQ=255)

                                                                                                                                           tAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGAGTAAGGGGAATTAATCTTTGCTTATTGGTGGATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa      <  2:559942/143‑1 (MQ=255)

                                                                                                                                              gTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGa   <  1:947816/143‑1 (MQ=255)

                                                                                                                                               ttCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATATGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAc  <  1:1046129/143‑1 (MQ=255)

GATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAA‑ATT‑AAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAC  >  NC_000913/4335513‑4335796

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