BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,998,290	+G	coding (379/486 nt)	b4566 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,998,289	1	.	G	99.5%	653.2 / ‑4.2	201	coding (380/486 nt)	b4566	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (90/110); total (90/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 225 total aligned reads displayed.

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑G‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTAATCTCCCGCGACATAAAACACGCCG‑CCa                                                                                           >  1:3635753/1‑93 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑g                                                                                         <  1:4129731/93‑1 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑g                                                                                         <  1:1409562/93‑1 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑g                                                                                         >  2:2564221/1‑93 (MQ=255)

  cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGACG‑CCA‑gg                                                                                        >  1:3011428/1‑92 (MQ=255)

   tGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑gga                                                                                       <  2:2511824/93‑3 (MQ=255)

   tGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑gga                                                                                       <  1:1258242/93‑3 (MQ=255)

    gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑ggac                                                                                      >  1:2879678/1‑90 (MQ=255)

    gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑ggac                                                                                      <  1:418724/93‑4 (MQ=255)

     cACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑ggacg                                                                                     <  2:3361321/93‑5 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGGCATAAAACACGCCG‑CCAGG‑ACGGa                                                                                   <  1:3784410/93‑1 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGa                                                                                   <  1:1603425/93‑1 (MQ=255)

        cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGAt                                                                                  <  2:3017052/93‑1 (MQ=255)

         gCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATa                                                                                 >  1:4544272/1‑93 (MQ=255)

         gCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATa                                                                                 <  2:4339779/93‑1 (MQ=255)

         gCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATa                                                                                 <  2:777276/93‑1 (MQ=255)

          cAACCTGACTCTCTTTCCTGGTTTAGAAGGCATCCAGACGAACGATGCTAACTTTATCTCCCGCGACAAAAACCACGCCG‑CCAGG‑ACGGATaa                                                                                <  2:4243705/93‑1 (MQ=255)

           aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATaaa                                                                               <  1:4489359/93‑1 (MQ=255)

             aCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGa                                                                             <  2:1600303/93‑1 (MQ=255)

             aCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGa                                                                             >  2:3878176/1‑93 (MQ=255)

              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGaa                                                                            <  1:4450074/93‑1 (MQ=255)

               tGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑TCGGATAAAGaaa                                                                           >  2:2935187/1‑93 (MQ=255)

               tGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGaaa                                                                           >  1:2945950/1‑93 (MQ=255)

                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATACAACACTCCGCCCG‑G‑ACGGACAAAGACAtc                                                                         >  2:2621262/1‑92 (MQ=255)

                  aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTc                                                                        <  2:2531690/93‑1 (MQ=255)

                   tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTcc                                                                       <  1:3774131/93‑1 (MQ=255)

                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTccc                                                                      <  2:4430413/93‑1 (MQ=255)

                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCt                                                                     >  2:1971582/1‑93 (MQ=255)

                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCt                                                                     <  1:3016707/93‑1 (MQ=255)

                      ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCtt                                                                    <  1:2901745/93‑1 (MQ=255)

                      ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCtt                                                                    >  2:2606063/1‑93 (MQ=255)

                      ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCtt                                                                    <  1:4140713/93‑1 (MQ=255)

                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTc                                                                   <  1:4714952/93‑1 (MQ=255)

                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCt                                                               <  1:2937584/93‑1 (MQ=255)

                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCt                                                               >  2:3579187/1‑93 (MQ=255)

                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTg                                                              >  2:4527521/1‑93 (MQ=255)

                             ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGa                                                             <  1:4657010/93‑1 (MQ=255)

                                tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGc                                                          <  1:1292123/93‑1 (MQ=255)

                                tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGc                                                          <  2:4075527/93‑1 (MQ=255)

                                 ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCg                                                         <  2:484370/93‑1 (MQ=255)

                                   gAAGGCATCGAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCgag                                                       >  2:3143203/1‑93 (MQ=255)

                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                                      <  1:1499657/93‑1 (MQ=255)

                                     aGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAt                                                     >  1:3456809/1‑93 (MQ=255)

                                     aGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAt                                                     >  1:1964336/1‑93 (MQ=255)

                                      ggCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAtt                                                    >  1:3111354/1‑93 (MQ=255)

                                       gCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAttt                                                   >  2:1791670/1‑93 (MQ=255)

                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 >  2:1898114/1‑93 (MQ=255)

                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 >  2:2023141/1‑93 (MQ=255)

                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 <  1:1000989/93‑1 (MQ=255)

                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 >  2:1421369/1‑93 (MQ=255)

                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 <  1:396383/93‑1 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:4089041/1‑93 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:2145818/1‑93 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                <  1:2342378/93‑1 (MQ=255)

                                           ccAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACg                                               >  1:2447391/1‑93 (MQ=255)

                                           ccAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACg                                               <  2:2448150/93‑1 (MQ=255)

                                             aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                                             >  1:1734909/1‑93 (MQ=255)

                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAGAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  1:629557/1‑93 (MQ=255)

                                                cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAcc                                          <  2:2650002/93‑1 (MQ=255)

                                                 gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         <  2:4013043/93‑1 (MQ=255)

                                                 gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         >  2:4652495/1‑93 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  2:1292420/1‑93 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  2:2718756/1‑93 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  1:1436093/1‑93 (MQ=255)

                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  2:219016/1‑93 (MQ=255)

                                                      acgCCGACTTTAGCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAATAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    <  1:3591824/91‑1 (MQ=255)

                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    >  1:339330/1‑93 (MQ=255)

                                                       tGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                   <  2:813768/93‑1 (MQ=255)

                                                           caCTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAcc                               <  1:1358401/92‑1 (MQ=255)

                                                             cTTTATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGt                             <  2:1940103/93‑1 (MQ=255)

                                                              tttATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            >  2:3831348/1‑93 (MQ=255)

                                                               ttATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           >  1:268528/1‑93 (MQ=255)

                                                               ttATCTCCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           <  2:226012/93‑1 (MQ=255)

                                                                  tctcCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCg                        >  2:3299796/1‑93 (MQ=255)

                                                                    tcCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  1:3909717/93‑1 (MQ=255)

                                                                    tcCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      >  2:3977337/1‑93 (MQ=255)

                                                                    tcCCGCGACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  2:3400746/93‑1 (MQ=255)

                                                                       cgcgACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTa                   >  2:204688/1‑93 (MQ=255)

                                                                        gcgACATAAAACACGCCT‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  <  2:4016080/93‑1 (MQ=255)

                                                                        ccgACATAAAACACGCCG‑CCA‑GTACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  <  1:189511/92‑1 (MQ=255)

                                                                         cgACATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACa                 >  2:2203931/1‑93 (MQ=255)

                                                                           aCATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               >  1:2591125/1‑93 (MQ=255)

                                                                           aCATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               >  2:1700493/1‑93 (MQ=255)

                                                                           aCATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               >  2:1188923/1‑93 (MQ=255)

                                                                            cATAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              <  2:1129335/93‑1 (MQ=255)

                                                                             aTAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             >  2:893662/1‑93 (MQ=255)

                                                                             aTAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             <  1:4390180/93‑1 (MQ=255)

                                                                              tAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            <  2:3046298/93‑1 (MQ=255)

                                                                              tAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            <  2:4164004/93‑1 (MQ=255)

                                                                              tAAAACACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            <  2:2244739/93‑1 (MQ=255)

                                                                               aaaaCACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGa           <  2:290578/93‑1 (MQ=255)

                                                                               aaaaCACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGa           <  1:2639739/93‑1 (MQ=255)

                                                                                aaaCACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg          >  2:3588722/1‑93 (MQ=255)

                                                                                aaaCACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg          >  1:2530727/1‑93 (MQ=255)

                                                                                aaaCACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg          >  1:2254241/1‑93 (MQ=255)

                                                                                 aaCACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc         <  2:35073/93‑1 (MQ=255)

                                                                                 aaCACGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc         >  1:968236/1‑93 (MQ=255)

                                                                                   cacGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcgg       <  1:1196776/93‑1 (MQ=255)

                                                                                    acGCCG‑CCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      >  1:3041124/1‑93 (MQ=255)

                                                                                       ccg‑ccAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   <  2:474486/93‑1 (MQ=255)

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCG‑CCA‑G‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

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