BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,487,208	(A)_6→5	coding (20/1185 nt)	yfdE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,487,203	0	A	.	98.4%	831.1 / 4.9	184	coding (25/1185 nt)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base . (103/78); total (105/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 203 total aligned reads displayed.

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCA  >  NC_000913/2487111‑2487291

gAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGa                                                                                          >  1:587892/1‑93 (MQ=255)

 aaaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaa                                                                                         >  1:3765521/1‑93 (MQ=255)

  aaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaa                                                                                        >  2:4237210/1‑93 (MQ=255)

  aaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaa                                                                                        <  1:2401202/93‑1 (MQ=255)

  aaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaa                                                                                        <  2:2205713/93‑1 (MQ=255)

    ggTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaa                                                                                      >  1:1953497/1‑93 (MQ=255)

    ggTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaa                                                                                      >  1:4643725/1‑93 (MQ=255)

    ggTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaa                                                                                      <  1:4266811/93‑1 (MQ=255)

        cATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagccg                                                                                  >  2:2559501/1‑89 (MQ=255)

        cATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagccg                                                                                  <  2:1993186/93‑5 (MQ=255)

        cATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagccg                                                                                  >  1:1972908/1‑89 (MQ=255)

        cATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagccg                                                                                  >  1:1312362/1‑89 (MQ=255)

          ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGaa                                                                               >  2:292115/1‑93 (MQ=255)

          ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGaa                                                                               <  2:3040236/93‑1 (MQ=255)

           tAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAg                                                                              >  1:4517559/1‑93 (MQ=255)

              ggACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTg                                                                           >  1:159295/1‑93 (MQ=255)

               gACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGa                                                                          >  1:1009796/1‑93 (MQ=255)

               gACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGa                                                                          >  2:4665645/1‑93 (MQ=255)

               gACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGa                                                                          <  1:1884744/93‑1 (MQ=255)

                  aTGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGaaaa                                                                       >  2:1589323/1‑93 (MQ=255)

                  aTGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGaaaa                                                                       >  2:1020616/1‑93 (MQ=255)

                  aTGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGaaaa                                                                       <  2:4422461/93‑1 (MQ=255)

                    gtgtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGa                                                                     >  1:3458709/1‑93 (MQ=255)

                    gtgtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGa                                                                     <  1:3736245/93‑1 (MQ=255)

                      gtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGaca                                                                   <  2:3212673/93‑1 (MQ=255)

                      gtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGaca                                                                   >  2:2428454/1‑93 (MQ=255)

                        cATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACt                                                                 <  1:1396753/93‑1 (MQ=255)

                         atatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACtt                                                                >  2:1559515/1‑93 (MQ=255)

                         atatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACtt                                                                >  2:364005/1‑93 (MQ=255)

                          tatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTc                                                               <  2:1075306/93‑1 (MQ=255)

                          tatCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTc                                                               >  1:4627483/1‑93 (MQ=255)

                           atCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCg                                                              >  2:592925/1‑93 (MQ=255)

                           atCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCg                                                              <  2:1424832/93‑1 (MQ=255)

                            tcaaTAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACCCTTCgg                                                             <  1:1714950/90‑1 (MQ=255)

                             cGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGc                                                            <  1:1454006/93‑1 (MQ=255)

                              gATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCa                                                           <  2:788145/93‑1 (MQ=255)

                               aTAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAt                                                          >  1:302800/1‑93 (MQ=255)

                                tAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAtt                                                         >  2:1689208/1‑93 (MQ=255)

                                tAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAtt                                                         <  2:2837138/93‑1 (MQ=255)

                                 aaCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCAttt                                                        <  1:2205612/93‑1 (MQ=255)

                                   cTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTAt                                                      <  2:647218/93‑1 (MQ=255)

                                   cTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTAt                                                      >  1:398185/1‑93 (MQ=255)

                                   cTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTAt                                                      >  1:918844/1‑93 (MQ=255)

                                    taataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATg                                                     >  1:3581740/1‑93 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                    >  2:3370062/1‑93 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                    >  1:2364007/1‑93 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                    >  2:1296625/1‑93 (MQ=255)

                                     aataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGt                                                    <  1:3959778/93‑1 (MQ=255)

                                      ataaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTa                                                   <  2:3324833/93‑1 (MQ=255)

                                       taaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaa                                                  >  2:2379019/1‑93 (MQ=255)

                                       taaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaa                                                  >  2:1918654/1‑93 (MQ=255)

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaa                                                 >  2:3305596/1‑93 (MQ=255)

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaa                                                 <  2:828162/93‑1 (MQ=255)

                                        aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaa                                                 <  1:480214/93‑1 (MQ=255)

                                          gCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAAt                                               >  1:412322/1‑93 (MQ=255)

                                            cTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGCAGGTGAAAAGACACTTCGGCATTTATGTAAAATCa                                             >  1:4352501/1‑93 (MQ=255)

                                             ttCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCaa                                            >  2:4545949/1‑93 (MQ=255)

                                              tCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAAt                                           <  2:3058300/93‑1 (MQ=255)

                                              tCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAAt                                           >  1:863323/1‑93 (MQ=255)

                                              tCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAAt                                           >  1:2015131/1‑93 (MQ=255)

                                              tCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAAt                                           <  2:650451/93‑1 (MQ=255)

                                               cAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATa                                          <  2:3957115/93‑1 (MQ=255)

                                                aaaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                         >  2:2951883/1‑93 (MQ=255)

                                                aaaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                         >  1:2840251/1‑93 (MQ=255)

                                                 aaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAAc                                        <  1:935639/93‑1 (MQ=255)

                                                  aCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACa                                       >  1:4078471/1‑93 (MQ=255)

                                                   cGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACa                                       <  1:1578478/93‑1 (MQ=255)

                                                   cGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACAt                                      <  2:163116/93‑1 (MQ=255)

                                                   cGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACAt                                      >  2:3278967/1‑93 (MQ=255)

                                                    ggCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATg                                     <  1:325202/93‑1 (MQ=255)

                                                     gCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGa                                    >  2:4138265/1‑93 (MQ=255)

                                                       ccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatt                                  >  2:4578254/1‑93 (MQ=255)

                                                          tttcTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaat                               <  1:1100678/90‑1 (MQ=255)

                                                          ttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaat                               <  1:802734/93‑1 (MQ=255)

                                                            gCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatta                             <  1:4597321/93‑1 (MQ=255)

                                                            gCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatta                             >  2:4186128/1‑93 (MQ=255)

                                                             cTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTAt                            <  2:3289099/93‑1 (MQ=255)

                                                               ttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGa                          <  1:4152126/93‑1 (MQ=255)

                                                               ttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGa                          <  2:4258362/93‑1 (MQ=255)

                                                                 cATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaaa                        <  1:2932064/93‑1 (MQ=255)

                                                                  attattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAc                       <  2:2098751/93‑1 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                     <  1:2081053/93‑1 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                     <  1:2470119/93‑1 (MQ=255)

                                                                     attTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAg                    <  2:3444077/93‑1 (MQ=255)

                                                                      tttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAgc                   >  1:3645309/1‑92 (MQ=255)

                                                                      tttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGt                   >  1:2225885/1‑93 (MQ=255)

                                                                      tttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGt                   >  2:4091706/1‑93 (MQ=255)

                                                                         gTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGa                >  2:2676007/1‑93 (MQ=255)

                                                                          tCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGAc               <  2:4185682/93‑1 (MQ=255)

                                                                            aTTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTg             >  2:3791940/1‑93 (MQ=255)

                                                                            aTTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTg             >  2:3751821/1‑93 (MQ=255)

                                                                             tttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGa            <  1:676361/93‑1 (MQ=255)

                                                                                gCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACAt         >  1:1749940/1‑93 (MQ=255)

                                                                                  ttATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaa       >  1:528331/1‑93 (MQ=255)

                                                                                  ttATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaa       >  1:2831311/1‑93 (MQ=255)

                                                                                  ttATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGAAATaa       <  2:3384218/93‑1 (MQ=255)

                                                                                   tATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaaa      >  2:1127303/1‑93 (MQ=255)

                                                                                   tATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaaa      >  2:170545/1‑93 (MQ=255)

                                                                                   tATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaaa      >  2:2263408/1‑93 (MQ=255)

                                                                                       ccTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCca  >  2:232002/1‑93 (MQ=255)

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCA  >  NC_000913/2487111‑2487291

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