BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,127,469	G→A	G254G (GGC→GGT)	b2981 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,127,469	0	G	A	99.5%	674.7 / NA	193	G254G (GGC→GGT)	b2981	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (119/73); minor base C (0/1); total (119/74)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

GTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAT  >  NC_000913/3127379‑3127558

gTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAcc                                                                                         >  1:1233821/1‑93 (MQ=255)

gTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAcc                                                                                         >  1:1608288/1‑93 (MQ=255)

 ttCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAccc                                                                                        <  2:306983/93‑1 (MQ=255)

  tCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCa                                                                                       <  1:3903415/93‑1 (MQ=255)

       aaGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTc                                                                                  >  2:1522665/1‑93 (MQ=255)

        aGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcg                                                                                 >  2:1984136/1‑93 (MQ=255)

         gCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgc                                                                                >  2:188891/1‑93 (MQ=255)

          ccGACAATAAACGCTCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgcg                                                                               >  1:2754061/1‑93 (MQ=255)

          ccGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCACATATTTAAGTTGATCACCCATGGTcgcg                                                                               >  1:3955472/1‑93 (MQ=255)

           cGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGa                                                                              >  1:3564101/1‑93 (MQ=255)

            gACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGaa                                                                             <  1:4338200/93‑1 (MQ=255)

              cAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAtt                                                                           <  2:695503/93‑1 (MQ=255)

               aaTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTc                                                                          <  1:459467/93‑1 (MQ=255)

               aaTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTc                                                                          >  1:3520800/1‑93 (MQ=255)

                aTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCa                                                                         >  1:3773583/1‑93 (MQ=255)

                 tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                        >  2:3469178/1‑93 (MQ=255)

                   aaCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGt                                                                      >  2:4652819/1‑93 (MQ=255)

                   aaCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGt                                                                      >  1:1420403/1‑93 (MQ=255)

                    aCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGtt                                                                     >  2:2132978/1‑93 (MQ=255)

                    aCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGtt                                                                     <  2:107065/93‑1 (MQ=255)

                    aCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGtt                                                                     <  1:4269893/93‑1 (MQ=255)

                     cgcgCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTc                                                                    <  2:2657085/93‑1 (MQ=255)

                     cgcgCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTc                                                                    <  2:3178562/93‑1 (MQ=255)

                      gcgcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCg                                                                   <  2:3988386/93‑1 (MQ=255)

                       cgcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCgg                                                                  >  1:2179855/1‑93 (MQ=255)

                        gcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGt                                                                 <  2:1924218/93‑1 (MQ=255)

                           gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                              >  2:1470957/1‑93 (MQ=255)

                             gTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACg                                                            >  1:2621368/1‑93 (MQ=255)

                                  gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       >  1:322812/1‑93 (MQ=255)

                                  gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       >  1:526802/1‑93 (MQ=255)

                                  gAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcaca                                                       >  1:2367092/1‑93 (MQ=255)

                                   aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                      >  1:3415143/1‑93 (MQ=255)

                                   aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                      >  2:1276723/1‑93 (MQ=255)

                                   aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                      >  2:4501125/1‑93 (MQ=255)

                                    atcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaaa                                                     <  2:1707509/93‑1 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                  >  2:1856046/1‑93 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                  >  2:2941689/1‑93 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 <  1:1948042/93‑1 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 >  2:2156589/1‑93 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 >  1:3395611/1‑93 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 <  2:4348884/93‑1 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 <  1:448356/93‑1 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACACATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 >  2:356590/1‑93 (MQ=255)

                                         tctcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcccc                                                >  1:2255831/1‑93 (MQ=255)

                                          ctcAGCGAGATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCa                                               <  2:181131/93‑1 (MQ=255)

                                           tcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGTTAAACGGCACAAATCCCCAg                                              >  2:2227264/1‑93 (MQ=255)

                                           tcAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAg                                              <  2:3352478/93‑1 (MQ=255)

                                               cGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCCTGGACGCGAATTCATGTTCGGTAAAAGGCACAAATCACAAGgtgt                                          >  2:926802/1‑93 (MQ=255)

                                                gATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgtg                                         >  1:608153/1‑93 (MQ=255)

                                                gATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgtg                                         >  1:1852191/1‑93 (MQ=255)

                                                 atatAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgc                                        >  2:489273/1‑93 (MQ=255)

                                                 atatAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgc                                        >  2:845618/1‑93 (MQ=255)

                                                 atatAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgc                                        >  2:2352910/1‑93 (MQ=255)

                                                 atatAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgc                                        <  1:2101285/93‑1 (MQ=255)

                                                 atatAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgc                                        <  1:2935142/93‑1 (MQ=255)

                                                  tataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgct                                       >  2:4316792/1‑93 (MQ=255)

                                                    taGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgctgc                                     >  1:1976889/1‑93 (MQ=255)

                                                      gatgatATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcca                                   <  1:4248572/93‑1 (MQ=255)

                                                      gATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcca                                   >  2:3925608/1‑93 (MQ=255)

                                                       aTCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccag                                  >  2:955772/1‑93 (MQ=255)

                                                        tCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagc                                 >  2:1860014/1‑93 (MQ=255)

                                                          atatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcca                               <  1:1672334/93‑1 (MQ=255)

                                                             tCGTCCGGCACCAGATATTTAAGTTGATTACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            >  2:2526693/1‑93 (MQ=255)

                                                             tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            >  1:1032695/1‑93 (MQ=255)

                                                             tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            >  1:1540701/1‑93 (MQ=255)

                                                             tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            >  1:1913753/1‑93 (MQ=255)

                                                              cGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGt                           <  1:2776260/93‑1 (MQ=255)

                                                              cGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGt                           <  1:128402/93‑1 (MQ=255)

                                                              cGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGt                           <  1:2596703/93‑1 (MQ=255)

                                                                tCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTg                         >  2:632597/1‑93 (MQ=255)

                                                                   ggCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaa                      <  1:1018817/93‑1 (MQ=255)

                                                                   ggCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaa                      >  1:1656265/1‑93 (MQ=255)

                                                                   ggCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaa                      <  1:663420/93‑1 (MQ=255)

                                                                    gCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaa                     >  1:2716267/1‑93 (MQ=255)

                                                                     cACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCGCCAGTTGTGCTGCCAGCCAGAGTTGAAAAt                    >  2:3614455/1‑93 (MQ=255)

                                                                      tccAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaagc                   <  2:590041/92‑3 (MQ=255)

                                                                      aCCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtc                   >  1:3736102/1‑93 (MQ=255)

                                                                       ccAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtct                  <  1:48981/93‑1 (MQ=255)

                                                                       ccAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtct                  >  1:4097445/1‑93 (MQ=255)

                                                                        cAGATATTTAAGTTGATCCCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                 <  1:1571728/93‑1 (MQ=255)

                                                                        cAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                 >  2:3833918/1‑93 (MQ=255)

                                                                        cAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                 >  2:974249/1‑93 (MQ=255)

                                                                        cAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                 >  1:872277/1‑93 (MQ=255)

                                                                         aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCa                >  2:4246632/1‑93 (MQ=255)

                                                                         aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCa                >  2:331704/1‑93 (MQ=255)

                                                                          gATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCAc               <  2:1146256/93‑1 (MQ=255)

                                                                          gATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCAc               <  1:4109628/93‑1 (MQ=255)

                                                                          gATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCAc               >  2:4463901/1‑93 (MQ=255)

                                                                          gATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCAc               <  2:898772/93‑1 (MQ=255)

                                                                          gATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCAc               >  2:2890491/1‑93 (MQ=255)

                                                                           atatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACg              <  1:2214012/93‑1 (MQ=255)

                                                                           atatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACg              >  2:2629671/1‑93 (MQ=255)

                                                                                 aaGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGAt        <  2:2287668/93‑1 (MQ=255)

                                                                                  aGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATc       <  2:957882/93‑1 (MQ=255)

                                                                                   gTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCt      <  2:1804567/93‑1 (MQ=255)

                                                                                     tGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGc    >  1:4202913/1‑93 (MQ=255)

                                                                                      gATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCa   >  1:3163827/1‑93 (MQ=255)

                                                                                       aTCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAt  <  1:3029032/93‑1 (MQ=255)

                                                                                       aTCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAt  <  1:717547/93‑1 (MQ=255)

                                                                                       aTCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAt  <  2:1641009/93‑1 (MQ=255)

GTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAT  >  NC_000913/3127379‑3127558

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