BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,265,670	T→G	17.6%	F453V (TTT→GTT)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,265,670	0	T	G	17.6%	299.0 / 65.9	137	F453V (TTT→GTT)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand): ref base T (60/52); new base G (12/12); total (73/64)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

Only 100 of 142 total aligned reads displayed.

GTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTAT  >  NC_000913/4265572‑4265769

gTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCGATCGGGACGGTACGCCTGACCtt                                                                                                    >  1:506881/1‑100 (MQ=255)

gTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCtt                                                                                                    >  2:1307092/1‑100 (MQ=255)

 tgatgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCttt                                                                                                   >  2:23301/1‑100 (MQ=255)

 tgatgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACATAACGGCCCAATCGGGACGGTACGCCTGACCttt                                                                                                   <  1:2431959/100‑1 (MQ=255)

   atgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTaa                                                                                                 >  2:2413842/1‑100 (MQ=255)

   atgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTaa                                                                                                 <  1:1954721/100‑1 (MQ=255)

   atgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGAAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTaa                                                                                                 >  1:1417537/1‑100 (MQ=255)

    tgaCGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAAc                                                                                                <  1:1596761/100‑1 (MQ=255)

        gtgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTc                                                                                            <  2:2216102/100‑1 (MQ=255)

         tgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGAACGGTACGCCTGACCTTTAACGGTCa                                                                                           >  1:991998/1‑100 (MQ=255)

           tATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCagg                                                                                         <  2:2226860/100‑3 (MQ=255)

            aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATTGGGACGGTACGCCTGACCGTTAAAGGTCaggc                                                                                        <  2:2458452/100‑4 (MQ=255)

            aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggc                                                                                        <  1:487925/100‑4 (MQ=255)

            aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggc                                                                                        <  2:1429741/100‑4 (MQ=255)

            aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggc                                                                                        <  2:1484162/100‑4 (MQ=255)

             tCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATgg                                                                                       >  2:1643048/1‑100 (MQ=255)

              cACGAAAACAGTGATTTAAAAGGCATCGCGGGAATTATTTTCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTCTAACGGTCAATGGt                                                                                      <  1:993929/100‑1 (MQ=255)

              cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGt                                                                                      >  2:326809/1‑100 (MQ=255)

              cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGt                                                                                      >  2:2084302/1‑100 (MQ=255)

                cGAGAACAGTGATGTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCTATCGGGACGGTACGCCTGTCCTTTAACGGTCAATGGTcg                                                                                    <  1:226249/100‑1 (MQ=255)

                cGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcg                                                                                    <  2:773892/100‑1 (MQ=255)

                cGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcg                                                                                    >  2:616101/1‑100 (MQ=255)

                 gAAAACAGTGATTTAAAGGGCAACAGGGAAATAATTATCGGTAGACAACGTAACGGGCCAATAGGGACGGTACGCCTGACCGTTAAAGGTCaggcgtacc                                                                                   <  1:618246/100‑9 (MQ=255)

                  aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg                                                                                  >  2:1149914/1‑100 (MQ=255)

                  aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg                                                                                  >  2:1149854/1‑100 (MQ=255)

                   aaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggcgtaccgt                                                                                 <  2:1277175/100‑11 (MQ=38)

                      cAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGGAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTc                                                                              >  1:2264472/1‑100 (MQ=255)

                        gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGa                                                                            >  1:777664/1‑100 (MQ=255)

                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAAGAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                          >  2:1824357/1‑100 (MQ=255)

                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                          >  2:581719/1‑100 (MQ=255)

                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                          >  1:870065/1‑100 (MQ=255)

                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                          >  1:1522319/1‑100 (MQ=255)

                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                          <  2:446791/100‑1 (MQ=255)

                          gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                          <  1:2245416/100‑1 (MQ=255)

                           aTTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACaa                                                                         >  1:855957/1‑100 (MQ=255)

                            tttAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAAc                                                                        >  2:185085/1‑100 (MQ=255)

                            tttAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTCACCTTTAACGGCCAATGGTCGCGCCTcaacaa                                                                         >  1:2153136/1‑99 (MQ=255)

                            tttAAAAGACATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAAc                                                                        <  2:2343788/100‑1 (MQ=255)

                             ttAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACt                                                                       >  1:716062/1‑100 (MQ=255)

                             ttAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACt                                                                       >  2:735975/1‑100 (MQ=255)

                                 aaGGGATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGTGACGGTAGGGCTATCCTTTAACGGACAATGGTCGCGGTTCGACAACTa                                                                      >  1:644353/1‑97 (MQ=255)

                                 aaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGc                                                                   >  2:1223829/1‑100 (MQ=255)

                                    gCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCggg                                                                <  2:2444662/100‑1 (MQ=255)

                                     cATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgggg                                                               <  1:951407/100‑1 (MQ=255)

                                       tCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGcc                                                             >  2:1008463/1‑100 (MQ=255)

                                        cgcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGCCAACTATGCGGGGCCg                                                            >  2:1622988/1‑100 (MQ=255)

                                        cgcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCg                                                            >  2:888553/1‑100 (MQ=255)

                                            gAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGt                                                        >  1:2006815/1‑100 (MQ=255)

                                                 cttccgatctAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgac                                                   <  2:1698684/91‑1 (MQ=255)

                                                     atCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacgacg                                               >  2:2417738/1‑100 (MQ=255)

                                                        ggTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGCCCTTTAACGGTCAATGGTCGCGCGTCGACAACGATGCGGGGCCGCAGTATGACGACGaat                                            <  1:232756/100‑1 (MQ=255)

                                                        ggTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaat                                            >  2:1140601/1‑100 (MQ=255)

                                                          tAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataa                                          <  2:1027200/100‑1 (MQ=255)

                                                           aaaCAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacagatcggaag                                         >  1:1698675/1‑90 (MQ=255)

                                                            aacaacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                        <  1:1141363/100‑1 (MQ=255)

                                                                  gTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAattatt                                  <  2:548561/100‑1 (MQ=255)

                                                                  gTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAattatt                                  >  1:325351/1‑100 (MQ=255)

                                                                   taccGGCCCAATCGGGACGGTACCCCTGTCCTTTAACGGTCAATGGTCGCGCTTAGAAAACTATGCGGGGCCGCAGTACGACGACGAATAATCattatt                                  <  1:248427/96‑1 (MQ=255)

                                                                     aCGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTa                               <  2:1455878/100‑1 (MQ=255)

                                                                       ggCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATg                             <  2:2483302/100‑1 (MQ=255)

                                                                        gCCCAATCGGGACGGTACGCCTGCCCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGa                            >  1:671420/1‑100 (MQ=255)

                                                                        gCCCAATCGGGACGGTACGCCTGACAATTAACGGTAAATGGTCGCGCTTAGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATg                             >  1:212054/1‑99 (MQ=255)

                                                                         cccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGaa                           >  2:1819006/1‑100 (MQ=255)

                                                                         cccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGaa                           >  1:1553771/1‑100 (MQ=255)

                                                                         cacAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTCCGACGCCGAATAATAATTATTTTATGaa                           >  1:36379/3‑100 (MQ=255)

                                                                          ccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                          >  1:1399015/1‑100 (MQ=255)

                                                                          ccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCACTCCGACGACGAATACTAATTATTTTATGaa                           >  1:1959453/1‑99 (MQ=255)

                                                                           cgattgggccGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAtt                         >  2:1708988/10‑100 (MQ=255)

                                                                           cAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAtt                         <  2:1788978/100‑1 (MQ=255)

                                                                              tCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAgg                      <  2:1921686/100‑1 (MQ=255)

                                                                               cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                     >  2:2268317/1‑100 (MQ=255)

                                                                               cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                     <  1:2413827/100‑1 (MQ=255)

                                                                               cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                     <  2:1311979/100‑1 (MQ=255)

                                                                               cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                     <  2:2205209/100‑1 (MQ=255)

                                                                                gggACGGTACGCCTGACATTTAACGGTCAATTGTCGCGATTCGACAACTATGCGTGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAgg                      <  1:811165/98‑1 (MQ=255)

                                                                                     ccgatctcTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTaa               <  1:1970442/93‑1 (MQ=255)

                                                                                      cgaccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTaaa              >  1:1913111/8‑100 (MQ=255)

                                                                                       gaccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAg             >  1:189778/7‑100 (MQ=255)

                                                                                       gaccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAg             >  2:2480193/7‑100 (MQ=255)

                                                                                       gaccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAg             <  1:1145291/94‑1 (MQ=255)

                                                                                        accattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc            >  2:907603/6‑100 (MQ=255)

                                                                                        aCGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc            <  2:1792151/100‑1 (MQ=255)

                                                                                        aCGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc            <  2:1790550/100‑1 (MQ=255)

                                                                                        aCGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc            <  2:540636/100‑1 (MQ=255)

                                                                                         ccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCa           >  2:2537855/5‑100 (MQ=255)

                                                                                         cGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCa           <  2:1900210/100‑1 (MQ=255)

                                                                                         cGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCa           <  2:1831553/100‑1 (MQ=255)

                                                                                          gCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaa          >  2:1972753/1‑100 (MQ=255)

                                                                                           ccTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACAACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaaa         >  1:241062/1‑100 (MQ=255)

                                                                                           attGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGCATAATCATTATTTTCTGAATTAGGTAATTAAAGCaaa         >  1:1744126/3‑100 (MQ=255)

                                                                                           attGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaaa         >  2:861150/3‑100 (MQ=255)

                                                                                           attGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaaa         >  2:1760654/3‑100 (MQ=255)

                                                                                            cTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAgatcgg        >  2:1970453/1‑95 (MQ=255)

                                                                                             tGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAaca       <  1:423197/100‑1 (MQ=255)

                                                                                              gACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAacac      >  1:851227/1‑100 (MQ=255)

                                                                                               aCCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACt     <  2:2379356/100‑1 (MQ=255)

                                                                                                ccTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACtt    >  1:841171/1‑100 (MQ=255)

                                                                                                 cTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTa   >  2:1220137/1‑100 (MQ=255)

                                                                                                 cTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTa   >  1:569213/1‑100 (MQ=255)

                                                                                                  tttAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTAt  >  1:1310676/1‑100 (MQ=255)

GTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTAT  >  NC_000913/4265572‑4265769

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