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	.	99.3%	665.6 / ‑2.7	146	coding (25/1185 nt)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (76/69); minor base T (1/0); total (77/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

Only 100 of 163 total aligned reads displayed.

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCA  >  NC_000913/2487111‑2487291

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

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

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

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

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

   aGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaa                                                                                       >  2:3994591/1‑93 (MQ=255)

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

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

      tCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagc                                                                                    <  2:1565190/93‑3 (MQ=255)

      tCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagc                                                                                    <  1:3630507/93‑3 (MQ=255)

      tCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaagc                                                                                    >  1:2029003/1‑91 (MQ=255)

         aTTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGa                                                                                >  2:20834/1‑93 (MQ=255)

         aTTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGa                                                                                >  1:2578636/1‑93 (MQ=255)

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

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

            aaGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAgg                                                                             <  1:4089335/93‑1 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAggg                                                                            >  2:46435/1‑92 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                            >  2:1019068/1‑93 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                            >  2:3651252/1‑93 (MQ=255)

             aGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGt                                                                            >  1:286371/1‑93 (MQ=255)

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

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

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

                   tgtgTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAg                                                                      <  1:1232236/93‑1 (MQ=255)

                     tgtCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGac                                                                    >  1:1048353/1‑93 (MQ=255)

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

                       tCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGacac                                                                  <  1:1342989/93‑1 (MQ=255)

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

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

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

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

                            tCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCgg                                                             >  1:3498514/1‑93 (MQ=255)

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

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

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTa                                                       >  2:1822315/1‑93 (MQ=255)

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTa                                                       >  1:275699/1‑93 (MQ=255)

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTa                                                       <  1:1012231/93‑1 (MQ=255)

                                  aCTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTa                                                       >  1:4196581/1‑93 (MQ=255)

                                       taaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaa                                                  <  1:1307712/93‑1 (MQ=255)

                                         aGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTaaaa                                                >  2:740886/1‑93 (MQ=255)

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

                                           ccTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATc                                              <  1:4487577/93‑1 (MQ=255)

                                           ccTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATc                                              <  1:888140/93‑1 (MQ=255)

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

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

                                                aaaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                         <  1:1443053/93‑1 (MQ=255)

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

                                                 aaCGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAAc                                        >  1:3837409/1‑93 (MQ=255)

                                                  aCGTCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACa                                       <  2:3873340/93‑1 (MQ=255)

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

                                                     gCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGa                                    <  2:3338760/93‑1 (MQ=255)

                                                      cccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGat                                   >  2:2815301/1‑93 (MQ=255)

                                                      cccTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGat                                   >  1:3608084/1‑93 (MQ=255)

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

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

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                 >  2:806553/1‑93 (MQ=255)

                                                        cTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGatta                                 >  1:3086006/1‑93 (MQ=255)

                                                         tttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                >  1:792909/1‑93 (MQ=255)

                                                         tttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                <  1:2429676/93‑1 (MQ=255)

                                                          ttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCGGGAGGTGAAAAGACACTTCGGCCTTTATGTAAAATCAATAACATGattaat                               >  1:2987834/1‑93 (MQ=255)

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

                                                          ttGCTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaat                               >  1:4682847/1‑93 (MQ=255)

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

                                                            tcTTTCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatta                             <  1:3864275/92‑1 (MQ=255)

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

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

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

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

                                                              tttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATg                           <  2:3600952/93‑1 (MQ=255)

                                                               ttCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGa                          >  2:15323/1‑93 (MQ=255)

                                                                tCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaa                         <  1:4618279/93‑1 (MQ=255)

                                                                tCATTATTTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaa                         <  1:3871329/93‑1 (MQ=255)

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

                                                                  attattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAc                       >  1:4264058/1‑93 (MQ=255)

                                                                   ttattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAcc                      >  1:720628/1‑93 (MQ=255)

                                                                   ttattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAcc                      <  1:3412510/93‑1 (MQ=255)

                                                                   ttattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAAca                      >  2:1082586/1‑92 (MQ=255)

                                                                    tattTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCa                     >  2:3269981/1‑93 (MQ=255)

                                                                     attTGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAg                    >  2:2964497/1‑93 (MQ=255)

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

                                                                      tttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGt                   <  1:3282397/93‑1 (MQ=255)

                                                                       ttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTc                  >  2:2101956/1‑93 (MQ=255)

                                                                       ttGTCATTTGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTc                  <  2:3842448/93‑1 (MQ=255)

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

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

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

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

                                                                             tttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGa            >  2:2467419/1‑93 (MQ=255)

                                                                             tttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGa            >  2:4435516/1‑93 (MQ=255)

                                                                              ttGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGAc           <  1:3228984/93‑1 (MQ=255)

                                                                               tGCTTATGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACa          <  2:2605086/93‑1 (MQ=255)

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

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

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

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

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

                                                                                    aTGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAAt     >  2:3215836/1‑93 (MQ=255)

                                                                                     tGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATc    <  2:234989/93‑1 (MQ=255)

                                                                                     tGCCTTG‑AAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATc    >  1:1690664/1‑93 (MQ=255)

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

GAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCA  >  NC_000913/2487111‑2487291

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