BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,246,073	C→T	100%	G435D (GGC→GAC)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,246,073	0	C	T	98.4%	385.1 / ‑4.5	114	G435D (GGC→GAC)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (59/53); total (60/54)
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.95e-01

Only 100 of 116 total aligned reads displayed.

GTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGCCG  >  NC_000913/1245976‑1246169

gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCGTCTCATTGTCGCCCACTTTTCACCGACGTGGTGTTCAGgtcg                                                                                                >  2:2281266/1‑100 (MQ=255)

   gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                             >  1:326418/1‑100 (MQ=255)

     tttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCggg                                                                                           <  1:1141217/100‑1 (MQ=255)

     tttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCggg                                                                                           >  2:1361293/1‑100 (MQ=255)

     tttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCggg                                                                                           >  2:1361892/1‑100 (MQ=255)

       tGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGtt                                                                                         >  2:3053498/1‑100 (MQ=255)

        gaaaTCCTTCTGTCGCGACCCACTGTAACGGAGCCCAGCCATGTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTg                                                                                        >  1:2194523/3‑100 (MQ=255)

        gTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCGCTTTTCACCGCCGTGGTGTTCAGGTCGCCGGGTTg                                                                                        >  2:1160054/1‑100 (MQ=255)

        gTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTg                                                                                        >  2:2372381/1‑100 (MQ=255)

        gTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTg                                                                                        <  1:2775060/100‑1 (MQ=255)

         tAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGtt                                                                                         >  1:2625772/1‑98 (MQ=255)

          aaTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgca                                                                                      >  1:688874/1‑100 (MQ=255)

          aaTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgca                                                                                      >  1:2225602/1‑100 (MQ=255)

           aTCCTTCTGTCGCGACCCACTGTATCGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgcag                                                                                     <  2:610548/100‑1 (MQ=255)

           aTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgca                                                                                      <  2:525419/99‑1 (MQ=255)

            tCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagc                                                                                    >  1:216339/1‑100 (MQ=255)

             ccTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagca                                                                                   >  2:220669/1‑100 (MQ=255)

              cTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGcagcag                                                                                  <  1:986084/100‑1 (MQ=255)

              cTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGcagcag                                                                                  >  1:316891/1‑100 (MQ=255)

                tCTGTCTCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAt                                                                                <  2:1157940/100‑1 (MQ=255)

                tCTGTCGCGACCCACTGTAACGGTGCCCATCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAt                                                                                <  1:1879594/100‑1 (MQ=255)

                  tGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgt                                                                              >  1:2847089/1‑100 (MQ=255)

                   gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgtg                                                                             >  1:194692/1‑100 (MQ=255)

                   gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgtg                                                                             >  1:2789902/1‑100 (MQ=255)

                   gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgtg                                                                             <  1:1838549/100‑1 (MQ=255)

                     cgcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGtt                                                                           >  1:616609/1‑100 (MQ=255)

                     cgcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGtt                                                                           <  2:2682398/100‑1 (MQ=255)

                         aCCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcg                                                                       >  1:90077/1‑100 (MQ=255)

                         aCCCACTGTAACGGTGCCAAGCCATTTAGCGCAACCCATTGATGCCCACTTTCCACCGACGTGGAGTTCAGGCCGCCGGGTTGCAGCAGAtgtgctatc                                                                        >  2:1822468/1‑94 (MQ=255)

                          cccACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgt                                                                      >  1:45594/1‑100 (MQ=255)

                           ccACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtc                                                                     >  2:1392094/1‑100 (MQ=255)

                            cgcTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcg                                                                    >  1:2796944/3‑100 (MQ=255)

                            cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcg                                                                    >  2:115500/1‑100 (MQ=255)

                            cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcg                                                                    >  2:2733512/1‑100 (MQ=255)

                            cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGACGTGTTTtcgtcg                                                                    <  1:2131063/100‑1 (MQ=255)

                             aCTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcac                                                                   >  2:1121069/1‑98 (MQ=255)

                               tGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccg                                                                 <  1:667412/100‑1 (MQ=255)

                                gTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgt                                                                >  1:2799541/1‑100 (MQ=255)

                                  aaCGGTGCCCCGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcg                                                              <  1:461555/100‑1 (MQ=255)

                                   aCGGTGCCCAGCCATTTGACGCATCCCATTGTTGCCCACCTTTCACCGACGTGCTGCTCAGGTCGCCGGGCTGCAGCAGATGTGTTTTcgtcgccgtcgc                                                             >  2:2038827/1‑100 (MQ=255)

                                     ggTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCa                                                           <  2:2590876/100‑1 (MQ=255)

                                     ggTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCa                                                           >  2:65796/1‑100 (MQ=255)

                                       tGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcgc                                                             <  2:2402434/96‑1 (MQ=255)

                                       tGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAtt                                                         >  1:2070141/1‑100 (MQ=255)

                                        gCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAttt                                                        >  1:2804134/1‑100 (MQ=255)

                                        gCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAttt                                                        >  1:452647/1‑100 (MQ=255)

                                         cccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAtttt                                                       <  1:459759/100‑1 (MQ=255)

                                         cccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAtttt                                                       >  2:3045414/1‑100 (MQ=255)

                                          ccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttg                                                      <  2:1276255/100‑1 (MQ=255)

                                          ccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttg                                                      <  2:2482830/100‑1 (MQ=255)

                                           cAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCATCGCCATTttgt                                                     <  2:2392447/100‑1 (MQ=255)

                                              ccATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTTCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgttg                                                   <  2:1173861/99‑1 (MQ=255)

                                              ccATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTgg                                                  <  2:1463611/100‑1 (MQ=255)

                                              ccATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTgg                                                  >  2:2786466/1‑100 (MQ=255)

                                              ccATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTgg                                                  <  1:2751642/100‑1 (MQ=255)

                                               cATTTGGCGCATCCCATTGTTACCCACTTTTCACCGACGTGGAGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGGCGCCATTTTGTTGgc                                                 >  2:1502371/1‑100 (MQ=255)

                                                agttGGCGCCGCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTgg                                                  <  2:2241772/96‑1 (MQ=255)

                                                     gcgcATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATc                                           <  1:2372397/100‑1 (MQ=255)

                                                      cgcATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCt                                          >  2:2188456/1‑100 (MQ=255)

                                                       gcATCCCATTGTTGCCCCCTTTTCCCCGACGTGGTGTTCAGCTCGCCGGGTTGCCGCAGACGTGTCTTCGTCCCCGTCCCCATTTTGTTCgc                                                 >  2:3230313/1‑92 (MQ=255)

                                                         aTCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTc                                       <  1:1593228/100‑1 (MQ=255)

                                                          tCCCATTGTTGCCTACTTTTCACCGACGGGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCTATTTTGTTTGCGCGATCTGTCg                                      <  2:750015/100‑1 (MQ=255)

                                                          tCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCg                                      >  1:336935/1‑100 (MQ=255)

                                                           cccATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGc                                     >  2:3001623/1‑100 (MQ=255)

                                                             cATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTg                                   >  2:484233/1‑100 (MQ=255)

                                                              aTTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGc                                  <  1:220671/100‑1 (MQ=255)

                                                               ttgttgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGcc                                 <  1:2507671/100‑1 (MQ=255)

                                                               ttgttgCCCACTTTTCACCGACGGGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCGGcc                                 <  2:1873629/100‑1 (MQ=255)

                                                                 gttgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGc                               >  1:2493019/1‑100 (MQ=255)

                                                                 gttgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGc                               <  1:2757362/100‑1 (MQ=255)

                                                                 attgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGc                               >  2:468790/2‑100 (MQ=255)

                                                                  ttgCCCGCTGCTCACCGACGTGCTGTGCAGGTCGCCGGGTGGCAGCAGATGTGTTTTCGTCCCCGTCGCCATTTTGTTGGCCCGATCTTTCGCTGCCg                                <  2:1877751/98‑1 (MQ=255)

                                                                    gCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCAtt                            >  1:2028005/1‑100 (MQ=255)

                                                                    gCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCAtt                            >  1:3146752/1‑100 (MQ=255)

                                                                    gCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCAtt                            >  2:2035938/1‑100 (MQ=255)

                                                                    gCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTACCGCAtt                            >  2:702899/1‑100 (MQ=255)

                                                                     cccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                           >  2:596611/1‑100 (MQ=255)

                                                                        aCTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACg                        <  1:1911626/100‑1 (MQ=255)

                                                                         gatctCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGt                       <  2:432675/96‑1 (MQ=255)

                                                                         cTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGt                       <  1:888741/100‑1 (MQ=255)

                                                                           tttCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTAc                     <  2:2225586/100‑1 (MQ=255)

                                                                              cACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTagatc                  >  1:432677/1‑96 (MQ=255)

                                                                              cACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACatc                  <  1:2188472/100‑3 (MQ=255)

                                                                              cACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGc                  <  2:140743/100‑1 (MQ=255)

                                                                              cACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGc                  <  2:183126/100‑1 (MQ=255)

                                                                               aCCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCg                 <  2:2236149/100‑1 (MQ=255)

                                                                                ccGTCGTGGTGGTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCg                 <  2:2789878/99‑1 (MQ=255)

                                                                                ccGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCgg                <  1:766709/100‑1 (MQ=255)

                                                                                 cGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg               >  1:3217064/1‑100 (MQ=255)

                                                                                 cGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg               <  1:480759/100‑1 (MQ=255)

                                                                                  gACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg               >  1:2853204/1‑99 (MQ=255)

                                                                                     gtggtgTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACa           >  2:2929633/1‑100 (MQ=255)

                                                                                      tggtgTTCAGGTCGACGGGTTGCAGCAGATGTGATTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAg          >  2:630500/1‑100 (MQ=255)

                                                                                        gtgtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAggg        <  1:218144/100‑1 (MQ=255)

                                                                                         tgtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggc       >  1:2204275/1‑100 (MQ=255)

                                                                                         tgtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggc       <  2:336933/100‑1 (MQ=255)

                                                                                          gtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcg      <  1:1788448/100‑1 (MQ=255)

                                                                                             cAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcggc    <  1:3053525/99‑1 (MQ=255)

                                                                                             cAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGcc   >  1:1273411/1‑100 (MQ=255)

                                                                                              aGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGCcg  >  1:1050049/1‑100 (MQ=255)

GTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGCCG  >  NC_000913/1245976‑1246169

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: `‑`