BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,247,916	G→A	A206V (GCG→GTG)	lysP ←	lysine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,247,916	0	G	A	99.6%	977.2 / NA	238	A206V (GCG→GTG)	lysP	lysine:H(+) symporter
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (138/99); minor base T (0/1); total (138/100)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

GAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCGCCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAGCACGC  >  NC_000913/2247818‑2248012

 aaTGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATcaccg                                                                                               >  2:3663823/1‑101 (MQ=255)

  aTGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATcaccgc                                                                                              <  1:3663427/101‑1 (MQ=255)

    gTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCaa                                                                                            >  2:3663827/1‑101 (MQ=255)

     tttttCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCaaa                                                                                           >  1:3663828/1‑101 (MQ=255)

       tttCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAac                                                                                         <  1:3663518/101‑1 (MQ=255)

          cGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAaccac                                                                                      <  1:3663177/101‑1 (MQ=255)

          cGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAaccac                                                                                      <  1:3663354/101‑1 (MQ=255)

          cGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAaccac                                                                                      >  2:3663833/1‑101 (MQ=255)

           gCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAaccacc                                                                                     >  1:3663835/1‑101 (MQ=255)

           gCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAaccacc                                                                                     >  2:3663836/1‑101 (MQ=255)

              ggATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGc                                                                                  <  2:3663767/101‑1 (MQ=255)

              ggATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGc                                                                                  >  1:3663840/1‑101 (MQ=255)

                aTCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCaa                                                                                <  2:3663469/101‑1 (MQ=255)

                aTCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCaa                                                                                <  1:3663363/101‑1 (MQ=255)

                aTCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCaa                                                                                >  2:3663842/1‑101 (MQ=255)

                 tCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCaaa                                                                               >  2:3663847/1‑101 (MQ=255)

                 tCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCaaa                                                                               >  1:3663846/1‑101 (MQ=255)

                   ttCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCTCCGCAAAACCACCAGCAAACg                                                                             <  2:3663199/101‑1 (MQ=255)

                    tCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAAc                                                                              >  1:3663854/1‑99 (MQ=255)

                    tCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACgg                                                                            <  2:3663444/101‑1 (MQ=255)

                    tCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACgg                                                                            >  2:3663853/1‑101 (MQ=255)

                     cGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGgc                                                                           <  1:3663303/101‑1 (MQ=255)

                      ggACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGgcg                                                                          >  2:3663858/1‑101 (MQ=255)

                       gACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGgcgc                                                                         <  2:3663666/101‑1 (MQ=255)

                        aCTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCt                                                                        >  1:3663859/1‑101 (MQ=255)

                         cTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCtt                                                                       <  2:3663269/101‑1 (MQ=255)

                         cTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCtt                                                                       <  2:3663489/101‑1 (MQ=255)

                           cGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCg                                                                     <  1:3663453/101‑1 (MQ=255)

                           cGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCg                                                                     <  2:3663457/101‑1 (MQ=255)

                            gCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGc                                                                    >  2:3663861/1‑101 (MQ=255)

                              cTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCg                                                                  >  2:3663865/1‑101 (MQ=255)

                               tGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAg                                                                                   <  2:3663308/83‑1 (MQ=255)

                               tGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGa                                                                 >  1:3663868/1‑101 (MQ=255)

                               tGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGa                                                                 >  2:3663869/1‑101 (MQ=255)

                                 cagcaATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGAtt                                                               <  1:3663338/101‑1 (MQ=255)

                                    caATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTc                                                            >  1:3663874/1‑101 (MQ=255)

                                    caATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTc                                                            >  2:3663873/1‑101 (MQ=255)

                                     aaTACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTcc                                                           >  1:3663878/1‑101 (MQ=255)

                                      aTACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCa                                                          <  2:3663376/101‑1 (MQ=255)

                                       tACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAg                                                         >  2:3663880/1‑101 (MQ=255)

                                        aCCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGt                                                        >  1:3663882/1‑101 (MQ=255)

                                           gATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGc                                                     >  1:3663886/1‑101 (MQ=255)

                                            aTCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCt                                                    >  2:3663889/1‑101 (MQ=255)

                                              cAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTcc                                                  >  2:3663893/1‑101 (MQ=255)

                                              cAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTcc                                                  >  1:3663891/1‑101 (MQ=255)

                                               aGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCa                                                 <  2:3663677/101‑1 (MQ=255)

                                               aGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCa                                                 <  1:3663344/101‑1 (MQ=255)

                                                gCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAg                                                <  1:3663456/101‑1 (MQ=255)

                                                gCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAg                                                >  1:3663897/1‑101 (MQ=255)

                                                 cTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGc                                               <  2:3663432/101‑1 (MQ=255)

                                                  tCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGcc                                              >  2:3663901/1‑101 (MQ=255)

                                                  tCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGcc                                              >  1:3663900/1‑101 (MQ=255)

                                                   cGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGccc                                             <  1:3663334/101‑1 (MQ=255)

                                                    ggTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGAtt                                                               <  2:3663517/82‑1 (MQ=255)

                                                      ttCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCa                                          >  1:3663904/1‑101 (MQ=255)

                                                       tCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAg                                         <  2:3663359/101‑1 (MQ=255)

                                                        cccTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAgg                                        >  2:3663907/1‑101 (MQ=255)

                                                        cccTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAg                                         >  1:3663908/1‑100 (MQ=255)

                                                          cTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCt                                      >  1:3663911/1‑101 (MQ=255)

                                                           tGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTg                                     >  2:3663917/1‑101 (MQ=255)

                                                           tGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTg                                     >  2:3663915/1‑101 (MQ=255)

                                                            ggAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTgc                                    <  2:3663602/101‑1 (MQ=255)

                                                               aaGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTgcgcg                                 >  2:3663921/1‑101 (MQ=255)

                                                               aaGAGAAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTgcgcg                                 >  1:3663922/1‑101 (MQ=255)

                                                                 gagaAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGcc                               <  1:3663549/101‑1 (MQ=255)

                                                                  agaAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGc                                       >  2:3663930/1‑92 (MQ=255)

                                                                  agaAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCt                              <  2:3663573/101‑1 (MQ=255)

                                                                  agaAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCt                              <  1:3663378/101‑1 (MQ=255)

                                                                  agaAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCt                              >  2:3663926/1‑101 (MQ=255)

                                                                  agaAGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCt                              >  1:3663927/1‑101 (MQ=255)

                                                                    aaGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCttt                            >  1:3663935/1‑101 (MQ=255)

                                                                    aaGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCttt                            >  2:3663934/1‑101 (MQ=255)

                                                                    aaGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCttt                            <  1:3663709/101‑1 (MQ=255)

                                                                     aGCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCt                                      >  2:3663938/1‑90 (MQ=255)

                                                                      gCCGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGa                          <  1:3663782/101‑1 (MQ=255)

                                                                        cGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAg                        <  1:3663576/101‑1 (MQ=255)

                                                                        cGACAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAg                        <  2:3663565/101‑1 (MQ=255)

                                                                          aCAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGAt                      <  1:3663539/101‑1 (MQ=255)

                                                                           cAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAgg                                        >  1:3663941/1‑82 (MQ=255)

                                                                           cAATCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATa                     <  2:3663492/101‑1 (MQ=255)

                                                                            aaTCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATAc                    >  2:3663942/1‑101 (MQ=255)

                                                                            aaTCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATAc                    <  2:3663788/101‑1 (MQ=255)

                                                                             aTCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATAcc                   <  1:3663606/101‑1 (MQ=255)

                                                                             aTCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATAcc                   >  2:3663944/1‑101 (MQ=255)

                                                                              tCATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCg                  <  2:3663540/101‑1 (MQ=255)

                                                                               cATAGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGa                 >  1:3663947/1‑101 (MQ=255)

                                                                                  aGCTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATaa              >  2:3663951/1‑101 (MQ=255)

                                                                                    cTACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAatc            >  1:3663953/1‑101 (MQ=255)

                                                                                     tACGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAatca           <  2:3663405/101‑1 (MQ=255)

                                                                                      aCGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAatcat          >  2:3663954/1‑101 (MQ=255)

                                                                                       cGCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAatcatc         <  1:3663661/101‑1 (MQ=255)

                                                                                        gCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAAtcatca        <  1:3663502/101‑1 (MQ=255)

                                                                                        gCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAAtcatca        >  1:3663957/1‑101 (MQ=255)

                                                                                        gCCGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAAtcatca        <  2:3663416/101‑1 (MQ=255)

                                                                                         ccGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAg       <  1:3663487/101‑1 (MQ=255)

                                                                                          cGATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAGc      >  2:3663960/1‑101 (MQ=255)

                                                                                           gATCATCACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAGCa     <  1:3663418/101‑1 (MQ=255)

                                                                                            atcatcACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAGCAc    >  1:3663963/1‑101 (MQ=255)

                                                                                             tcatcaCCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAGCACg   >  2:3663964/1‑101 (MQ=255)

                                                                                              catcACCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAGCACGc  >  2:3663968/1‑101 (MQ=255)

GAATGTTTTTCGCCGGATCTTCGGACTCGCCTGCAGCAATACCGATCAGCTCGGTTCCCTGGAAAGAGAAGCCGACAATCATAGCTACGCCGATCATCGCCGCAAAACCACCAGCAAACGGCGCTTCGCCGATTGTCCAGTTGCTCCAGCCCGCAGGCTGCGCGCCTTTGAAGATACCGATAATCATCAGCACGC  >  NC_000913/2247818‑2248012

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