BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,222	A→C	I156L (ATC→CTC)	rhtC →	threonine efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,222	0	A	C	98.6%	479.5 / NA	141	I156L (ATC→CTC)	rhtC	threonine efflux pump
Reads supporting (aligned to +/- strand): ref base A (0/0); major base C (71/68); minor base T (2/0); total (73/68)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

CGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGATCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATT  >  NC_000913/4008125‑4008317

cGCTAATCCGAAAGCGATTATCTACTTTGGCTCTGTGTTCTCATTGTTTGTCGCTGATAACGTTGGCACTACCGCGTGCTGGGGCATTTTTGCGCTGCTc                                                                                               <  2:628976/100‑1 (MQ=255)

cGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTc                                                                                               <  1:2135310/100‑1 (MQ=255)

 gCTCATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCCGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCa                                                                                              >  2:786183/1‑100 (MQ=255)

   tAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCAtt                                                                                            <  1:550027/100‑1 (MQ=255)

    aaTCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTTATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTg                                                                                           >  2:3626562/1‑100 (MQ=255)

    aaTCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTg                                                                                           <  1:3914001/100‑1 (MQ=255)

    aaTCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTg                                                                                           <  2:62848/100‑1 (MQ=255)

    aaTCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTg                                                                                           <  1:3717508/100‑1 (MQ=255)

     aTCCGAAAGCGATTATCTACTTTGGCTCCGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGt                                                                                          <  1:3705442/100‑1 (MQ=255)

      tCCTAAAGCGATTATCGTCTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATCTTTGCGCTGCTTATTGTc                                                                                         <  1:3407277/100‑1 (MQ=255)

      tCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTc                                                                                         <  2:829689/100‑1 (MQ=255)

      tCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTc                                                                                         >  1:2403701/1‑100 (MQ=255)

      tCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTc                                                                                         <  1:1785776/100‑1 (MQ=255)

       ccGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCg                                                                                        <  1:2110620/100‑1 (MQ=255)

         gAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGaa                                                                                      >  1:2063932/1‑100 (MQ=255)

         gAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGaa                                                                                      <  2:221609/100‑1 (MQ=255)

         gAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGaa                                                                                      >  2:3300985/1‑100 (MQ=255)

         gAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTCCGCTGCTCCTTGTCGaa                                                                                      >  2:3959866/1‑100 (MQ=255)

            aGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACg                                                                                   <  2:1130573/100‑1 (MQ=255)

               gATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTg                                                                                <  1:415320/100‑1 (MQ=255)

                aTTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGGGCAGCTCATTGTCGAAACGCTgg                                                                               >  2:3406241/1‑100 (MQ=255)

                aTTATCTACTTTGGCTCGGTGTTCTCATTATTTGTCGGTGATAACGTTGGCACTACCGCGGGCTGGGGCATTTTTGCGGTGCTCATTGTCGAAACGCTgg                                                                               >  1:3804368/1‑100 (MQ=255)

                    tCTACTTGGGCTCGGTGGGGTCAGTGTTTGTCGGTTATAACGTTGTCACTACCGCGCACTGGGGCATGTTTGCGCTGCTCATTGTCGAAACGCTTGCGTg                                                                           <  2:3020282/100‑1 (MQ=255)

                      tACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGt                                                                         <  2:2485627/100‑1 (MQ=255)

                      tACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGt                                                                         >  2:2892384/1‑100 (MQ=255)

                      tACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGt                                                                         >  1:2742889/1‑100 (MQ=255)

                      tACTTTGGCTCGGTGTTCTCATTGTCTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGt                                                                         >  2:3187938/1‑100 (MQ=255)

                        cTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGttt                                                                       >  2:1938456/1‑100 (MQ=255)

                          gtGGCTCGGTGTGCTCATTGTTTGTCGGTGATAACGTGGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTAGCGTGGTTTAc                                                                     <  2:1510117/99‑1 (MQ=255)

                            ggCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACcg                                                                   <  2:3465402/100‑1 (MQ=255)

                             gCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACcgt                                                                  >  2:1075525/1‑100 (MQ=255)

                             gCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACcgt                                                                  >  1:3190191/1‑100 (MQ=255)

                               tCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACcgtcg                                                                <  1:4052352/100‑1 (MQ=255)

                                    gtTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGcc                                                           <  2:1817987/100‑1 (MQ=255)

                                     ttCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCa                                                          >  1:2616664/1‑100 (MQ=255)

                                     ttCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCa                                                          <  2:2390180/100‑1 (MQ=255)

                                      tctcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAg                                                         >  1:3346395/1‑100 (MQ=255)

                                      tctcATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAg                                                         >  1:19346/1‑100 (MQ=255)

                                         cATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCt                                                      >  2:1773841/1‑100 (MQ=255)

                                         cATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCt                                                      >  1:184714/1‑100 (MQ=255)

                                         cATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCt                                                      >  2:2616230/1‑100 (MQ=255)

                                         cATTGTTTGTCGGTGATAACGTTGCCACTACCGCGCGCTGTGGCATTTTTGCGCGGCGCAGTGTCGAAACGCTCGCGTGGTTTACCGTCGTTGCCAGCCt                                                      >  1:3811422/1‑100 (MQ=255)

                                          aTTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTg                                                     >  1:2508370/1‑100 (MQ=255)

                                           ttgtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTTTTTTTTTTTTTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGGCTGt                                                    >  2:1655176/1‑100 (MQ=255)

                                           ttgtttatCGGTGATAACGTTGGCACTACAGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGt                                                    >  2:1654291/1‑100 (MQ=255)

                                            tgtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGtt                                                   <  2:2853620/100‑1 (MQ=255)

                                             gtttgtCGGTGATAACGTTGGCACTACTGCGCGCTGGGGCATTTTTGCGCTGCTCATTGCCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGttt                                                  <  1:3534991/100‑1 (MQ=255)

                                             gtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGttt                                                  >  1:741720/1‑100 (MQ=255)

                                             gtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGttt                                                  >  1:522259/1‑100 (MQ=255)

                                             gtttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGttt                                                  >  2:2529242/1‑100 (MQ=255)

                                              tttgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTg                                                 <  1:2348218/100‑1 (MQ=255)

                                                tgtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGcc                                               <  1:3898448/100‑1 (MQ=255)

                                                 gtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGccc                                              >  2:2231755/1‑100 (MQ=255)

                                                 gtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGccc                                              >  1:3577477/1‑100 (MQ=255)

                                                 gtCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGccc                                              <  1:1325125/100‑1 (MQ=255)

                                                  tCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCt                                             <  2:2251311/100‑1 (MQ=255)

                                                   cGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTg                                            <  1:366381/100‑1 (MQ=255)

                                                   cGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTg                                            >  1:4062292/1‑100 (MQ=255)

                                                    ggTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGc                                           <  2:1799142/100‑1 (MQ=255)

                                                     gTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGcc                                          <  2:1063964/100‑1 (MQ=255)

                                                           aCGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCTTTGTCGAAACGCTGGCGTGCTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAAt                                    >  2:2524598/1‑100 (MQ=255)

                                                            cGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATg                                   <  2:3326957/100‑1 (MQ=255)

                                                              ttGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATgcg                                 >  2:2420772/1‑100 (MQ=255)

                                                              ctGGCTCTACCGCTCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATgcg                                 <  2:3724000/99‑1 (MQ=255)

                                                               tGGCACTACCGCGCGCTGGGGCATTTTTTCGCTTCTCATTGTCTAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATgcgc                                >  2:1521230/1‑100 (MQ=255)

                                                                ggCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGcc                               <  2:3797095/100‑1 (MQ=255)

                                                                ggCACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGcc                               >  1:2112719/1‑100 (MQ=255)

                                                                  cACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGt                             >  2:1306085/1‑100 (MQ=255)

                                                                  cACTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGt                             >  1:1837096/1‑100 (MQ=255)

                                                                   aCTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTg                            >  1:3303530/1‑100 (MQ=255)

                                                                    cTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTgg                           <  1:1097603/100‑1 (MQ=255)

                                                                    cTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTgg                           >  2:429378/1‑100 (MQ=255)

                                                                    cTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTgg                           >  1:4001370/1‑100 (MQ=255)

                                                                    cTACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGACGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGCgg                           >  2:2499176/1‑100 (MQ=255)

                                                                     tACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGt                          <  2:1416223/100‑1 (MQ=255)

                                                                     tACCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGt                          >  2:2565984/1‑100 (MQ=255)

                                                                      aCCGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGtt                         <  1:3423255/100‑1 (MQ=255)

                                                                       ccGCGCGCTGGGGCATTTTTTTTTTTTTTTTTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTa                        >  2:2291729/1‑100 (MQ=255)

                                                                       ccGCGCGCTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTa                        >  1:1233102/1‑100 (MQ=255)

                                                                           gcgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCaa                    >  2:2752820/1‑100 (MQ=255)

                                                                           gcgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCaa                    <  1:437633/100‑1 (MQ=255)

                                                                           gcgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCaa                    <  1:1051121/100‑1 (MQ=255)

                                                                            cgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAAc                   <  1:1592201/100‑1 (MQ=255)

                                                                            cgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCATGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAAc                   <  2:3888159/100‑1 (MQ=255)

                                                                            cgcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGAGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTTCCGCAACTGCGCCGTGGTTATCAAc                   >  2:3316784/1‑100 (MQ=255)

                                                                             gcTGGGGCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACg                  >  2:535727/1‑100 (MQ=255)

                                                                                ggggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCt               <  1:2617485/100‑1 (MQ=255)

                                                                                 gggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTg              >  2:2290249/1‑100 (MQ=255)

                                                                                 gggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTg              >  1:3171937/1‑100 (MQ=255)

                                                                                  ggCATTTTTGCGCTGCTCATTGTGGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTgg             >  1:2990012/1‑100 (MQ=255)

                                                                                  ggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCtcg             >  2:1907478/1‑98 (MQ=255)

                                                                                  ggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTgg             >  1:328921/1‑100 (MQ=255)

                                                                                  ggCATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTgg             >  1:912419/1‑100 (MQ=255)

                                                                                    cATTTTTGCGCTGCTCGTTGTCGAAACGCTGGCGTGGTTTACCTTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCg           <  1:218088/100‑1 (MQ=255)

                                                                                    cATTTTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCg           >  1:2963096/1‑100 (MQ=255)

                                                                                     aTTGTTGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGa          <  1:2711107/100‑1 (MQ=255)

                                                                                       ttttGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAg        <  2:3556360/100‑1 (MQ=255)

                                                                                        tttGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGt       >  2:3561150/1‑100 (MQ=255)

                                                                                          tGCGCTGCTCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTgg     >  1:1114995/1‑100 (MQ=255)

                                                                                             gctgctCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGAtt  <  1:1639463/100‑1 (MQ=255)

CGCTAATCCGAAAGCGATTATCTACTTTGGCTCGGTGTTCTCATTGTTTGTCGGTGATAACGTTGGCACTACCGCGCGCTGGGGCATTTTTGCGCTGATCATTGTCGAAACGCTGGCGTGGTTTACCGTCGTTGCCAGCCTGTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATT  >  NC_000913/4008125‑4008317

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