BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,182,806	C→A	15.2%	P86Q (CCA→CAA)	zupT →	zinc transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,806	0	C	A	15.2%	454.2 / 58.3	175	P86Q (CCA→CAA)	zupT	zinc transporter
Reads supporting (aligned to +/- strand): ref base C (90/57); new base A (13/14); total (103/72)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 194 total aligned reads displayed.

TTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCA  >  NC_000913/3182707‑3182905

ttCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc                                                                                                     <  2:1395231/100‑1 (MQ=255)

ttCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc                                                                                                     >  1:3510167/1‑100 (MQ=255)

 tCCTGCCTCACTAGCATCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGcc                                                                                                     <  2:3170768/99‑1 (MQ=255)

 tCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGTTTGGCTATTTTGGCCTGGACCGCATGTTGCca                                                                                                    <  1:3978142/100‑1 (MQ=255)

  ccTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcac                                                                                                   <  1:3655205/100‑1 (MQ=255)

   cTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcaca                                                                                                  <  2:1067880/100‑1 (MQ=255)

   cTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcaca                                                                                                  >  1:1291512/1‑100 (MQ=255)

    tGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACAt                                                                                                 >  1:3535792/1‑100 (MQ=255)

    tGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACAt                                                                                                 >  1:2239621/1‑100 (MQ=255)

     gCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATg                                                                                                >  1:2303457/1‑100 (MQ=255)

      ccGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGc                                                                                               >  1:3571701/1‑100 (MQ=255)

       cGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGcn                                                                                              <  1:652904/100‑2 (MQ=255)

        gCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTc                                                                                             >  1:1122424/1‑100 (MQ=255)

        gCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTc                                                                                             >  2:1966309/1‑100 (MQ=255)

         cACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                            <  2:68941/100‑4 (MQ=255)

         cACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                            <  1:414522/100‑4 (MQ=255)

             agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCg                                                                                        >  1:3447655/1‑100 (MQ=255)

             agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCg                                                                                        >  1:159753/1‑100 (MQ=255)

             agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggtccag                                                                                        <  2:3506387/100‑8 (MQ=255)

               cagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggtccaggc                                                                                      <  1:2367208/100‑10 (MQ=255)

               cagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCAGGACCGCATGTGGCAACATGcggt                                                                                            <  2:4018376/94‑4 (MQ=255)

                 gcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAgg                                                                                    >  2:2039343/1‑100 (MQ=255)

                 gcTGAAGGAATGGCGCCTGTGTTGGGTTATGGGATGTTTATCGTCGGTCGGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggtcc                                                                                          <  2:3590726/94‑6 (MQ=255)

                   tGAAGGAATGTTGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCGTGGCTATTTTGGCCTGGACCGCATGTTGCCAAAAGCTCATCAGAAGGa                                                                                   >  2:2612632/1‑99 (MQ=255)

                   tGAAGGAATGTCGCCTGTGTTGGGTTGTGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAt                                                                                  >  1:2544070/1‑100 (MQ=255)

                    gAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTCTTTTGGACTGGCCTCAATGTTGTCACATGCTCATTCGCAGGa                                                                                   >  2:1855929/1‑98 (MQ=255)

                     aaGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAttt                                                                                >  1:2758295/1‑100 (MQ=255)

                     aaGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAttt                                                                                >  2:82094/1‑100 (MQ=255)

                       ggAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTaa                                                                              >  1:3593902/1‑100 (MQ=255)

                       ggAATGTCGCCTGTGTTAGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTAGCATTGACCGAATGTTGCAACATGCTCATCCGCAGGATTTaa                                                                              >  2:3399607/1‑100 (MQ=255)

                        gAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAt                                                                             <  2:1710703/100‑1 (MQ=255)

                        gAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAt                                                                             <  2:1634783/100‑1 (MQ=255)

                         aaTGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATg                                                                            >  2:3854192/1‑100 (MQ=255)

                          aTGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGc                                                                           >  1:3572295/1‑100 (MQ=255)

                          aTGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGc                                                                           >  2:2998008/1‑100 (MQ=255)

                           tGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGCTTTAATGCa                                                                          >  1:3538702/1‑100 (MQ=255)

                              cGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCaaaa                                                                       >  2:2134001/1‑100 (MQ=255)

                                 cTGTGTTGGGTTATGGAATGTTTATCTTCGGGCTGCTTGGCTAGTTTGGCCTGGTCCGCATGGTGCCACATGCTCATCCGCAGGATTTAATGCaaaa                                                                       <  2:899060/97‑1 (MQ=255)

                                   gtgtTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCgg                                                                  >  1:136588/1‑100 (MQ=255)

                                   gtgtTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCgg                                                                  <  1:2586149/100‑1 (MQ=255)

                                    tgtTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGt                                                                 >  1:290795/1‑100 (MQ=255)

                                           ttATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCg                                                          >  1:1112929/1‑100 (MQ=255)

                                           ttATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGAAAAAATCGGTGCAGCCg                                                          <  2:3640269/100‑1 (MQ=255)

                                              tGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTg                                                       >  1:2259755/1‑100 (MQ=255)

                                               ggAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAGCACATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAACCGGTGCAGCCCTTGc                                                      >  2:4189809/1‑100 (MQ=255)

                                               ggAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGc                                                      >  1:389742/1‑100 (MQ=255)

                                                gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGCTTTAATGCAAAAATCGGTGCAGCCGTTGcc                                                     >  1:1956964/1‑100 (MQ=255)

                                                gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAATTAATGAAAAAATCGGTGCAGCCGTTGcc                                                     <  2:95525/100‑1 (MQ=255)

                                                gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCAGGAACGCATGTTGCCACATGCTCATCCGCAGAAAAAAAAAAAAAAATCGGTGCAGCCGTTGc                                                      <  2:755115/99‑1 (MQ=255)

                                                   tGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGcc                                                     <  2:67345/97‑1 (MQ=255)

                                                     tttATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAAt                                                >  2:1065658/1‑100 (MQ=255)

                                                     tttATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATAGGTGCAGCCGTTGCCaaaa                                                 >  2:3336856/1‑99 (MQ=255)

                                                      ttATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGAATGTTGCCACATGCTCATCCGCAGAATAAAAAAAAAAAATCGGTGCAGCCGTTGCCAAAAt                                                <  2:216490/99‑1 (MQ=255)

                                                          cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           <  2:2520404/100‑1 (MQ=255)

                                                          cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           <  1:1731324/100‑1 (MQ=255)

                                                          cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           >  2:433777/1‑100 (MQ=255)

                                                           ttCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCCCATGCTCATCCGCAGGATTTAATGCAACAATCGGTGCAGCCGTTGCCAAAATCGATCa                                          >  2:1233211/1‑100 (MQ=255)

                                                           ttCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCa                                          >  1:3289449/1‑100 (MQ=255)

                                                            tCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           <  2:4070516/98‑1 (MQ=255)

                                                            tCGGTCTGCTTGGCTATTATGGACTGGACCGCATGTTGCCACATGCTCATCCGCAAAATTTAAAACAAAAATCGGTGCAGCCGTTGCCAAAATCGATCaa                                         <  2:204708/100‑1 (MQ=255)

                                                              ggTCTGCTT‑GCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAgcg                                      <  1:1964091/100‑1 (MQ=255)

                                                                  tGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGcac                                   >  1:2893532/1‑100 (MQ=255)

                                                                    cTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAg                                 >  1:870228/1‑100 (MQ=255)

                                                                     ttGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGc                                >  1:2138039/1‑100 (MQ=255)

                                                                        gCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCAt                             <  2:51257/100‑1 (MQ=255)

                                                                        gCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCAt                             >  1:107814/1‑100 (MQ=255)

                                                                         gtATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCAtt                            >  2:831561/2‑100 (MQ=255)

                                                                          tATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTc                           >  1:492569/1‑100 (MQ=255)

                                                                          tATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTc                           >  1:2816092/1‑100 (MQ=255)

                                                                          tATTTTGGCCTGGACCACATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTc                           <  1:4189838/100‑1 (MQ=255)

                                                                           aTTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCt                          <  2:595947/100‑1 (MQ=255)

                                                                           aTTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCt                          >  2:3028044/1‑100 (MQ=255)

                                                                            ttttGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTg                         >  1:504006/1‑100 (MQ=255)

                                                                              ttGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCt                       <  1:3695329/100‑1 (MQ=255)

                                                                               tGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTc                      <  1:3149681/100‑1 (MQ=255)

                                                                                   cggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCActc                  >  2:2846777/8‑100 (MQ=255)

                                                                                    tGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTCGCCAAAATCGATCAAGCGCACAGCCATTCTGCTc                      <  2:561299/95‑1 (MQ=255)

                                                                                    ggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCActct                 >  1:2321701/7‑100 (MQ=255)

                                                                                     gatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACtctc                >  2:1619661/6‑100 (MQ=255)

                                                                                       tgtgCATGTCGCATCATACTCATCCGCAGGATTTAATTCAAAAATCGGTGCAGCCGTTGCCATAATCGAACAATCGCACAGCCCTTCTACTCAcgccgta              >  2:1954371/4‑94 (MQ=255)

                                                                                       tgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              >  1:2467322/4‑100 (MQ=255)

                                                                                       aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              >  2:1158474/1‑100 (MQ=255)

                                                                                       aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              <  1:4063593/100‑1 (MQ=255)

                                                                                       aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              <  1:1682289/100‑1 (MQ=255)

                                                                                       aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              >  1:1203640/1‑100 (MQ=255)

                                                                                       aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              >  2:1020803/1‑100 (MQ=255)

                                                                                       aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCg               >  2:3914677/1‑99 (MQ=255)

                                                                                         cGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAATTCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCa            <  2:2239616/100‑1 (MQ=255)

                                                                                          gCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCAt           <  1:1874852/100‑1 (MQ=255)

                                                                                           cATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATc          >  1:3788769/1‑100 (MQ=255)

                                                                                           cATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATc          <  2:3275738/100‑1 (MQ=255)

                                                                                             tGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAg        >  2:1332498/1‑100 (MQ=255)

                                                                                              gTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAg        >  1:2440204/1‑99 (MQ=255)

                                                                                               ttGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTc      >  1:3915939/1‑100 (MQ=255)

                                                                                                tGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCCCAGCCATTCTGCTCACTCTCGGCATCAGTCt     >  2:3851875/1‑100 (MQ=255)

                                                                                                 gCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTg    >  1:1909772/1‑100 (MQ=255)

                                                                                                 gCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTg    >  2:2218031/1‑100 (MQ=255)

                                                                                                  ccACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGc   <  2:2674634/100‑1 (MQ=255)

                                                                                                   cacaTGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCa  >  2:276875/1‑100 (MQ=255)

                                                                                                   cacaTGCTCATCCGCAGGATTTAATGCAAAAATCCGTGCAGCCGTTGCCAAAATCGATCAAGCGCAAAGCCATTCTGCTCACTCTCGGCATCAGTCTGCa  >  2:3757833/1‑100 (MQ=255)

TTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCA  >  NC_000913/3182707‑3182905

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