BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	115,788	A→G	C438R (TGC→CGC)	hofB ←	HofB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	115,788	0	A	G	99.3%	459.5 / NA	148	C438R (TGC→CGC)	hofB	HofB
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (86/61); minor base T (0/1); total (86/62)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAA  >  NC_000913/115698‑115877

ccaGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAggcggc                                                                                         >  2:2662810/1‑93 (MQ=255)

 caGCGCCAGAGTTGCTTTCTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGcc                                                                                        <  1:3875762/93‑1 (MQ=255)

  aGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCg                                                                                       >  1:884764/1‑93 (MQ=255)

   gctccAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATAAACTCTTCAAGGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGt                                                                                      <  2:2387807/90‑1 (MQ=255)

    cgcCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGtt                                                                                     <  2:3554679/93‑1 (MQ=255)

     gcCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGttt                                                                                    <  1:3598917/93‑1 (MQ=255)

     gcCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGttt                                                                                    <  1:4118704/93‑1 (MQ=255)

      ccAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGtttt                                                                                   >  2:4085414/1‑93 (MQ=255)

       cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTc                                                                                  >  2:359374/1‑93 (MQ=255)

       cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTc                                                                                  <  1:905994/93‑1 (MQ=255)

         gagTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaa                                                                                >  2:879692/1‑93 (MQ=255)

           gTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaa                                                                              <  1:3147149/93‑1 (MQ=255)

           gTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaa                                                                              >  1:3991104/1‑93 (MQ=255)

            ttGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaaa                                                                             <  2:3046162/93‑1 (MQ=255)

            ttGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaaa                                                                             <  2:1949310/93‑1 (MQ=255)

            ttGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCaaaaa                                                                             <  1:1073421/93‑1 (MQ=255)

                ttACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGt                                                                         >  1:1846204/1‑93 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  2:3246955/93‑1 (MQ=255)

                 tACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTa                                                                        <  1:1328294/93‑1 (MQ=255)

                   cTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCTTACg                                                                      <  2:666341/93‑1 (MQ=255)

                   cTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTCGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACg                                                                      >  2:482147/1‑93 (MQ=255)

                    tCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGc                                                                     >  2:1921741/1‑93 (MQ=255)

                    tCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGc                                                                     >  1:4117029/1‑93 (MQ=255)

                     cGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCa                                                                    <  1:3733670/93‑1 (MQ=255)

                     cGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCa                                                                    <  1:3892538/93‑1 (MQ=255)

                      gCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGACGGCCGTTTTCAAAAAGCGTACGCAt                                                                   >  1:4526077/1‑93 (MQ=255)

                         aTGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAcc                                                                >  2:1452621/1‑93 (MQ=255)

                         aTGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAcc                                                                <  1:3455044/93‑1 (MQ=255)

                         aTGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAcc                                                                <  1:2279046/93‑1 (MQ=255)

                          tGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATAccc                                                               >  1:52536/1‑93 (MQ=255)

                             ggCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCGGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGcc                                                            <  1:4609633/93‑1 (MQ=255)

                              gCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCt                                                           >  2:801693/1‑93 (MQ=255)

                              gCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCt                                                           >  1:2483280/1‑93 (MQ=255)

                                aTCCCCAGTACGCGGATTAGCTCTTCAACGGCGGTTAAGCCTTGCTCTACGGCCACGCGGCCGTTGTCAAAAAGCCTACGCATACCCGCCTGt                                                         >  2:1637297/1‑93 (MQ=255)

                                aTCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGt                                                         <  2:1827844/93‑1 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        >  2:2820480/1‑93 (MQ=255)

                                 tCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTc                                                        >  1:1561946/1‑93 (MQ=255)

                                    ccAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtg                                                     >  1:4174331/1‑93 (MQ=255)

                                       gTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgt                                                  >  1:4249249/1‑93 (MQ=255)

                                       gTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgt                                                  >  2:715042/1‑93 (MQ=255)

                                       gTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgg                                                  >  2:3886920/1‑92 (MQ=255)

                                        tACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgcgtg                                                 >  2:1677006/1‑93 (MQ=255)

                                           gcgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTgcgtgcgt                                              >  2:2256129/1‑93 (MQ=255)

                                            cgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGtt                                             <  1:72653/93‑1 (MQ=255)

                                            cgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGtt                                             <  1:38182/93‑1 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            >  2:4483917/1‑93 (MQ=255)

                                             ggATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGttt                                            >  1:2612871/1‑93 (MQ=255)

                                              gATTACCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTc                                           <  2:3790888/93‑1 (MQ=255)

                                              gATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTc                                           <  2:3705518/93‑1 (MQ=255)

                                              gATTAACTCGTCAAAGGTGGTTAAGCCTGGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTc                                           <  1:1117650/93‑1 (MQ=255)

                                                ttAACTCTTCAAAGGTGGTTAAGCTTTGCTCTACGGCCAGGCGGCCTTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCa                                         <  1:1533682/93‑1 (MQ=255)

                                                ttAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCa                                         >  1:2897608/1‑93 (MQ=255)

                                                 tAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAg                                        >  1:4292121/1‑93 (MQ=255)

                                                  aaCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGc                                       <  2:2144548/93‑1 (MQ=255)

                                                  aaCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGc                                       <  2:1746143/93‑1 (MQ=255)

                                                  aaCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCAGCCTGTCGTGCGTGCGTTTCCAGc                                       >  2:1832478/1‑93 (MQ=255)

                                                   aCTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCg                                      >  2:1713098/1‑93 (MQ=255)

                                                    ctctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCATTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGa                                     >  1:2519689/1‑93 (MQ=255)

                                                    ctctACAAGGGTGTTAACGCTTTGCTCTGCGGCAAGGCGCCCGTTTTAAAAAACCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGa                                     <  2:2514805/93‑1 (MQ=255)

                                                     tctTCAAAGGTGTTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCATGCGTGCGTTGCCAGCGAt                                    <  2:4574321/93‑1 (MQ=255)

                                                     tctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAt                                    >  1:1777223/1‑93 (MQ=255)

                                                     tctTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAt                                    >  2:1282820/1‑93 (MQ=255)

                                                      ctTCAAAGGTGGTTAAGCCTTGCTCTACGCCCAGGCGGCCGTTTTCAAAAAGCGTACCCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGAtt                                   <  1:221200/93‑1 (MQ=255)

                                                       ttCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                  >  1:1501498/1‑93 (MQ=255)

                                                       ttCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                  >  1:700905/1‑93 (MQ=255)

                                                       ttCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGCCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTc                                  <  2:4300503/93‑1 (MQ=255)

                                                        tCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCa                                 >  1:2800221/1‑93 (MQ=255)

                                                         cAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCaa                                <  2:2538246/93‑1 (MQ=255)

                                                         cAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCaa                                >  2:440933/1‑93 (MQ=255)

                                                          aaaGGTGTTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                               <  2:2533882/93‑1 (MQ=255)

                                                          aaaGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                               >  2:4464361/1‑93 (MQ=255)

                                                          aaaGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAAc                               <  2:2117874/93‑1 (MQ=255)

                                                            aGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGt                             >  2:2007896/1‑93 (MQ=255)

                                                            aGGTGGTTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGt                             <  2:4551254/93‑1 (MQ=255)

                                                                ggtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGt                         >  2:2385748/1‑93 (MQ=255)

                                                                 gtTAAGCCTTGCTCTCGGTCCAGGTGACGTTTTTAAAAAAGCCTACGCATTCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  2:4250122/93‑1 (MQ=255)

                                                                 gtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAGCGTCGGTa                        >  2:2198031/1‑93 (MQ=255)

                                                                 gtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  1:3361251/93‑1 (MQ=255)

                                                                 gtTAAGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  2:3472544/93‑1 (MQ=255)

                                                                 gtTAAGCCTTCCTTTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACCCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTa                        <  2:3430876/93‑1 (MQ=255)

                                                                  taaaGCGTTTCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAt                       <  2:1109185/91‑1 (MQ=255)

                                                                   tAAGCTTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTAtt                      <  1:55847/93‑1 (MQ=255)

                                                                    aaGCTTTGCTCTGCGCCCGGGCGGCCGTTTTCAAAAAGCGTACCCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTa                     <  1:916223/93‑1 (MQ=255)

                                                                    aaGCCTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTa                     >  2:1215061/1‑93 (MQ=255)

                                                                        cTTGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCgg                 <  1:4579434/93‑1 (MQ=255)

                                                                         ttGCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGa                >  1:3055986/1‑93 (MQ=255)

                                                                           tctctAGGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              <  2:2940872/92‑1 (MQ=255)

                                                                           tctctACGCCCAGGCGGCTGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              <  2:1952631/92‑1 (MQ=255)

                                                                           gCTCTACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGaaa              >  1:3094058/1‑93 (MQ=255)

                                                                            ctctACGGCCAGGCGGCGTTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAAt             <  1:1073611/93‑1 (MQ=255)

                                                                             tctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATa            <  1:935125/93‑1 (MQ=255)

                                                                             tctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATa            >  1:1302769/1‑93 (MQ=255)

                                                                             tctACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATa            >  1:390738/1‑93 (MQ=255)

                                                                              tgaCGGCCAGGCGGCCTTTTACAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATaa           <  2:2924593/91‑1 (MQ=255)

                                                                               tACGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAg          >  2:4356598/1‑93 (MQ=255)

                                                                                 cGGCCAGGTGGCCGTTTTCTAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt        >  2:4585965/1‑93 (MQ=255)

                                                                                 cGGCCAGGCGGCCGTTTTCAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt        <  1:3737801/93‑1 (MQ=255)

                                                                                     cAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACg    >  2:392982/1‑93 (MQ=255)

                                                                                     cAGGCGGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACg    >  2:32190/1‑93 (MQ=255)

                                                                                       ggcggcCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGat  >  1:3385342/1‑92 (MQ=255)

CCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAA  >  NC_000913/115698‑115877

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