BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,896,324	+G	coding (698/711 nt)	b3718 ←	putative isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,896,324	1	.	G	99.6%	931.4 / ‑4.9	250	V233A (GTG→GCG)	b3718	putative isomerase
Reads supporting (aligned to +/- strand): ref base . (0/0); major base G (113/136); minor base T (1/0); total (114/136)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 266 total aligned reads displayed.

TTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410

ttCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCa                                                                                         >  2:481864/1‑93 (MQ=255)

ttCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCa                                                                                         <  1:2489907/93‑1 (MQ=255)

ttCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCa                                                                                         >  1:4559455/1‑93 (MQ=255)

 tCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCag                                                                                        <  1:5066284/93‑2 (MQ=255)

  cGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagc                                                                                       <  2:1934228/93‑3 (MQ=255)

  cGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagc                                                                                       <  2:624972/93‑3 (MQ=255)

  cGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagc                                                                                       <  2:5316594/93‑3 (MQ=255)

   gTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGCGTGCTAATTCTGTGGCagag                                                                                      >  1:5253844/1‑90 (MQ=255)

   gTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagcg                                                                                      >  2:3173976/1‑90 (MQ=255)

      gTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCt                                                                                   <  1:2454463/93‑1 (MQ=255)

      gTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCt                                                                                   >  2:3746642/1‑93 (MQ=255)

       ttATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCtt                                                                                  <  1:4388228/93‑1 (MQ=255)

            tAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCa                                                                             >  1:1274962/1‑93 (MQ=255)

             aaaTAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAg                                                                            >  1:324469/1‑93 (MQ=255)

              aataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGc                                                                           >  1:2114675/1‑93 (MQ=255)

              aataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGc                                                                           <  2:173391/93‑1 (MQ=255)

               ataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCg                                                                          >  1:5478170/1‑93 (MQ=255)

                  atGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAtt                                                                       <  1:2039662/93‑1 (MQ=255)

                  atGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAtt                                                                       <  1:2823258/93‑1 (MQ=255)

                  atGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAtt                                                                       >  2:1093839/1‑93 (MQ=255)

                   tGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTa                                                                      <  1:2615513/93‑1 (MQ=255)

                   tGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTa                                                                      >  1:2710174/1‑93 (MQ=255)

                        aCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATc                                                                 <  2:3037802/93‑1 (MQ=255)

                        aCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATc                                                                 >  2:4376820/1‑93 (MQ=255)

                          aGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAg                                                               <  2:144539/93‑1 (MQ=255)

                           ggCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGc                                                              >  2:3137657/1‑93 (MQ=255)

                            gCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCg                                                             >  2:3345648/1‑93 (MQ=255)

                            gCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCg                                                             <  2:4011142/93‑1 (MQ=255)

                              cGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCATCGGCTTTATCAGCGATTACCATCAGCGaa                                                           >  2:6612861/1‑93 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  2:252299/1‑93 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  2:4958243/1‑93 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  1:3855044/1‑93 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  1:3416943/1‑93 (MQ=255)

                                gAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAgg                                                         <  1:2031098/93‑1 (MQ=255)

                                 aGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGa                                                        >  1:4625978/1‑93 (MQ=255)

                                 aGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGa                                                        <  2:754783/93‑1 (MQ=255)

                                 aGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGa                                                        <  1:1457696/93‑1 (MQ=255)

                                  gCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAt                                                       <  2:1154826/93‑1 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                     <  1:987516/93‑1 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                     <  1:1025611/93‑1 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                     >  1:4844597/1‑93 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                     >  1:5167344/1‑93 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                     <  1:3169860/93‑1 (MQ=255)

                                     tAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTTCCATCAGCGAAGGATGCa                                                    <  1:6172011/93‑1 (MQ=255)

                                     tAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCa                                                    <  2:5354560/93‑1 (MQ=255)

                                     tAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCa                                                    >  1:3646309/1‑93 (MQ=255)

                                      aaTTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAg                                                   <  1:101097/93‑1 (MQ=255)

                                        ttCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCt                                                 <  2:2368219/93‑1 (MQ=255)

                                          cacacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGc                                               >  2:245508/1‑93 (MQ=255)

                                          cacacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGc                                               <  1:3484685/93‑1 (MQ=255)

                                            cacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCaa                                             <  1:1832070/93‑1 (MQ=255)

                                             acaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAac                                            <  2:5864587/93‑1 (MQ=255)

                                             acaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAac                                            <  1:1905292/93‑1 (MQ=255)

                                               aTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAacac                                          <  1:5268715/93‑1 (MQ=255)

                                                   ggCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGaa                                      <  2:6457371/93‑1 (MQ=255)

                                                     ccTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGc                                    >  2:505118/1‑93 (MQ=255)

                                                       ttATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTg                                  <  1:3115859/93‑1 (MQ=255)

                                                        tATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTgg                                 <  2:2994071/93‑1 (MQ=255)

                                                         aTTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGa                                >  1:293369/1‑93 (MQ=255)

                                                         aTTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGa                                <  1:4313543/93‑1 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATc                           >  2:5735374/1‑93 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATc                           >  2:4461603/1‑93 (MQ=255)

                                                               tAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCt                          <  2:382654/93‑1 (MQ=255)

                                                               tAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCt                          >  2:915346/1‑93 (MQ=255)

                                                                aaGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCtt                         <  1:120702/93‑1 (MQ=255)

                                                                aaGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCtt                         <  1:6579259/93‑1 (MQ=255)

                                                                aaGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCtt                         <  1:5699969/93‑1 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTc                        <  1:286065/93‑1 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTc                        <  2:1468815/93‑1 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTc                        <  1:6415123/93‑1 (MQ=255)

                                                                  gCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCg                       <  2:5263845/93‑1 (MQ=255)

                                                                   cTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCgg                      <  1:400848/93‑1 (MQ=255)

                                                                   cTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCgg                      >  1:6187416/1‑93 (MQ=255)

                                                                   cTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCgg                      >  2:4056005/1‑93 (MQ=255)

                                                                    ttAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGt                     <  2:2574224/93‑1 (MQ=255)

                                                                     tAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCCAAACGGAAGCTGGAACATCTTCGGTa                    >  1:1745169/1‑93 (MQ=255)

                                                                     tAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTa                    <  1:3839595/93‑1 (MQ=255)

                                                                     tAACCGAGTGC‑AATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTaa                   <  2:2264210/93‑1 (MQ=255)

                                                                       aCCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTGCCATCAGCGAAGGATGCAGCTGCAACGCGGAAGCTGGAACATCTTCGGTAAc                  <  1:3302755/93‑1 (MQ=255)

                                                                       aCCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAc                  <  1:5724006/93‑1 (MQ=255)

                                                                       aCCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAc                  <  1:4967808/93‑1 (MQ=255)

                                                                        ccGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAcc                 <  2:4873143/93‑1 (MQ=255)

                                                                          gAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               >  1:612086/1‑93 (MQ=255)

                                                                          gAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               >  1:3804844/1‑93 (MQ=255)

                                                                          gAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               <  2:4597810/93‑1 (MQ=255)

                                                                          gAGTGCTAATGCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               <  2:3220873/93‑1 (MQ=255)

                                                                           aGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGGAACACGGAAGCTGGCACATCTTCGGTAACCGGa              >  1:4195334/1‑93 (MQ=255)

                                                                           aGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGa              <  2:5623997/93‑1 (MQ=255)

                                                                            gTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAc             <  2:1513834/93‑1 (MQ=255)

                                                                             tGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAcc            <  2:4969078/93‑1 (MQ=255)

                                                                             tGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAcc            >  2:2151375/1‑93 (MQ=255)

                                                                             tGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAcc            <  2:5551712/93‑1 (MQ=255)

                                                                             tGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAcc            >  1:6057507/1‑93 (MQ=255)

                                                                              gCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCt           >  2:420372/1‑93 (MQ=255)

                                                                               cTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCtt          <  1:3626492/93‑1 (MQ=255)

                                                                               cTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCtt          <  2:4434216/93‑1 (MQ=255)

                                                                               cTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCtt          <  1:2324222/93‑1 (MQ=255)

                                                                                tAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTg         <  2:5996091/93‑1 (MQ=255)

                                                                                tAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTg         >  2:2948587/1‑93 (MQ=255)

                                                                                tAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTg         <  1:2285017/93‑1 (MQ=255)

TTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410

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