BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	27,903	T→C	L204P (CTG→CCG)	yaaF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	27,903	0	T	C	99.5%	570.5 / ‑5.1	181	L204P (CTG→CCG)	yaaF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (97/83); total (97/84)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 188 total aligned reads displayed.

CCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACG  >  NC_000913/27813‑27992

 cTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGcc                                                                                        >  2:4400174/1‑93 (MQ=255)

 cTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGcc                                                                                        <  1:1009975/93‑1 (MQ=255)

  tgtgTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACAccgccg                                                                                       >  2:7562/1‑93 (MQ=255)

  tgtgTCTTCCGCAGTGGCATTGAAATCGTCATGTGCGGGTTGGATGTCACCAATCAGGCAATATTACTTCCTGACTATCTCTCTACACCGcct                                                                                       >  2:4010879/1‑92 (MQ=255)

   gtgtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACcgccgc                                                                                      >  1:310573/1‑93 (MQ=255)

    tgtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCa                                                                                     >  2:3165854/1‑93 (MQ=255)

     gtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAg                                                                                    >  2:1115139/1‑93 (MQ=255)

     gtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAg                                                                                    >  1:58627/1‑93 (MQ=255)

     gtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAg                                                                                    >  1:4298488/1‑93 (MQ=255)

      tCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGt                                                                                   <  1:752910/93‑1 (MQ=255)

      tCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGt                                                                                   <  2:1012983/93‑1 (MQ=255)

       cTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGtt                                                                                  <  1:1967795/93‑1 (MQ=255)

        ttCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCGCCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTa                                                                                 >  1:2508413/1‑93 (MQ=255)

          ccGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTaaa                                                                               >  1:645400/1‑93 (MQ=255)

          ccGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTaaa                                                                               <  1:3027238/93‑1 (MQ=255)

            gCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAcc                                                                             >  2:1429356/1‑93 (MQ=255)

             cAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCg                                                                            >  2:3461628/1‑93 (MQ=255)

             cAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCg                                                                            >  1:1297214/1‑93 (MQ=255)

              aGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGt                                                                           >  2:3505137/1‑93 (MQ=255)

               gTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTa                                                                          <  2:1991128/93‑1 (MQ=255)

               gTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTa                                                                          <  2:2625643/93‑1 (MQ=255)

                  gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCg                                                                       >  1:2398660/1‑93 (MQ=255)

                  gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCg                                                                       >  1:1145063/1‑93 (MQ=255)

                     ttGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGa                                                                    <  2:488781/93‑1 (MQ=255)

                         aaTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAAt                                                                >  2:3286042/1‑93 (MQ=255)

                          aTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAAtc                                                               >  1:3016848/1‑92 (MQ=255)

                          aTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATg                                                               <  1:789309/93‑1 (MQ=255)

                           tCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGc                                                              >  2:3592284/1‑93 (MQ=255)

                           tCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGc                                                              >  2:3103523/1‑93 (MQ=255)

                           tCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGc                                                              <  1:3604080/93‑1 (MQ=255)

                               cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                          <  2:1949769/93‑1 (MQ=255)

                                aTGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCAc                                                         <  1:2572536/93‑1 (MQ=255)

                                 tgtgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACg                                                        <  1:3490426/93‑1 (MQ=255)

                                    gCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGccc                                                     <  1:2373327/93‑1 (MQ=255)

                                      ggTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTg                                                   >  2:4337860/1‑93 (MQ=255)

                                       gTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGt                                                  >  2:4191768/1‑93 (MQ=255)

                                        tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                 <  1:4466579/93‑1 (MQ=255)

                                        tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                 >  1:3333738/1‑93 (MQ=255)

                                        tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                 >  1:3802203/1‑93 (MQ=255)

                                        tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                 <  1:22291/93‑1 (MQ=255)

                                          tgcaTGTCACCAATCAGGCAATATTAGCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTa                                               <  1:2756613/90‑1 (MQ=255)

                                                   ccAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAcc                                      >  1:3068978/1‑93 (MQ=255)

                                                    cAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCg                                     >  1:4583559/1‑93 (MQ=255)

                                                     aaTCAGGCATTATTAGCTCCTGACTATCTCACGACTCGGCCGAAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGt                                    <  2:2708719/93‑1 (MQ=255)

                                                     aaTCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAGAATGCTTCACGCCCTGTTTAGCCACTACCga                                    >  1:710394/1‑92 (MQ=255)

                                                     aaTCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGt                                    >  1:1225091/1‑93 (MQ=255)

                                                     aaTCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGt                                    >  1:4410633/1‑93 (MQ=255)

                                                      aTCAGTCAATATTAACTCCTGACTATCTCTCTACCCCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTa                                   <  1:1040983/93‑1 (MQ=255)

                                                      aTCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTa                                   >  2:168328/1‑93 (MQ=255)

                                                       tCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAg                                  >  2:3791735/1‑93 (MQ=255)

                                                        cAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGc                                 <  1:2644801/93‑1 (MQ=255)

                                                         aGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCg                                >  2:4156019/1‑93 (MQ=255)

                                                          ggCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCgg                               <  2:2594378/93‑1 (MQ=255)

                                                          ggCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCgg                               >  1:119623/1‑93 (MQ=255)

                                                              atatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGt                           >  1:773627/1‑93 (MQ=255)

                                                               tatTACCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTa                          <  2:3991414/93‑1 (MQ=255)

                                                               tatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTa                          >  1:4271009/1‑93 (MQ=255)

                                                               tatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTa                          >  1:200346/1‑93 (MQ=255)

                                                                atTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTAt                         >  2:2412942/1‑93 (MQ=255)

                                                                atTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTAt                         >  2:2069981/1‑93 (MQ=255)

                                                                atTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGCAGCGGCAGTAt                         <  1:1344397/93‑1 (MQ=255)

                                                                  tAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGc                       <  2:1203318/93‑1 (MQ=255)

                                                                  tAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGc                       >  1:4237314/1‑93 (MQ=255)

                                                                  tAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGc                       <  2:4579075/93‑1 (MQ=255)

                                                                  tAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGc                       <  1:3720132/93‑1 (MQ=255)

                                                                    aCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaa                     >  1:163705/1‑93 (MQ=255)

                                                                    aCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaa                     <  1:1790969/93‑1 (MQ=255)

                                                                    aCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaa                     <  1:506730/93‑1 (MQ=255)

                                                                    aCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaa                     >  1:1070577/1‑93 (MQ=255)

                                                                     cTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                    <  2:2024090/93‑1 (MQ=255)

                                                                     cTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                    <  1:3874928/93‑1 (MQ=255)

                                                                     cTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                    <  2:1440406/93‑1 (MQ=255)

                                                                     cTCCTGACGATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                    <  1:2300051/93‑1 (MQ=255)

                                                                      tCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTCTGCAAAg                   >  2:3374382/1‑93 (MQ=255)

                                                                        cTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCg                 >  1:3937331/1‑93 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGTGGCAGTATGCAAAGCgg                >  1:2889325/1‑93 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                <  1:4053576/93‑1 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                >  2:460114/1‑93 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                >  1:2120765/1‑93 (MQ=255)

                                                                          gACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGc               <  1:2852280/93‑1 (MQ=255)

                                                                          gACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGc               >  2:1308437/1‑93 (MQ=255)

                                                                           aCTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCt              <  2:3155006/93‑1 (MQ=255)

                                                                           aCTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCt              <  2:1864300/93‑1 (MQ=255)

                                                                           aCTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCt              >  2:3844721/1‑93 (MQ=255)

                                                                           aCTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCt              <  1:1105608/93‑1 (MQ=255)

                                                                            cTATCTCTCTACCCCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt             <  2:3666607/93‑1 (MQ=255)

                                                                            cTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt             <  2:1805818/93‑1 (MQ=255)

                                                                            cTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt             <  2:3216311/93‑1 (MQ=255)

                                                                              aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtttc           >  1:3175999/1‑91 (MQ=255)

                                                                              aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtttc           >  2:2148658/1‑91 (MQ=255)

                                                                              aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGc           >  2:1398523/1‑93 (MQ=255)

                                                                              aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGc           <  1:544616/93‑1 (MQ=255)

                                                                                ctctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGa         <  2:192462/93‑1 (MQ=255)

                                                                                 tctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGTGaa        >  2:4269647/1‑93 (MQ=255)

                                                                                 tctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGaa        <  2:4257313/93‑1 (MQ=255)

                                                                                   tctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAGTg      >  2:374172/1‑93 (MQ=255)

                                                                                     tACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATggc    >  1:1411220/1‑91 (MQ=255)

                                                                                     tACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCa    <  1:3047952/93‑1 (MQ=255)

                                                                                      acacCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATTCAc   >  2:592745/1‑93 (MQ=255)

                                                                                      acacCGCCGCAGTTAAACCGTACCGAGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCAc   <  2:1767210/93‑1 (MQ=255)

CCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACG  >  NC_000913/27813‑27992

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