BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,754,708	G→A	L249L (CTC→CTT)	yiaY ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,754,708	0	G	A	83.5%	922.3 / 125.1	345	L249L (CTC→CTT)	yiaY	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (34/23); new base A (145/143); total (179/166)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 345 total aligned reads displayed.

TACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCAC  >  NC_000913/3754618‑3754797

tacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAg                                                                                         <  1:4149941/93‑1 (MQ=255)

tacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAg                                                                                         >  1:3006627/1‑93 (MQ=255)

tacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAg                                                                                         <  2:6866458/93‑1 (MQ=255)

tacaCCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAg                                                                                         >  2:6780455/1‑93 (MQ=255)

 acacCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAgg                                                                                        >  1:1181235/1‑93 (MQ=255)

 acacCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAgg                                                                                        <  1:5152190/93‑1 (MQ=255)

  cacCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGa                                                                                       >  2:4487368/1‑93 (MQ=255)

   acCGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGaa                                                                                      >  1:1617478/1‑93 (MQ=255)

    ccGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAAc                                                                                     <  1:6445193/93‑1 (MQ=255)

     cGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACt                                                                                    <  1:5371816/93‑1 (MQ=255)

     cGTGTGGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACt                                                                                    <  1:2904104/93‑1 (MQ=255)

        gtgGCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTggg                                                                                 <  1:6128230/93‑1 (MQ=255)

          ggCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCa                                                                               <  2:5297614/93‑1 (MQ=255)

           gCAGGCTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCAt                                                                              >  1:6264024/1‑93 (MQ=255)

               gctGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGc                                                                          <  2:6936217/91‑1 (MQ=255)

                  gTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCAt                                                                       <  2:7504935/93‑1 (MQ=255)

                  gTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCAt                                                                       >  1:2633697/1‑93 (MQ=255)

                   tAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGCCATAAGCCAtt                                                                      <  1:2538362/93‑1 (MQ=255)

                   tAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGAAACTGGGCATAAGCCAtt                                                                      >  1:6344263/1‑93 (MQ=255)

                    aGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTg                                                                     >  1:5638140/1‑93 (MQ=255)

                     gAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGc                                                                    <  2:3948218/93‑1 (MQ=255)

                     gAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGc                                                                    >  2:1321282/1‑93 (MQ=255)

                     gAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGc                                                                    >  2:1735208/1‑93 (MQ=255)

                       aaCCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCtt                                                                  <  1:2271797/93‑1 (MQ=255)

                        aCCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTc                                                                 <  2:1642833/93‑1 (MQ=255)

                        aCCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTc                                                                 <  1:6881720/93‑1 (MQ=255)

                         ccGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCa                                                                >  2:5672207/1‑93 (MQ=255)

                         ccGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCa                                                                >  2:6014403/1‑93 (MQ=255)

                            cccAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCAcgc                                                             <  2:5234631/93‑1 (MQ=255)

                             ccAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCAcgcg                                                            >  1:4382752/1‑93 (MQ=255)

                              cAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACgcgc                                                           >  2:6700873/1‑93 (MQ=255)

                              cAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACgcgc                                                           >  2:4210998/1‑93 (MQ=255)

                                 cTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCttt                                                        <  1:1387409/93‑1 (MQ=255)

                                    gTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGc                                                     >  2:2250512/1‑93 (MQ=255)

                                       cgcCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCAtt                                                  <  2:6284392/93‑1 (MQ=255)

                                       cgcCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCAtt                                                  <  2:5075269/93‑1 (MQ=255)

                                         ccATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTAc                                                <  1:1037296/93‑1 (MQ=255)

                                          cATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACt                                               <  2:5177043/93‑1 (MQ=255)

                                          cATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACt                                               >  2:5467532/1‑93 (MQ=255)

                                          cATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACt                                               <  1:4460988/93‑1 (MQ=255)

                                           aTCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTg                                              >  2:2895196/1‑93 (MQ=255)

                                           aTCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTg                                              >  2:1974405/1‑93 (MQ=255)

                                           aTCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTg                                              >  2:6548922/1‑93 (MQ=255)

                                            tCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGc                                             >  1:5507931/1‑93 (MQ=255)

                                            tCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGc                                             <  1:6691450/93‑1 (MQ=255)

                                              gCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCa                                           <  1:179872/93‑1 (MQ=255)

                                               cATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCAt                                          <  2:1045932/93‑1 (MQ=255)

                                               cATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCAt                                          <  1:1515201/93‑1 (MQ=255)

                                                aTGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATc                                         >  1:3355383/1‑93 (MQ=255)

                                                aTGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATc                                         >  1:4075549/1‑93 (MQ=255)

                                                aTGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATc                                         >  2:1814685/1‑93 (MQ=255)

                                                aTGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATc                                         <  2:3105194/93‑1 (MQ=255)

                                                  gAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCtt                                       >  1:3179090/1‑93 (MQ=255)

                                                   aaCATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTc                                      >  1:3538091/1‑93 (MQ=255)

                                                   aaCATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTc                                      <  2:6592372/93‑1 (MQ=255)

                                                    aCATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCa                                     >  1:3912421/1‑93 (MQ=255)

                                                     cATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCaa                                    <  2:150782/93‑1 (MQ=255)

                                                     cATAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCaa                                    <  2:1780463/93‑1 (MQ=255)

                                                      aTAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAAc                                   <  2:4221756/93‑1 (MQ=255)

                                                       tAACCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACg                                  >  1:4158604/1‑93 (MQ=255)

                                                        aaCCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCCTTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACgg                                 >  1:5785989/1‑93 (MQ=255)

                                                        aaCCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACgg                                 <  1:4408916/93‑1 (MQ=255)

                                                         aCCCAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGc                                >  1:855285/1‑93 (MQ=255)

                                                          cccAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCt                               >  1:4742648/1‑93 (MQ=255)

                                                          cccAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCt                               <  1:1129173/93‑1 (MQ=255)

                                                           ccAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTa                              >  2:2451392/1‑93 (MQ=255)

                                                           ccAGAGAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTa                              >  1:7266424/1‑93 (MQ=255)

                                                            cAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTaa                             <  1:6519065/93‑1 (MQ=255)

                                                             agagAAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAAc                            <  2:4582470/93‑1 (MQ=255)

                                                              gagaAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACg                           <  2:921596/93‑1 (MQ=255)

                                                                gaAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGc                         >  2:381001/1‑93 (MQ=255)

                                                                gaAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGc                         >  2:4045211/1‑93 (MQ=255)

                                                                gaAGCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCCTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACggt                         >  2:4921362/1‑92 (MQ=255)

                                                                   gCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAgg                      >  2:7148318/1‑93 (MQ=255)

                                                                   gCATTATTGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAgg                      >  2:5626287/1‑93 (MQ=255)

                                                                    cATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAggg                     >  1:6812199/1‑92 (MQ=255)

                                                                      ttattGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGttt                   >  1:3726066/1‑93 (MQ=255)

                                                                       tattGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGtttt                  <  2:555191/93‑1 (MQ=255)

                                                                        attGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTc                 <  2:3218473/93‑1 (MQ=255)

                                                                        attGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTc                 <  1:2066342/93‑1 (MQ=255)

                                                                        attGAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTc                 <  2:6174447/93‑1 (MQ=255)

                                                                           gAACGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGc              >  1:3976000/1‑93 (MQ=255)

                                                                            aaCGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCa             <  2:7423659/93‑1 (MQ=255)

                                                                            aaCGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGCTTTCGGCa             >  1:1958275/1‑93 (MQ=255)

                                                                             aCGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCaa            <  2:1879795/93‑1 (MQ=255)

                                                                              cGCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAAt           <  2:5849075/93‑1 (MQ=255)

                                                                               gCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATc          <  1:4966103/93‑1 (MQ=255)

                                                                               gCCATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATc          >  1:717550/1‑93 (MQ=255)

                                                                                 cATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCAt        <  2:3586796/93‑1 (MQ=255)

                                                                                 cATACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTGTTCGGCGATCAt        >  1:7487272/1‑93 (MQ=255)

                                                                                   tACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATgg      >  1:7052348/1‑93 (MQ=255)

                                                                                   tACCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATgg      >  2:1189015/1‑93 (MQ=255)

                                                                                    aCCGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATggg     >  1:5277563/1‑92 (MQ=255)

                                                                                     ccGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTc    >  2:2706255/1‑93 (MQ=255)

                                                                                      cGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCa   >  2:2649049/1‑93 (MQ=255)

                                                                                      cGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCa   <  2:3006512/93‑1 (MQ=255)

                                                                                      cGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCa   <  2:7024868/93‑1 (MQ=255)

                                                                                      cGGCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCa   <  1:1677473/93‑1 (MQ=255)

                                                                                       ggCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCAc  >  1:258571/1‑93 (MQ=255)

                                                                                       ggCAAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCAc  >  2:3311375/1‑93 (MQ=255)

TACACCGTGTGGCAGGTTGTAGAAACCGCCCAGCTGGTGCGCCATCGCATGAACATAACCCAGAGAAGCATTATTGAACGCCATACCGGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCAC  >  NC_000913/3754618‑3754797

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