BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,474	T→A	19.6%	L61* (TTG→TAG) ‡	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,474	0	T	A	19.6%	407.4 / 92.9	185	L61* (TTG→TAG)	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base T (79/68); new base A (19/17); total (98/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

Only 100 of 204 total aligned reads displayed.

ACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113376‑4113573

aCATAATGACTCTCTCCTGACGTTAATTCATATTGAAGATTGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                     <  1:2639520/99‑1 (MQ=255)

aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                    <  2:3741756/100‑1 (MQ=255)

aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                    <  2:1048351/100‑1 (MQ=255)

aCATAATGACGCTCACCTGACGTTAATTCAGATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                    <  1:3692464/100‑1 (MQ=255)

 cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                   >  1:1129418/1‑100 (MQ=255)

  aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                  <  1:1952550/100‑1 (MQ=255)

   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 >  2:3592963/1‑100 (MQ=255)

    aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGCTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                                                                                                >  1:22702/1‑100 (MQ=255)

      tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                              <  2:83308/100‑1 (MQ=255)

      tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                              >  1:3391039/1‑100 (MQ=255)

          gCTCACCTGACGTTAAGTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                          <  1:2420147/100‑1 (MQ=255)

            tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatc                                                                                        <  2:3888754/100‑3 (MQ=255)

            tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatc                                                                                        <  1:1032100/100‑3 (MQ=255)

             cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                       <  2:2670604/100‑4 (MQ=38)

               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                     <  2:3206801/100‑1 (MQ=255)

               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  2:2749478/100‑6 (MQ=38)

               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  1:3684702/100‑6 (MQ=38)

               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  2:906825/100‑6 (MQ=38)

               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  1:2052242/100‑6 (MQ=38)

               ccTGACGTCAATTCATATTGATGATGGCTTAAGCGAGTTCTATCTGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCCATTTTTAAAGAATTAgcc                                                                                      >  1:1531662/1‑97 (MQ=255)

                cTGACTTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTCCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                    <  1:1235717/100‑7 (MQ=37)

                cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                    <  2:1352979/100‑7 (MQ=38)

                 tgaacgtAATTCATATTGAGGATGGCTTAAGCGATTTGTATCCGGGTATCTACGCACCTGCAACAGAAGCTATTCTTCGATTGTTGAAGAATAAGTc                                                                                      <  1:3312338/91‑1 (MQ=255)

                 tGACGTTAAGTCATATTGATGGTGGCTTAAGCGAGTTGTATCCGGGTATCTACTGCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGa                                                                                   <  1:356582/100‑1 (MQ=255)

                  gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAt                                                                                  <  1:13572/100‑1 (MQ=255)

                  aactttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGa                                                                                   <  1:4218126/95‑1 (MQ=255)

                   aCGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa                                                                                 <  2:2649213/100‑1 (MQ=255)

                          ttCATATTGATGATGGCTTAAGCCAGTTGTAACCTAGTATCTACATCCCTCCAACAAAACATATTCTTCAATAGTTCAAGAATAAATCAGATAACa                                                                              >  1:3495093/1‑96 (MQ=255)

                           tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                         <  2:2324763/100‑1 (MQ=255)

                           tCATATTGATGACGGCTTAAGCGAGGTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATCGTTGAAGAATAAGTCGGATAACAAGCTg                                                                         <  1:1579149/100‑1 (MQ=255)

                             atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                       >  2:1526874/1‑100 (MQ=255)

                              tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                      >  2:1315716/1‑100 (MQ=255)

                              tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                      >  2:3697210/1‑100 (MQ=255)

                               atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                     <  1:963255/100‑1 (MQ=255)

                                 tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                   >  1:1160554/1‑100 (MQ=255)

                                  gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                  >  2:3421124/1‑100 (MQ=255)

                                  gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                  <  2:2243306/100‑1 (MQ=255)

                                  gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                  <  1:3274220/100‑1 (MQ=255)

                                     gatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGGTATTCATCAAATGTTGAAGAATAAGCCGGATAAAAAACTGTATAAACTGa                                                               >  1:1394611/1‑100 (MQ=255)

                                      atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                              >  2:2515550/1‑100 (MQ=255)

                                         gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                           >  2:86194/1‑100 (MQ=255)

                                          cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                          >  1:880920/1‑100 (MQ=255)

                                            tAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaaa                                                        <  2:3046170/100‑1 (MQ=255)

                                            tAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaaa                                                        >  2:1181201/1‑100 (MQ=255)

                                             aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAat                                                       >  2:1205171/1‑100 (MQ=255)

                                              aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                      >  2:3176670/1‑100 (MQ=255)

                                               gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                     >  2:2122224/1‑100 (MQ=255)

                                               gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                     >  1:3224753/1‑100 (MQ=255)

                                                 gggTTGTATCCGGGTATCTACTTCTCTTTAACAGAATATATTCTTTAATTGTTGAAGAATATGTCGTATAACAAGCTGTATAAACTGACGAAAAATAtt                                                    <  1:2348279/97‑1 (MQ=255)

                                                 gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                   <  2:101765/100‑1 (MQ=255)

                                                 gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGGTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                   <  1:2692061/100‑1 (MQ=255)

                                                  aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                  >  1:331940/1‑100 (MQ=255)

                                                  aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                  <  1:586625/100‑1 (MQ=255)

                                                   gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                 >  1:3119025/1‑100 (MQ=255)

                                                        aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATTGc                                             >  1:4103115/1‑99 (MQ=255)

                                                         tCCGGGTATCTACTTCCCTGCAACAGCAGATATTCTGCAATTGTTGAAGACTAAGTCGGATACCAAGCTGTATAAACTCACGAAAAATATTCAATGGcc                                            >  1:426622/1‑99 (MQ=255)

                                                         tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                           <  2:3003011/100‑1 (MQ=255)

                                                         tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGACCg                                           >  2:790205/1‑100 (MQ=255)

                                                          ccGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGa                                          >  2:707888/1‑100 (MQ=255)

                                                          ccGGGTATCTACTTCCCTGCAACAGAAGAGATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGTCGAAAAATATTCAATGGCCg                                           <  1:1602121/99‑1 (MQ=255)

                                                            gggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                        <  2:2690960/100‑1 (MQ=255)

                                                             ggTATCTACTTCCCTGCAACAGAAGATATTCGTCAATTGTTGCAGCATAAGTCGGATAAAAAGCTGTATACACTGACGAAAAATATTCAATGGCCGAAGa                                       >  1:806885/1‑100 (MQ=255)

                                                                 tCTACTTCCCTGAAACATAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTCTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaa                                   <  1:2959411/100‑1 (MQ=255)

                                                                  cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATACACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                  >  2:434600/1‑100 (MQ=255)

                                                                  cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                  >  2:3177988/1‑100 (MQ=255)

                                                                   tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                 >  1:4159388/1‑100 (MQ=255)

                                                                   tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                 <  1:2761935/100‑1 (MQ=255)

                                                                    aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACt                                >  2:409975/1‑100 (MQ=255)

                                                                     cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                               >  2:1670825/1‑100 (MQ=255)

                                                                     cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                               >  2:3762574/1‑100 (MQ=255)

                                                                       tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                              <  1:4388731/99‑1 (MQ=255)

                                                                        cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACGGCGt                            <  2:258900/100‑1 (MQ=255)

                                                                         ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGATGACAAAACTGCGt                            <  1:2727959/99‑1 (MQ=255)

                                                                         ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                           <  1:1315711/100‑1 (MQ=255)

                                                                          cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAGGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                          >  2:4205153/1‑100 (MQ=255)

                                                                          cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                          >  2:2415709/1‑100 (MQ=255)

                                                                           tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                         >  2:4325597/1‑100 (MQ=255)

                                                                           tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                         >  2:1481987/1‑100 (MQ=255)

                                                                            gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                        <  1:1819411/100‑1 (MQ=255)

                                                                              aaCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                      >  2:2007008/1‑100 (MQ=255)

                                                                              aaCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                      <  2:3633412/100‑1 (MQ=255)

                                                                               aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                     >  2:1822235/1‑100 (MQ=255)

                                                                                cAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                    <  1:1171832/100‑1 (MQ=255)

                                                                                  gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  1:2096953/1‑100 (MQ=255)

                                                                                  cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  2:1067124/6‑100 (MQ=38)

                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  2:1981649/5‑100 (MQ=38)

                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  2:3559952/5‑100 (MQ=38)

                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  2:2138063/5‑100 (MQ=38)

                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  1:3976413/5‑100 (MQ=38)

                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATACAATGGCCGAATACAAAACTGCGTATTGaa                      >  1:3995243/5‑95 (MQ=37)

                                                                                    acttatTCTTCAACAATTGAAGGCTAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  1:2567306/4‑100 (MQ=255)

                                                                                    acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                <  2:4228360/97‑1 (MQ=38)

                                                                                     cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag               >  2:3256119/3‑100 (MQ=255)

                                                                                      atatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgaga              <  2:3494237/100‑1 (MQ=255)

                                                                                            ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc        <  2:2699894/100‑1 (MQ=255)

                                                                                             tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc        <  1:3431644/99‑1 (MQ=255)

                                                                                              cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      >  1:2127554/1‑100 (MQ=255)

                                                                                                aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa    >  2:3496560/1‑100 (MQ=255)

                                                                                                aTTGTTGAAGAATAAGCCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa    >  2:606986/1‑100 (MQ=255)

                                                                                                 ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaaa   <  2:1632617/100‑1 (MQ=255)

ACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113376‑4113573

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