BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,361,323	T→A	Q31L (CAG→CTG)	rhmD ←	L‑rhamnonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,361,323	0	T	A	98.5%	452.1 / ‑3.5	136	Q31L (CAG→CTG)	rhmD	L‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base A (73/61); total (74/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 136 total aligned reads displayed.

AAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGT  >  NC_000913/2361226‑2361419

aaGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAgg                                                                                                <  1:2529158/100‑1 (MQ=255)

aaGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAgg                                                                                                <  1:2249885/100‑1 (MQ=255)

 aGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGt                                                                                               >  1:1955080/1‑100 (MQ=255)

 aGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATAGCGGCATTTACTCATCGGGGTGGCAATATGACCGTCTATCCAGTGATTCGCCCCCAgg                                                                                                >  1:1294137/1‑99 (MQ=255)

    aCGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGt                                                                                            <  1:2388414/100‑1 (MQ=255)

     cGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGt                                                                                            <  2:1400149/99‑1 (MQ=255)

      gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                          >  2:1455838/1‑100 (MQ=255)

      gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                          <  1:1452774/100‑1 (MQ=255)

      gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                          >  2:597171/1‑100 (MQ=255)

      gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                          <  1:148403/100‑1 (MQ=255)

         aaTGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATaa                                                                                       <  2:1122735/100‑1 (MQ=255)

         aaTGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATaa                                                                                       <  2:280936/100‑1 (MQ=255)

          aTGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAAt                                                                                      <  2:953536/100‑1 (MQ=255)

           tGCCAAATGACTTGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCTAGTGATTTGCTCCCAGGTCGTGATAAt                                                                                      >  1:1525255/1‑99 (MQ=255)

           tGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATc                                                                                     <  2:84627/100‑1 (MQ=255)

            gCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCa                                                                                    <  2:1239437/100‑1 (MQ=255)

            gCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCa                                                                                    <  2:331046/100‑1 (MQ=255)

               aaaTGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCAccg                                                                                 >  2:541828/1‑100 (MQ=255)

                     cTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCa                                                                           >  1:1844220/1‑100 (MQ=255)

                      tGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAg                                                                          >  1:3082800/1‑100 (MQ=255)

                      tGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAg                                                                          <  2:727001/100‑1 (MQ=255)

                       ggCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAgc                                                                         >  2:1588794/1‑100 (MQ=255)

                       ggCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAgc                                                                         >  2:2978171/1‑100 (MQ=255)

                         cGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAgcgc                                                                       >  2:1173839/1‑100 (MQ=255)

                         cGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAgcgc                                                                       >  1:2475702/1‑100 (MQ=255)

                           tGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCt                                                                     >  1:1602141/1‑100 (MQ=255)

                              cTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCtttt                                                                  >  2:702796/1‑100 (MQ=255)

                               tGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCttttt                                                                 >  2:1754103/1‑100 (MQ=255)

                                gCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTc                                                                >  1:1887439/1‑100 (MQ=255)

                                gCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTc                                                                >  1:10185/1‑100 (MQ=255)

                                gCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTc                                                                >  2:951217/1‑100 (MQ=255)

                                 cTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                               >  1:1871429/1‑100 (MQ=255)

                                 cTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                               >  1:2758040/1‑100 (MQ=255)

                                  tCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctg                                                              >  2:1458404/1‑100 (MQ=255)

                                  tCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctg                                                              >  2:59335/1‑100 (MQ=255)

                                  tCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCACGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctg                                                              >  1:2061828/1‑100 (MQ=255)

                                   cGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                             >  2:3143406/1‑100 (MQ=255)

                                   cGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                             <  1:1173824/100‑1 (MQ=255)

                                   cGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                             <  1:811154/100‑1 (MQ=255)

                                   cGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                             >  2:1501810/1‑100 (MQ=255)

                                       aTCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTc                                                         <  2:382430/100‑1 (MQ=255)

                                       aTCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTc                                                         <  1:339136/100‑1 (MQ=255)

                                       aTCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTc                                                         <  1:2727288/100‑1 (MQ=255)

                                       aTCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTc                                                         <  2:2758084/100‑1 (MQ=255)

                                         cgcgGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGc                                                       >  2:921947/1‑100 (MQ=255)

                                          gcgGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCa                                                      <  2:867429/100‑1 (MQ=255)

                                          gcgGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCa                                                      >  2:2797631/1‑100 (MQ=255)

                                           cgGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTCATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAc                                                     <  1:2224412/100‑1 (MQ=255)

                                            ggTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAc                                                     >  1:1945628/1‑99 (MQ=255)

                                             gTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAcc                                                    >  1:1055818/1‑99 (MQ=255)

                                             gTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCAcc                                                    <  2:74728/99‑1 (MQ=255)

                                             gTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCg                                                   <  2:2786297/100‑1 (MQ=255)

                                                 ttACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGTTGTCTCACCGCCAg                                               >  1:592855/1‑100 (MQ=255)

                                                  tACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATTACCACCTCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGt                                              >  1:1534781/1‑100 (MQ=255)

                                                     tCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTaaa                                           <  1:2453109/100‑1 (MQ=255)

                                                      cATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAc                                          >  1:747659/1‑100 (MQ=255)

                                                      cATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAc                                          >  2:1105995/1‑100 (MQ=255)

                                                       aTCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAcc                                         >  1:3236496/1‑100 (MQ=255)

                                                       aTCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAcc                                         >  2:1414746/1‑100 (MQ=255)

                                                        tCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCa                                        >  2:2080593/1‑100 (MQ=255)

                                                        tCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCa                                        <  1:1754078/100‑1 (MQ=255)

                                                          ggggTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAgg                                      >  1:2826947/1‑100 (MQ=255)

                                                          ggggTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAgg                                      >  1:609505/1‑100 (MQ=255)

                                                            ggTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcg                                    <  2:3137047/100‑1 (MQ=255)

                                                             gTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcgc                                   >  2:1340463/1‑100 (MQ=255)

                                                                gCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTTGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGa                                 >  1:26945/1‑99 (MQ=255)

                                                                gCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                                >  2:1294559/1‑100 (MQ=255)

                                                                 cAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc                               >  1:787724/1‑100 (MQ=255)

                                                                 cAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc                               >  2:2551597/1‑100 (MQ=255)

                                                                  aaTATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAcc                              <  2:2971855/100‑1 (MQ=255)

                                                                  aaTATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAcc                              >  1:2693384/1‑100 (MQ=255)

                                                                   atatGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCt                             >  1:2202806/1‑100 (MQ=255)

                                                                     atGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGt                           <  2:257436/100‑1 (MQ=255)

                                                                      tGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGtt                          >  1:146120/1‑100 (MQ=255)

                                                                      tGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGtt                          >  1:1953721/1‑100 (MQ=255)

                                                                          cGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAt                      >  1:2389413/1‑100 (MQ=255)

                                                                           gTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtt                     <  1:1125254/100‑1 (MQ=255)

                                                                           gTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAtt                     <  1:3375996/100‑1 (MQ=255)

                                                                              tATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAttttt                  <  2:2954372/100‑1 (MQ=255)

                                                                               aTCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAttgtt                  >  1:3136595/1‑99 (MQ=255)

                                                                               aTCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg                 <  2:1955109/100‑1 (MQ=255)

                                                                               aTCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg                 <  2:1953750/100‑1 (MQ=255)

                                                                               aTCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg                 <  2:3131067/100‑1 (MQ=255)

                                                                                tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg                <  1:1775157/100‑1 (MQ=255)

                                                                                tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg                >  1:3297212/1‑100 (MQ=255)

                                                                                 ccAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGt               <  2:2194700/100‑1 (MQ=255)

                                                                                   aGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAg             >  2:135468/1‑100 (MQ=255)

                                                                                    gTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAgg            >  1:2546132/1‑100 (MQ=255)

                                                                                     tGATTCGCCCCCATGTCGTGATAATCACCGCCGCCAGCGGCTTTTTCTGCTGTCGCTCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAg             <  1:2391281/98‑1 (MQ=255)

                                                                                       aTTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc         >  2:1594495/1‑100 (MQ=255)

                                                                                        ttCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCa        >  1:2454894/1‑100 (MQ=255)

                                                                                          cGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatn      <  1:2057271/100‑2 (MQ=255)

                                                                                          cGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatg      >  1:1044196/1‑100 (MQ=255)

                                                                                          cGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatg      >  2:1273468/1‑100 (MQ=255)

                                                                                           gCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga     >  2:1100905/1‑100 (MQ=255)

                                                                                            cccccAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat    <  1:2614106/100‑1 (MQ=255)

                                                                                            cccccAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat    >  1:1291073/1‑100 (MQ=255)

                                                                                            cccccAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat    >  1:1289316/1‑100 (MQ=255)

                                                                                             ccccAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgatg   >  2:2316946/1‑100 (MQ=255)

                                                                                              cccAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGt  <  1:1146161/100‑1 (MQ=255)

AAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGT  >  NC_000913/2361226‑2361419

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