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.7%	1140.6 / ‑4.2	362	coding (380/486 nt)	b4566	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (203/158); total (203/159)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 392 total aligned reads displayed.

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCA‑GACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCa                                                                                          <  1:2971416/93‑1 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCA‑g                                                                                        <  2:1202547/93‑1 (MQ=255)

  cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCA‑gg                                                                                       >  2:4086115/1‑92 (MQ=255)

   tGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCA‑gga                                                                                      >  1:6444214/1‑91 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGa                                                                                   >  1:6145021/1‑93 (MQ=255)

       ccGCAAACTCAATCTCTTTACTGGTTTTGAAGCCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGa                                                                                   <  1:3476789/93‑1 (MQ=255)

          cAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATaa                                                                                >  2:2391825/1‑93 (MQ=255)

           aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATaaa                                                                               >  1:7833500/1‑93 (MQ=255)

           aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATaaa                                                                               >  1:7150561/1‑93 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAg                                                                              >  2:4704844/1‑93 (MQ=255)

                gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAAt                                                                          >  2:618556/1‑93 (MQ=255)

                gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAAt                                                                          >  2:6314006/1‑93 (MQ=255)

                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAAtt                                                                         >  2:5454981/1‑93 (MQ=255)

                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAAtt                                                                         >  2:6033417/1‑93 (MQ=255)

                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAAtt                                                                         >  1:474784/1‑93 (MQ=255)

                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAAAACG‑CCGCCAGGACGGATAAAGAAAtt                                                                         >  2:2775474/1‑93 (MQ=255)

                   tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTcc                                                                       >  2:3584520/1‑93 (MQ=255)

                   tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTcc                                                                       <  1:1018525/93‑1 (MQ=255)

                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTccc                                                                      <  1:1497250/93‑1 (MQ=255)

                      ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCtt                                                                    >  2:5903860/1‑93 (MQ=255)

                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   >  1:5241052/1‑93 (MQ=255)

                        ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCTCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCa                                                                  >  1:3758434/1‑93 (MQ=255)

                        ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCa                                                                  <  2:807367/93‑1 (MQ=255)

                         tACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCAc                                                                 >  2:7221260/1‑93 (MQ=255)

                         tACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCAc                                                                 >  1:4799056/1‑93 (MQ=255)

                          aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCAcc                                                                <  2:1412085/93‑1 (MQ=255)

                          aCTGGTATTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCAcc                                                                <  1:386935/93‑1 (MQ=255)

                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCt                                                               >  2:7476789/1‑93 (MQ=255)

                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                              >  2:2276985/1‑93 (MQ=255)

                              tttttAAAGGAATCAAGAGGAACATGGCCAACTTAACCCCCCCCG‑CCTAAAACCCGCCCGCAG‑GACGGATAAAGAAATTCCCTTCACCTGAAg                                                           <  1:5536458/92‑1 (MQ=255)

                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                                            <  2:4467392/93‑1 (MQ=255)

                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                                            >  1:6796096/1‑93 (MQ=255)

                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                                            <  1:7725683/93‑1 (MQ=255)

                                  tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCga                                                        <  2:3112771/93‑1 (MQ=255)

                                   gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgag                                                       >  1:5068640/1‑93 (MQ=255)

                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                                      >  1:6938269/1‑93 (MQ=255)

                                      ggCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAtt                                                    <  1:2321414/93‑1 (MQ=255)

                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 <  1:3108708/93‑1 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:897221/1‑93 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:181651/1‑93 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                <  1:5901283/93‑1 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  2:3449729/1‑93 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:7410395/1‑93 (MQ=255)

                                           ccAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACg                                               >  1:2458154/1‑93 (MQ=255)

                                            cAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGc                                              >  1:568809/1‑93 (MQ=255)

                                               acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAc                                           <  2:3849494/93‑1 (MQ=255)

                                                cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAcc                                          >  1:6258520/1‑93 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  1:3192118/1‑93 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  1:4307628/1‑93 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  1:5129020/1‑93 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  1:5337463/1‑93 (MQ=255)

                                                   acgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTc                                       >  2:3753509/1‑93 (MQ=255)

                                                    cgaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                      <  1:2906944/93‑1 (MQ=255)

                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGTATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  2:2989279/1‑93 (MQ=255)

                                                       tGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                   <  2:5905429/93‑1 (MQ=255)

                                                       tGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                   >  1:4071966/1‑93 (MQ=255)

                                                       tGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                   <  1:1356940/93‑1 (MQ=255)

                                                       tGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                                   <  2:3358234/93‑1 (MQ=255)

                                                          cGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAGGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAc                                <  2:705192/93‑1 (MQ=255)

                                                          cGACTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAc                                <  1:1195216/93‑1 (MQ=255)

                                                            aCTTTATCTCCCGCGACATAAAAAACG‑CCGCCAGGACGGCTAAAGAAATTCCCTTCGGCAGAAGCGAGACCTAACGCCACCTTCCGCTGACCg                              >  2:2257895/1‑93 (MQ=255)

                                                             cTTTATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGt                             >  2:6389301/1‑93 (MQ=255)

                                                              tttATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            >  2:1669025/1‑93 (MQ=255)

                                                              tttATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            >  2:5291266/1‑93 (MQ=255)

                                                               ttATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           >  1:6673267/1‑93 (MQ=255)

                                                                tATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:1476310/93‑1 (MQ=255)

                                                                tATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:5642263/93‑1 (MQ=255)

                                                                tATCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          >  2:5422920/1‑93 (MQ=255)

                                                                 aTCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         <  1:3565781/93‑1 (MQ=255)

                                                                 aTCTCCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         >  2:4423222/1‑93 (MQ=255)

                                                                   ctcCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGa                       >  1:2500287/1‑93 (MQ=255)

                                                                   cccccGCGACATAAAGCCCG‑CCGCCAGGGCGGATAAAGAAATTCCCTTCACCTGACGCGAGATTTACCGCCACCTTCGGCTAACCGTTTTCGa                       <  2:4342720/91‑1 (MQ=255)

                                                                    tcCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  1:1115662/93‑1 (MQ=255)

                                                                    tcCCGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  1:1899067/93‑1 (MQ=255)

                                                                    ccccGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATGCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  1:4433941/92‑1 (MQ=255)

                                                                     cccGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATg                     <  1:2749652/93‑1 (MQ=255)

                                                                      ccGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt                    <  2:2070741/93‑1 (MQ=255)

                                                                      ccGCGACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt                    <  1:493714/93‑1 (MQ=255)

                                                                       cgcgACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTa                   <  2:4549580/93‑1 (MQ=255)

                                                                        gcgACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  <  1:1735462/93‑1 (MQ=255)

                                                                        gcgACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  <  1:6895517/93‑1 (MQ=255)

                                                                          gACATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                >  1:2602610/1‑93 (MQ=255)

                                                                            cATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              >  2:4617844/1‑93 (MQ=255)

                                                                            cATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              <  2:1981746/93‑1 (MQ=255)

                                                                            cATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              >  2:6686824/1‑93 (MQ=255)

                                                                            cATAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              <  1:449192/93‑1 (MQ=255)

                                                                             aTAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             >  2:1060423/1‑93 (MQ=255)

                                                                              tAAAACACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  1:4181192/1‑93 (MQ=255)

                                                                               aaaaCACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCGCCTTCTGCTGACCGTTTTCGATGTACAAACCGa           <  1:2162990/93‑1 (MQ=255)

                                                                                aaaCCCG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGCCCGTTTTCGATGTACAAACCGAg          <  2:429981/93‑1 (MQ=255)

                                                                                aaaCACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg          >  1:4561219/1‑93 (MQ=255)

                                                                                 aaCACG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAATCTAgc         >  2:684071/1‑93 (MQ=255)

                                                                                  acacG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcg        >  2:5748165/1‑93 (MQ=255)

                                                                                    acG‑CCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      >  1:5781050/1‑93 (MQ=255)

                                                                                      g‑ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc    <  1:6603993/93‑1 (MQ=255)

                                                                                      g‑ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc    <  2:2052989/93‑1 (MQ=255)

                                                                                      g‑ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc    >  2:1497094/1‑93 (MQ=255)

                                                                                        ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   >  1:2096141/1‑93 (MQ=255)

                                                                                        ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   >  2:7724550/1‑93 (MQ=255)

                                                                                         cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAt  <  2:5039347/93‑1 (MQ=255)

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACG‑CCGCCA‑GACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  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: `‑`