BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	784,012	T→A	15.2%	Q271L (CAG→CTG)	zitB ←	zinc efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	784,012	0	T	A	15.2%	412.4 / 60.2	165	Q271L (CAG→CTG)	zitB	zinc efflux system
Reads supporting (aligned to +/- strand): ref base T (65/75); new base A (9/16); total (74/91)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

Only 100 of 170 total aligned reads displayed.

CCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACAT  >  NC_000913/783913‑784111

ccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGctct                                                                                                     >  2:2334439/1‑100 (MQ=255)

ccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCt                                                                                                     <  2:1749610/100‑1 (MQ=255)

   tcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGAt                                                                                                  <  1:4007138/100‑1 (MQ=255)

    cattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATc                                                                                                 <  2:4017367/100‑1 (MQ=255)

     attcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATcc                                                                                                >  2:4259121/1‑100 (MQ=255)

     attcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATcc                                                                                                <  2:3105051/100‑1 (MQ=255)

     attcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATcc                                                                                                <  2:1455515/100‑1 (MQ=255)

      ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCca                                                                                               <  1:3335230/100‑1 (MQ=255)

       tcaGGTGGCAGTCTGGCCCATGACAAGGTTGATCTTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGGAGTGTTGGATCTGATCcaa                                                                                              <  1:989801/100‑1 (MQ=255)

         aGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaaca                                                                                            <  1:2053476/100‑1 (MQ=255)

          ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaacaa                                                                                           >  1:2112011/1‑100 (MQ=255)

          ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaacaa                                                                                           >  2:2730272/1‑100 (MQ=255)

          ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaacaa                                                                                           >  1:3838072/1‑100 (MQ=255)

          ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacac                                                                                           <  2:867540/100‑2 (MQ=255)

          ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCaaacac                                                                                           <  1:1019454/100‑2 (MQ=255)

           gTGGCAGTCTGGCCCATGACAAGGTTTATATTCCATCTGAATGGTGGCGTGCTCAATCTGCTAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacact                                                                                          <  1:1601425/100‑3 (MQ=255)

           gTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacact                                                                                          <  2:2460383/100‑3 (MQ=255)

           gTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacact                                                                                          <  2:93361/100‑3 (MQ=255)

            tGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacacta                                                                                         <  2:3136343/100‑4 (MQ=255)

             ggCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATATTGATCCATCAGGTAGTGTTGGATCAGATCcaacactac                                                                                        <  1:2706314/100‑5 (MQ=255)

              gCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacactacc                                                                                       <  1:3418427/100‑6 (MQ=255)

              gCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacactacc                                                                                       <  2:4301069/100‑6 (MQ=255)

               cAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCat                                                                                      >  2:243462/1‑100 (MQ=255)

               cAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCat                                                                                      >  2:2196556/1‑100 (MQ=255)

               cAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCat                                                                                      >  1:165525/1‑100 (MQ=255)

               cAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacactacct                                                                                      <  2:1090055/100‑7 (MQ=255)

                 gTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGTATCTGATCCAACAAGGCatcg                                                                                    <  1:4329486/100‑1 (MQ=255)

                 gTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcg                                                                                    <  2:1346648/100‑1 (MQ=255)

                 gTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcg                                                                                    <  2:3397560/100‑1 (MQ=255)

                  tCTGGCGCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacactacctgat                                                                                   <  1:1296109/100‑10 (MQ=38)

                     ggCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgat                                                                                <  1:274267/100‑1 (MQ=255)

                       cccATGACAAGGTTGATCTTCCATCTGAATGTTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcg                                                                              >  1:228765/1‑100 (MQ=255)

                       cccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcg                                                                              >  2:2520655/1‑100 (MQ=255)

                       cccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcg                                                                              >  1:39642/1‑100 (MQ=255)

                          aTGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCgtgt                                                                           <  2:68416/100‑1 (MQ=255)

                            gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTgg                                                                         <  1:847854/100‑1 (MQ=255)

                             aCAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGc                                                                        >  2:267766/1‑100 (MQ=255)

                              cAAGGTTGATATTCCCTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCg                                                                       <  1:492521/100‑1 (MQ=255)

                                aGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCggg                                                                     >  2:2759777/1‑100 (MQ=255)

                                aGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCggg                                                                     <  2:1793611/100‑1 (MQ=255)

                                 ggTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGa                                                                    <  2:2849441/100‑1 (MQ=255)

                                 ggTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGa                                                                    <  1:4365583/100‑1 (MQ=255)

                                  gTTGATATTCGATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCGGGTAGTGTTGGATCTGAGCCAACAAGGCAACGTGATCGTGTGGCGGGa                                                                    <  1:937108/99‑1 (MQ=255)

                                   ttGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATc                                                                  >  2:2653963/1‑100 (MQ=255)

                                    tGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCa                                                                 >  1:1132111/1‑100 (MQ=255)

                                    tGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCa                                                                 >  1:4301800/1‑100 (MQ=255)

                                    tGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCa                                                                 <  1:1917037/100‑1 (MQ=255)

                                     gATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAc                                                                >  1:2538075/1‑100 (MQ=255)

                                      atatTCCATCGTAATTGTGGCGTGCTCAATCTTATAGTGATCCATCAGTTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAc                                                                <  2:1903564/99‑1 (MQ=255)

                                        atTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTg                                                             >  1:403483/1‑100 (MQ=255)

                                        atTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTg                                                             >  2:4134590/1‑100 (MQ=255)

                                         ttCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGc                                                            >  2:1101180/1‑100 (MQ=255)

                                             aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACAt                                                        <  1:2857458/100‑1 (MQ=255)

                                               cTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGc                                                      <  2:3259052/100‑1 (MQ=255)

                                                tGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                     <  1:4128162/100‑1 (MQ=255)

                                                 gAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCaa                                                    >  1:4403888/1‑100 (MQ=255)

                                                   aTGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACg                                                  <  2:3857432/100‑1 (MQ=255)

                                                    tggtggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGt                                                 <  2:2687306/100‑1 (MQ=255)

                                                      gtggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtca                                               <  2:3978843/100‑1 (MQ=255)

                                                      gtggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtca                                               <  1:3045752/100‑1 (MQ=255)

                                                      gtggCGCGGTCAATCTGAGAGTGAGCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATAACCTGCACATGCAACGtca                                               <  1:661699/100‑1 (MQ=255)

                                                       tggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcat                                              <  2:740250/100‑1 (MQ=255)

                                                       tggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcat                                              <  1:3260994/100‑1 (MQ=255)

                                                       tggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcat                                              <  1:1985247/100‑1 (MQ=255)

                                                         gCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcatca                                            <  2:666120/100‑1 (MQ=255)

                                                            tGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCg                                         <  1:3812198/100‑1 (MQ=255)

                                                             gCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCgg                                        >  1:1418164/1‑100 (MQ=255)

                                                              cTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGc                                       <  1:3991274/100‑1 (MQ=255)

                                                              cTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGc                                       <  1:2520643/100‑1 (MQ=255)

                                                               tCAATCTGATAGTGTTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGc                                       <  1:3718289/99‑1 (MQ=255)

                                                                  aTCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTtct                                   >  2:3696853/1‑100 (MQ=255)

                                                                    cTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCg                                 <  1:3868080/100‑1 (MQ=255)

                                                                     tGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGc                                >  1:2761816/1‑100 (MQ=255)

                                                                     tGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGc                                <  1:848564/100‑1 (MQ=255)

                                                                      gATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGcc                               <  2:2761831/100‑1 (MQ=255)

                                                                       aTAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCt                              <  2:3500184/100‑1 (MQ=255)

                                                                        tAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTa                             >  1:1822366/1‑100 (MQ=255)

                                                                         aGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                            <  1:2653950/100‑1 (MQ=255)

                                                                         aGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                            <  2:1132115/100‑1 (MQ=255)

                                                                         aGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                            <  1:892888/100‑1 (MQ=255)

                                                                          gTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAcc                           >  2:3698358/1‑100 (MQ=255)

                                                                           tGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAcc                           <  1:1661874/99‑1 (MQ=255)

                                                                             aTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatc                        >  1:526493/1‑100 (MQ=255)

                                                                              tCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcc                       >  1:1921994/1‑100 (MQ=255)

                                                                              tCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcc                       <  2:1582699/100‑1 (MQ=255)

                                                                               ccATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcca                      >  1:1403878/1‑100 (MQ=255)

                                                                                 aTCAGGTAGTGTTGGTTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAta                    <  2:4329230/100‑1 (MQ=255)

                                                                                 aTCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAta                    >  2:34410/1‑100 (MQ=255)

                                                                                      gTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgt               >  1:3814279/1‑100 (MQ=255)

                                                                                       ccttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGACTTCTCGCCTACCATCCAtacatgta              >  2:3009569/4‑100 (MQ=255)

                                                                                        cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtac             >  2:1783047/3‑100 (MQ=255)

                                                                                        cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtac             >  1:3412018/3‑100 (MQ=255)

                                                                                        cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCAGGATCACCTGCACATGCAAAGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtac             >  1:2394254/3‑100 (MQ=255)

                                                                                        aGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtac             <  1:3995344/100‑1 (MQ=255)

                                                                                         ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtaca            >  1:2783951/2‑100 (MQ=255)

                                                                                         gtgtTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtaca            <  1:3175677/100‑1 (MQ=255)

                                                                                             tGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGt        >  2:2129531/1‑100 (MQ=255)

                                                                                              ggATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTCCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTg       >  1:1756145/1‑100 (MQ=255)

                                                                                                aTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGca     <  1:2267639/100‑1 (MQ=255)

                                                                                                 tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGcac    >  1:92127/1‑100 (MQ=255)

CCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACAT  >  NC_000913/783913‑784111

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