BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	272,747	T→C	L221P (CTA→CCA)	ykfC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	272,747	0	T	C	78.5%	523.6 / 167.1	284	L221P (CTA→CCA)	ykfC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (30/31); new base C (101/122); total (131/153)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-04

Only 100 of 310 total aligned reads displayed.

GCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGA  >  NC_000913/272655‑272836

  tgtgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCcat                                                                                         <  1:4176491/93‑1 (MQ=255)

  tgtgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCcat                                                                                         <  1:7238694/93‑1 (MQ=255)

   gtgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGGTGTTATATCGCCGCcat                                                                                         >  1:2364599/1‑93 (MQ=255)

   gtgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCcatt                                                                                        <  1:1591735/93‑1 (MQ=255)

    tgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCtatta                                                                                       <  1:2825310/93‑1 (MQ=255)

    tgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCcatta                                                                                       >  1:1801455/1‑93 (MQ=255)

    tgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCcatta                                                                                       >  2:3380157/1‑93 (MQ=255)

    tgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCAGCcatta                                                                                       <  1:6910111/93‑1 (MQ=255)

        aaaCCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGa                                                                                   <  1:5334474/93‑1 (MQ=255)

         aaCCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGaa                                                                                  <  2:217383/93‑1 (MQ=255)

           ccATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAca                                                                                >  2:6762681/1‑93 (MQ=255)

           ccATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAca                                                                                >  1:3509581/1‑93 (MQ=255)

             aTCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAcatc                                                                              >  2:5173070/1‑93 (MQ=255)

             aTCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAcatc                                                                              <  2:6591615/93‑1 (MQ=255)

              tCAGGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAcatca                                                                             >  1:7502208/1‑93 (MQ=255)

              tCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAcatca                                                                             <  2:2258208/93‑1 (MQ=255)

               cAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAcatcat                                                                            <  1:5768035/93‑1 (MQ=255)

               cAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAAcatcat                                                                            >  2:2581813/1‑93 (MQ=255)

                aaGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATg                                                                           >  1:497128/1‑93 (MQ=255)

                aaGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATg                                                                           >  1:4389304/1‑93 (MQ=255)

                  ggCGGGACATATCGATGTCGGTCTCTATCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCt                                                                         >  2:6626625/1‑93 (MQ=255)

                    cGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGa                                                                       >  1:4834851/1‑93 (MQ=255)

                     gggACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGaa                                                                      >  2:1340495/1‑93 (MQ=255)

                     gggACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGaa                                                                      <  1:2993978/93‑1 (MQ=255)

                     gggACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGaa                                                                      <  2:4193484/93‑1 (MQ=255)

                      ggACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAAt                                                                     <  1:3801268/93‑1 (MQ=255)

                      ggACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAAt                                                                     <  2:327399/93‑1 (MQ=255)

                       gACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATg                                                                    >  2:1104671/1‑93 (MQ=255)

                       gACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATg                                                                    <  1:2891483/93‑1 (MQ=255)

                       gACATATCGATGTCGGCCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATg                                                                    <  1:7064500/93‑1 (MQ=255)

                        aCATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGa                                                                   >  2:7392764/1‑93 (MQ=255)

                           tatCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGtt                                                                >  1:2131852/1‑93 (MQ=255)

                           tatCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGtt                                                                <  2:3481531/93‑1 (MQ=255)

                            atCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTc                                                               <  2:485127/93‑1 (MQ=255)

                               gATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGAt                                                            <  1:5886872/93‑1 (MQ=255)

                                aTGTCGGTCTCTTTCGGGCGGCCAGGGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATc                                                           <  1:1040820/93‑1 (MQ=255)

                                  gTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCaa                                                         <  2:831057/93‑1 (MQ=255)

                                  gTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCaa                                                         <  2:6177574/93‑1 (MQ=255)

                                   tCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAAt                                                        <  2:4569076/93‑1 (MQ=255)

                                    cGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACCTCATGCTGAATGAGTTCGATCAATa                                                       >  1:5982614/1‑93 (MQ=255)

                                    cGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATa                                                       <  2:706669/93‑1 (MQ=255)

                                     ggTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATAc                                                      <  1:2672400/93‑1 (MQ=255)

                                      gTCTCTTTCGGGGGGCAAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATAcc                                                     <  2:7236707/93‑1 (MQ=255)

                                      gTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATAcc                                                     >  1:1668333/1‑93 (MQ=255)

                                         tctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGc                                                  <  2:2936883/93‑1 (MQ=255)

                                         tctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACATGc                                                  >  2:5797935/1‑93 (MQ=255)

                                          ctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCa                                                 <  2:6507194/93‑1 (MQ=255)

                                          ctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCa                                                 <  2:1882099/93‑1 (MQ=255)

                                           tttCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCAt                                                <  1:4291228/93‑1 (MQ=255)

                                           tttCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCAt                                                <  2:1456332/93‑1 (MQ=255)

                                           tttCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCAt                                                <  2:6263498/93‑1 (MQ=255)

                                            ttCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATg                                               <  1:869190/93‑1 (MQ=255)

                                             tCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGa                                              <  2:1027765/93‑1 (MQ=255)

                                             tCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGa                                              <  2:5980826/93‑1 (MQ=255)

                                             tCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGa                                              <  1:1144478/93‑1 (MQ=255)

                                             tCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGa                                              <  2:25200/93‑1 (MQ=255)

                                                ggcggcCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAgcg                                           >  2:5571197/1‑93 (MQ=255)

                                                 gcggcCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAgcgc                                          <  2:5453452/93‑1 (MQ=255)

                                                  cggcCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCt                                         <  1:5265907/93‑1 (MQ=255)

                                                   ggcCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTa                                        <  1:5469319/93‑1 (MQ=255)

                                                     ccAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTAcc                                      <  2:3529182/93‑1 (MQ=255)

                                                      cAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCt                                     >  2:2708726/1‑93 (MQ=255)

                                                       aGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTg                                    >  2:6381735/1‑93 (MQ=255)

                                                       aGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTg                                    >  2:6131955/1‑93 (MQ=255)

                                                          gAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGc                                 >  1:736786/1‑93 (MQ=255)

                                                           aaGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCg                                >  2:4710199/1‑93 (MQ=255)

                                                           aaGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCg                                <  2:7467874/93‑1 (MQ=255)

                                                            aGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCgg                               <  1:3653367/93‑1 (MQ=255)

                                                            aGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCgg                               >  1:3064304/1‑93 (MQ=255)

                                                             ggTGTGCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCggg                              >  2:5230513/1‑93 (MQ=255)

                                                                 tgCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAg                          <  1:662063/93‑1 (MQ=255)

                                                                  gCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGc                         >  1:2470905/1‑93 (MQ=255)

                                                                  gCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGc                         >  2:3123332/1‑93 (MQ=255)

                                                                  gCCACAGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGc                         <  2:4342550/93‑1 (MQ=255)

                                                                    cacaGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCa                       >  1:6473618/1‑93 (MQ=255)

                                                                    cacaGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCa                       >  1:5123374/1‑93 (MQ=255)

                                                                    cacaGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCa                       <  1:5582881/93‑1 (MQ=255)

                                                                      caGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGa                     <  2:936164/93‑1 (MQ=255)

                                                                      caGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGa                     <  1:6347336/93‑1 (MQ=255)

                                                                       aGGGC‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGaa                    <  2:608267/93‑1 (MQ=255)

                                                                           c‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGa                <  2:5304798/93‑1 (MQ=255)

                                                                           c‑GGTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGa                <  1:1927909/93‑1 (MQ=255)

                                                                             ggTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGAt               >  1:3904308/1‑93 (MQ=255)

                                                                             ggTGTTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGAt               <  1:4597744/93‑1 (MQ=255)

                                                                              gtgtTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATc              >  2:361521/1‑93 (MQ=255)

                                                                              gtgtTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATc              >  1:6194055/1‑93 (MQ=255)

                                                                               tgtTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCg             <  1:4045123/93‑1 (MQ=255)

                                                                               tgtTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCg             <  2:1682780/93‑1 (MQ=255)

                                                                                gtTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCCGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCgg            >  1:1989788/1‑93 (MQ=255)

                                                                                gtTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCgg            >  2:3985418/1‑93 (MQ=255)

                                                                                gtTATATCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCgg            <  1:181728/93‑1 (MQ=255)

                                                                                 ttATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggt           >  2:4994066/1‑93 (MQ=255)

                                                                                 ttATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggt           >  2:4441940/1‑93 (MQ=255)

                                                                                  tataTCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggtg          >  1:1332212/1‑93 (MQ=255)

                                                                                   atatCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggtgg         >  2:6154171/1‑93 (MQ=255)

                                                                                   atatCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggtgg         <  1:7156889/93‑1 (MQ=255)

                                                                                    tatCGCCGCCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggtggt        >  2:1603880/1‑93 (MQ=255)

                                                                                        gccgccATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTg    >  2:4069466/1‑93 (MQ=255)

                                                                                         ccgcCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTgg   <  2:268883/93‑1 (MQ=255)

                                                                                          cgcCATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGa  <  1:2934159/93‑1 (MQ=255)

GCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGC‑GGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGA  >  NC_000913/272655‑272836

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