BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,748,385	A→G	D115D (GAT→GAC)	grpE ←	GrpE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,748,385	0	A	G	99.6%	938.8 / ‑4.6	252	D115D (GAT→GAC)	grpE	GrpE
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (146/105); total (147/105)
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.29e-01

Only 100 of 255 total aligned reads displayed.

AGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAA  >  NC_000913/2748295‑2748474

aGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTc                                                                                         <  2:5365271/93‑1 (MQ=255)

 gTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTcc                                                                                        >  1:2331001/1‑93 (MQ=255)

 gTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTcc                                                                                        >  2:1822419/1‑93 (MQ=255)

  tttCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCg                                                                                       <  1:1255699/93‑1 (MQ=255)

  tttCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCg                                                                                       <  2:4137518/93‑1 (MQ=255)

    tCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCggg                                                                                     >  2:1069559/1‑93 (MQ=255)

     cGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGt                                                                                    >  1:2178847/1‑93 (MQ=255)

       gCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTa                                                                                  <  1:1762052/93‑1 (MQ=255)

       gCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTa                                                                                  <  1:2036195/93‑1 (MQ=255)

        cGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAg                                                                                 <  1:290490/93‑1 (MQ=255)

         gATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGc                                                                                >  1:4261587/1‑93 (MQ=255)

         gATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGc                                                                                <  1:4964411/93‑1 (MQ=255)

         gATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGc                                                                                <  2:6370016/93‑1 (MQ=255)

           tCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCtt                                                                              >  2:6595355/1‑93 (MQ=255)

            cACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCttt                                                                             <  2:5491678/93‑1 (MQ=255)

            cACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCttt                                                                             >  2:4018016/1‑93 (MQ=255)

             aCTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTa                                                                            <  1:2853777/93‑1 (MQ=255)

              cTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTAt                                                                           <  2:2212565/93‑1 (MQ=255)

               ttCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATc                                                                          <  2:2364243/93‑1 (MQ=255)

                  aaCGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGc                                                                       <  1:692227/93‑1 (MQ=255)

                   aCGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGcc                                                                      >  2:6220388/1‑93 (MQ=255)

                   aCGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGcc                                                                      <  2:6139468/93‑1 (MQ=255)

                   aCGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGcc                                                                      >  1:5239093/1‑93 (MQ=255)

                    cGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCa                                                                     <  1:2510780/93‑1 (MQ=255)

                    cGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCa                                                                     <  2:342843/93‑1 (MQ=255)

                     gCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCAc                                                                    >  2:839805/1‑93 (MQ=255)

                      ccAAACTTACGCACACCATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACt                                                                   <  2:44018/93‑1 (MQ=255)

                      ccAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACt                                                                   <  2:4449753/93‑1 (MQ=255)

                      ccAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACt                                                                   <  1:991836/93‑1 (MQ=255)

                        aaaCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTc                                                                 >  1:2125257/1‑93 (MQ=255)

                         aaCTTACGCACAACATCCAGCATCGACTTCAGGGGCAGCTCAATGCCTTCAACCATCGCAGACAGGTCCGGGTTAGCTTTATCAGCCACTTcc                                                                <  1:418906/93‑1 (MQ=255)

                         aaCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTcc                                                                >  2:2716574/1‑93 (MQ=255)

                         aaCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTcc                                                                >  1:2629937/1‑93 (MQ=255)

                         aaCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTcc                                                                <  1:5434680/93‑1 (MQ=255)

                          aCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCa                                                               <  2:3946942/93‑1 (MQ=255)

                          aCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCa                                                               <  1:789306/93‑1 (MQ=255)

                          aCTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCa                                                               >  1:503604/1‑93 (MQ=255)

                           cTTACGCACAACGTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGCTAGCTTTATCAGCCACTTCCAg                                                              >  1:5657131/1‑93 (MQ=255)

                           cTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAg                                                              <  1:877269/93‑1 (MQ=255)

                           cTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAg                                                              >  2:2140072/1‑93 (MQ=255)

                            ttACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAgc                                                             <  1:4882722/93‑1 (MQ=255)

                            ttACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAgc                                                             >  2:3676841/1‑93 (MQ=255)

                                gCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCAc                                                         >  2:5093256/1‑93 (MQ=255)

                                 cacaACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACg                                                        <  1:1324393/93‑1 (MQ=255)

                                  acaacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGa                                                       >  2:6434037/1‑93 (MQ=255)

                                  acaacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGa                                                       <  2:1591343/93‑1 (MQ=255)

                                  acaacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGa                                                       <  1:1545878/93‑1 (MQ=255)

                                   caacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGAt                                                      <  2:6490075/93‑1 (MQ=255)

                                    aacaTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATc                                                     >  1:6529922/1‑93 (MQ=255)

                                       aTCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAg                                                  <  2:6100745/93‑1 (MQ=255)

                                         ccAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGc                                                <  2:958035/93‑1 (MQ=255)

                                           agaaTCAACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCTCACGATCCAGGCTa                                              <  2:1224574/90‑1 (MQ=255)

                                           aGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTa                                              >  1:2572866/1‑93 (MQ=255)

                                           aGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTa                                              >  2:4607074/1‑93 (MQ=255)

                                              aTCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatca                                           <  1:6508007/93‑1 (MQ=255)

                                              aTCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatca                                           <  1:138913/93‑1 (MQ=255)

                                              aTCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatca                                           >  2:2436633/1‑93 (MQ=255)

                                               tCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaa                                          >  1:4411161/1‑93 (MQ=255)

                                               tCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaa                                          >  1:1871491/1‑93 (MQ=255)

                                               tCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaa                                          >  1:4467997/1‑93 (MQ=255)

                                               tCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaa                                          >  1:57346/1‑93 (MQ=255)

                                                cGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaat                                         <  2:1139340/93‑1 (MQ=255)

                                                cGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaat                                         <  1:4858578/93‑1 (MQ=255)

                                                 gACTTCAGCGACAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaatc                                        <  1:3756812/93‑1 (MQ=255)

                                                 gACTTCAACGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCATGATCCAGGCTatcaatc                                        >  1:6311726/1‑93 (MQ=255)

                                                  aCTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGCCATGTCCGGGTTAGCTTTATCAGCCCCTTCCACCGCACGATCCAGCCTatcaatac                                       >  2:716445/1‑91 (MQ=255)

                                                  aCTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTatcaatca                                       >  2:3216657/1‑93 (MQ=255)

                                                   cTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCAc                                      <  2:5357335/93‑1 (MQ=255)

                                                   cTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCAc                                      >  1:3533332/1‑93 (MQ=255)

                                                         cGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGgcag                                >  2:3814563/1‑93 (MQ=255)

                                                         cGTCAGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGGTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGTCTATCACGCACCGgcgt                                >  1:6452102/1‑91 (MQ=255)

                                                             aGCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAAt                            >  2:4091655/1‑93 (MQ=255)

                                                              gCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAAtt                           >  2:6049182/1‑93 (MQ=255)

                                                              gCTCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAAtt                           >  2:5068898/1‑93 (MQ=255)

                                                                tCAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCg                         >  1:2090598/1‑93 (MQ=255)

                                                                 cAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGt                        <  2:4192415/93‑1 (MQ=255)

                                                                 cAATGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGt                        <  2:5172930/93‑1 (MQ=255)

                                                                  aaTGCCTTCAACCATCGCAGACATGTCCGGGTTCGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGtt                       <  1:6104094/93‑1 (MQ=255)

                                                                  aaTGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGtt                       >  2:2538856/1‑93 (MQ=255)

                                                                   aTGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtg                      >  1:1484486/1‑93 (MQ=255)

                                                                    tGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtga                     <  1:4169025/93‑1 (MQ=255)

                                                                    tGCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtga                     <  2:3115417/93‑1 (MQ=255)

                                                                     gCCTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtgat                    >  1:838484/1‑93 (MQ=255)

                                                                       cTTCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTtgatga                  >  1:6027349/1‑93 (MQ=255)

                                                                        ttCAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGaa                 >  1:3393311/1‑93 (MQ=255)

                                                                          cAACCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAAtt               >  1:270816/1‑93 (MQ=255)

                                                                           aaCCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAAttt              >  1:4949136/1‑93 (MQ=255)

                                                                           aaCCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAAttt              >  1:1132892/1‑93 (MQ=255)

                                                                           aaCCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAAttt              <  1:2240924/93‑1 (MQ=255)

                                                                            aCCATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTtc             >  1:5569092/1‑93 (MQ=255)

                                                                              cATCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTtctc           <  2:3395020/93‑1 (MQ=255)

                                                                               aTCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTcc          <  1:292844/93‑1 (MQ=255)

                                                                               aTCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTcc          <  2:4950441/93‑1 (MQ=255)

                                                                                tCGCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCa         <  2:1346412/93‑1 (MQ=255)

                                                                                  gCAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAgc       <  1:5002448/93‑1 (MQ=255)

                                                                                   cAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAgcg      >  1:2354119/1‑93 (MQ=255)

                                                                                   cAGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAgcg      <  1:4750406/93‑1 (MQ=255)

                                                                                    aGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAgcgc     <  2:4825948/93‑1 (MQ=255)

                                                                                    aGACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAgcgc     <  2:2409506/93‑1 (MQ=255)

                                                                                     gACATGTCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGcgcg    >  2:5936642/1‑93 (MQ=255)

AGTTTCGGCGATCACTTCAACGCCAAACTTACGCACAACATCCAGCATCGACTTCAGCGTCAGCTCAATGCCTTCAACCATCGCAGACATATCCGGGTTAGCTTTATCAGCCACTTCCAGCGCACGATCCAGGCTATCAATCACCGGCAGCAATTCGTTGATGAATTTCTCCAGCGCGAA  >  NC_000913/2748295‑2748474

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