BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,007,515	T→C	V5A (GTT→GCT)	yedF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,007,515	0	T	C	99.2%	415.7 / NA	133	V5A (GTT→GCT)	yedF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (77/55); minor base A (1/0); total (78/55)
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.00e+00

Only 100 of 139 total aligned reads displayed.

CTGCGCTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGTTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAGGGG  >  NC_000913/2007426‑2007604

ctgcGCTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTcc                                                                                        >  2:1973312/1‑93 (MQ=255)

 tgcGCTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCt                                                                                       <  1:2504575/93‑1 (MQ=255)

   cgcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGa                                                                                     <  1:2063067/93‑1 (MQ=255)

   cgcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGa                                                                                     <  2:2828432/93‑1 (MQ=255)

   cgcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGa                                                                                     <  2:1085771/93‑1 (MQ=255)

   cgcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGa                                                                                     <  2:437821/93‑1 (MQ=255)

    gcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGAt                                                                                    >  2:289604/1‑93 (MQ=255)

    gcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGAt                                                                                    <  2:1983500/93‑1 (MQ=255)

    gcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGAt                                                                                    >  1:1834071/1‑93 (MQ=255)

    gcTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTACGCCGTGCGGCCGCGCAGACTGCTAAGGAGATCGCATGGACAATATCGCTCCTGAt                                                                                    >  2:746116/1‑93 (MQ=255)

      tGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTa                                                                                  >  2:1059016/1‑93 (MQ=255)

      tGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTa                                                                                  >  1:4067410/1‑93 (MQ=255)

       gATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTAc                                                                                 >  1:3884137/1‑93 (MQ=255)

       gATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTAc                                                                                 >  1:412430/1‑93 (MQ=255)

        aggTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTAcc                                                                                <  2:4127113/91‑1 (MQ=255)

         tGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCg                                                                               >  2:2460833/1‑93 (MQ=255)

          gTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGc                                                                              >  1:4332415/1‑93 (MQ=255)

           ttGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGcc                                                                             >  1:2225986/1‑93 (MQ=255)

            tGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCt                                                                            <  2:175024/93‑1 (MQ=255)

              aTTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTg                                                                          >  1:1002037/1‑93 (MQ=255)

              aTTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTg                                                                          >  1:915463/1‑93 (MQ=255)

               ttATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGa                                                                         >  2:4070778/1‑93 (MQ=255)

               ttATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGa                                                                         <  2:3601711/93‑1 (MQ=255)

                 aTCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGata                                                                       >  2:2073594/1‑93 (MQ=255)

                  tCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGatat                                                                      <  2:1806251/93‑1 (MQ=255)

                     gCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGCAAAATATCGCTCCTGATTACCGCCTTGATAtggg                                                                   >  1:2599572/1‑92 (MQ=255)

                      cTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATAtggtg                                                                  <  2:941696/93‑1 (MQ=255)

                         ggAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATggtgggt                                                               <  1:3086894/93‑1 (MQ=255)

                         ggAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATggtgggg                                                               >  2:4550861/1‑92 (MQ=255)

                          gAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATggtgggtg                                                              >  1:3639212/1‑93 (MQ=255)

                            aaaaCGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGaa                                                            >  2:2959880/1‑93 (MQ=255)

                            aaaaCGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGaa                                                            >  2:4381578/1‑93 (MQ=255)

                            aaaaCGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTCaa                                                            <  2:650597/93‑1 (MQ=255)

                            aaaaCGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGGGaa                                                            >  2:4255604/1‑93 (MQ=255)

                            aaaaCGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTACTGATTACAGCCTTGATATGGTGGGTGaa                                                            >  1:2006939/1‑93 (MQ=255)

                              aaCGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAAcc                                                          >  1:1157299/1‑93 (MQ=255)

                              aaCGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGGGAAcc                                                          >  1:1289386/1‑93 (MQ=255)

                                cGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCt                                                        >  1:3366242/1‑93 (MQ=255)

                                 gCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTg                                                       >  2:268139/1‑93 (MQ=255)

                                     cttcCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCt                                                   <  1:4240672/93‑1 (MQ=255)

                                          gccgTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATcc                                              >  1:3603948/1‑93 (MQ=255)

                                           ccgTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCg                                             <  1:4584374/93‑1 (MQ=255)

                                                 ggCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTc                                       <  1:4467073/93‑1 (MQ=255)

                                                 ggCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTc                                       <  1:1563332/93‑1 (MQ=255)

                                                   cACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGc                                     >  1:516567/1‑93 (MQ=255)

                                                     ccGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCaa                                   <  2:1965600/93‑1 (MQ=255)

                                                     ccGCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCaa                                   <  2:3372127/93‑1 (MQ=255)

                                                      cGCAGACTGCTAAGGAGATCGCATGAAAAATATTGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAAc                                  >  1:1899040/1‑93 (MQ=255)

                                                       gCAGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAAcc                                 <  1:740114/93‑1 (MQ=255)

                                                         aGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCt                               >  1:2783238/1‑93 (MQ=255)

                                                         aGACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCt                               <  2:1694394/93‑1 (MQ=255)

                                                          gACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCtt                              >  2:950346/1‑93 (MQ=255)

                                                          gACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCtt                              <  2:3914758/93‑1 (MQ=255)

                                                          gACTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCtt                              >  1:3878679/1‑93 (MQ=255)

                                                          gACTGCTAAGGAGATCGCATGAAAAATATAGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCtt                              >  1:1107979/1‑93 (MQ=255)

                                                            cTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCGACCCTTGa                            <  1:3816266/93‑1 (MQ=255)

                                                            cTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGa                            <  2:128518/93‑1 (MQ=255)

                                                            cTGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCCCCTTGATTTGGTGGGTGAACCCTGCCCTTATCGGGACGCCCTAACCCTTgg                            >  2:1259620/1‑92 (MQ=255)

                                                             tGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAg                           >  2:3870783/1‑93 (MQ=255)

                                                             tGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAg                           <  2:3652873/93‑1 (MQ=255)

                                                             tGCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAg                           >  1:427134/1‑93 (MQ=255)

                                                              gCTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAgg                          >  2:2250591/1‑93 (MQ=255)

                                                               cTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGc                         <  2:2717321/93‑1 (MQ=255)

                                                               cTAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGc                         <  1:41823/93‑1 (MQ=255)

                                                                tAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCg                        >  2:121110/1‑93 (MQ=255)

                                                                tAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCg                        >  2:2922230/1‑93 (MQ=255)

                                                                tAAGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCg                        >  2:4425228/1‑93 (MQ=255)

                                                                 aaGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGa                       <  2:3543304/93‑1 (MQ=255)

                                                                  aGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGAt                      <  1:4125007/93‑1 (MQ=255)

                                                                  aGGAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGAt                      >  2:2913117/1‑93 (MQ=255)

                                                                   ggAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATg                     >  2:2669555/1‑93 (MQ=255)

                                                                   ggAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATg                     >  2:261483/1‑93 (MQ=255)

                                                                   ggAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATg                     >  1:1893856/1‑93 (MQ=255)

                                                                   ggAGATCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATg                     >  2:1766581/1‑93 (MQ=255)

                                                                     agaTCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGcc                   <  1:1608343/93‑1 (MQ=255)

                                                                      gaTCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCg                  >  1:1118474/1‑93 (MQ=255)

                                                                       aTCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGc                 >  1:3224085/1‑93 (MQ=255)

                                                                       aTCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGc                 <  2:1094984/93‑1 (MQ=255)

                                                                        tCGCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCa                <  2:4349587/93‑1 (MQ=255)

                                                                          gCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGt              <  1:2638023/93‑1 (MQ=255)

                                                                          gCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGt              <  1:1633197/93‑1 (MQ=255)

                                                                          gCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGt              >  1:1401088/1‑93 (MQ=255)

                                                                          gCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGt              >  1:1794218/1‑93 (MQ=255)

                                                                          gCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGt              >  1:1168452/1‑93 (MQ=255)

                                                                          gCATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGt              <  1:3113760/93‑1 (MQ=255)

                                                                           cATGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGtt             >  1:1458746/1‑93 (MQ=255)

                                                                            aTGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGGCGCAGGTa            >  1:469829/1‑93 (MQ=255)

                                                                            aTGAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTa            >  2:1305930/1‑93 (MQ=255)

                                                                              gAAAAATATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTaaa          >  2:2854514/1‑93 (MQ=255)

                                                                                aaaaTATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTaaaaa        >  1:2434168/1‑93 (MQ=255)

                                                                                aaaaTATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTaaaaa        >  1:72358/1‑93 (MQ=255)

                                                                                 aaaTATCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCATCCCTTGAGGCGATGCCGCAGTTaaaaaa       <  1:2863034/93‑1 (MQ=255)

                                                                                  aaTATCGCTCCTGATTTCCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTaaaaaaa      <  2:3168518/93‑1 (MQ=255)

                                                                                  aaTATCGCTCCTGATTACCGCCTTGATATGGTGGGTGACCCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTaaaaaaa      <  2:4229569/93‑1 (MQ=255)

                                                                                    tatCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAgg    >  1:1740100/1‑93 (MQ=255)

                                                                                    tatCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAgg    <  2:70002/93‑1 (MQ=255)

                                                                                    tatCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAgg    >  2:161881/1‑93 (MQ=255)

                                                                                     atCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAggg   <  2:4325702/93‑1 (MQ=255)

                                                                                      tCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAgggg  <  1:2575485/93‑1 (MQ=255)

                                                                                      tCGCTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAgggg  <  1:3091034/93‑1 (MQ=255)

CTGCGCTGATGTTGATTATCGGCTGGGAAAAACGCTTCTTCCGCCGTGCGGCACCGCAGACTGCTAAGGAGATCGCATGAAAAATATCGTTCCTGATTACCGCCTTGATATGGTGGGTGAACCCTGCCCTTATCCGGCAGTCGCAACCCTTGAGGCGATGCCGCAGTTAAAAAAAGGGG  >  NC_000913/2007426‑2007604

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