BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,199,968	G→A	A155A (GCG→GCA)	hydH →	HydH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,199,968	0	G	A	86.3%	748.5 / 80.3	263	A155A (GCG→GCA)	hydH	HydH
Reads supporting (aligned to +/- strand): ref base G (16/20); new base A (122/105); total (138/125)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 268 total aligned reads displayed.

AGCCAGAGGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTATCGCCGTT  >  NC_000913/4199877‑4200058

  ccAGAGGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCAccc                                                                                         <  2:720503/93‑1 (MQ=255)

  ccAGAGGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCAccc                                                                                         <  1:650848/93‑1 (MQ=255)

  ccAGAGGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCAccc                                                                                         <  1:2372197/93‑1 (MQ=255)

   cAGAGGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCt                                                                                        <  2:5758739/93‑1 (MQ=255)

   cAGAGGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCt                                                                                        >  1:2226564/1‑93 (MQ=255)

     gagGAAAACTCCCGCTGTCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTgg                                                                                      <  2:22379/93‑1 (MQ=255)

     gagGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTgg                                                                                      >  1:5278332/1‑93 (MQ=255)

     gagGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTgg                                                                                      <  1:4492923/93‑1 (MQ=255)

     gagGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTgg                                                                                      >  1:3422038/1‑93 (MQ=255)

       ggAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTggcg                                                                                    <  1:1216929/93‑1 (MQ=255)

         aaaaCTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTggcggc                                                                                  <  2:3310850/93‑1 (MQ=255)

          aaaCTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCa                                                                                 <  2:2176901/93‑1 (MQ=255)

           aaCTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCAt                                                                                <  2:4390932/93‑1 (MQ=255)

              tCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGa                                                                             >  2:2943582/1‑93 (MQ=255)

               cccGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGaa                                                                            <  1:1718906/93‑1 (MQ=255)

                ccGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAAt                                                                           >  2:3170720/1‑93 (MQ=255)

                ccGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAAt                                                                           >  2:3855853/1‑93 (MQ=255)

                  gCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGc                                                                         <  2:2610880/93‑1 (MQ=255)

                  gCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGc                                                                         <  1:1601050/93‑1 (MQ=255)

                  gCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGc                                                                         <  2:4071933/93‑1 (MQ=255)

                    tGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGcac                                                                       >  2:2327035/1‑93 (MQ=255)

                    tGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGcac                                                                       >  1:2871234/1‑93 (MQ=255)

                        ggcTGCTTGGGAAAACGTAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCGCAata                                                                   <  1:6319904/93‑1 (MQ=255)

                        ggcTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAata                                                                   >  1:2133886/1‑93 (MQ=255)

                        ggcTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAata                                                                   >  2:3818051/1‑93 (MQ=255)

                         gctgctTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAatat                                                                  <  2:151235/93‑1 (MQ=255)

                          ctgctTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATg                                                                 >  2:2859030/1‑93 (MQ=255)

                            gctTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATgtc                                                               >  1:5401185/1‑91 (MQ=255)

                            gctTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGcc                                                               <  1:954819/93‑1 (MQ=255)

                            gctTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGcc                                                               <  1:3470016/93‑1 (MQ=255)

                             ctTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCcg                                                              >  2:5468451/1‑93 (MQ=255)

                               tGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCcgcg                                                            <  2:1125991/93‑1 (MQ=255)

                               tGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCcgcg                                                            <  1:5579189/93‑1 (MQ=255)

                                gggAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCgcgc                                                           >  1:1808691/1‑93 (MQ=255)

                                gggAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCgcgc                                                           <  2:102357/93‑1 (MQ=255)

                                gggAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCgcgc                                                           <  2:3901079/93‑1 (MQ=255)

                                 ggAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCt                                                          >  2:1830645/1‑93 (MQ=255)

                                 ggAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCt                                                          >  2:4278209/1‑93 (MQ=255)

                                   aaaaCGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGt                                                        <  1:1829614/93‑1 (MQ=255)

                                   aaaaCGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGt                                                        >  1:1967525/1‑93 (MQ=255)

                                   aaaaCGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGt                                                        <  1:5034762/93‑1 (MQ=255)

                                     aaCGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTaa                                                      <  2:4629083/93‑1 (MQ=255)

                                      aCGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAAc                                                     <  2:4039108/93‑1 (MQ=255)

                                      aCGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAAc                                                     <  2:5575305/93‑1 (MQ=255)

                                       cGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACg                                                    >  2:3734852/1‑93 (MQ=255)

                                       cGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACg                                                    <  2:1296823/93‑1 (MQ=255)

                                       cGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACg                                                    <  2:607375/93‑1 (MQ=255)

                                        ggAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACgg                                                   >  2:4250092/1‑93 (MQ=255)

                                        ggAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACgg                                                   >  2:2779044/1‑93 (MQ=255)

                                           aaCTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCaa                                                >  2:3435215/1‑93 (MQ=255)

                                           aaCTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCaa                                                >  2:2083329/1‑93 (MQ=255)

                                            aCTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCaaa                                               >  1:4632719/1‑93 (MQ=255)

                                            aCTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCaaa                                               <  2:2345308/93‑1 (MQ=255)

                                             cTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAg                                              <  2:4476901/93‑1 (MQ=255)

                                                cGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCtg                                           <  2:5629189/93‑1 (MQ=255)

                                                 gCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCtgt                                          >  1:1164630/1‑93 (MQ=255)

                                                  ccTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCtgtg                                         >  2:4559647/1‑93 (MQ=255)

                                                   cTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGc                                        >  2:6590897/1‑93 (MQ=255)

                                                   cTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGc                                        >  1:657987/1‑93 (MQ=255)

                                                    tGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGcc                                       <  1:2456963/93‑1 (MQ=255)

                                                    tGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGcc                                       <  1:1074989/93‑1 (MQ=255)

                                                    tGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGcc                                       <  1:4253316/93‑1 (MQ=255)

                                                    tGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGcc                                       >  1:2115394/1‑93 (MQ=255)

                                                       aCTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCcaca                                    <  2:1230778/93‑1 (MQ=255)

                                                         ttGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAg                                  >  2:5029839/1‑93 (MQ=255)

                                                           gAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTa                                <  1:4362459/93‑1 (MQ=255)

                                                            aGGTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAg                               <  1:3737794/93‑1 (MQ=255)

                                                            aGGTCTATCGTTTGTCCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAg                               <  2:2048432/93‑1 (MQ=255)

                                                              gTCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGGAGAt                             >  1:1303607/1‑93 (MQ=255)

                                                               tCTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATg                            >  2:3656623/1‑93 (MQ=255)

                                                                cTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGc                           >  1:3075291/1‑93 (MQ=255)

                                                                cTATCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGc                           >  1:2983303/1‑93 (MQ=255)

                                                                  aTCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGcac                         >  1:2687566/1‑93 (MQ=255)

                                                                   tCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGcaca                        >  2:6418090/1‑93 (MQ=255)

                                                                   tCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGcaca                        >  2:6447381/1‑93 (MQ=255)

                                                                   tCGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGcaca                        <  2:4753320/93‑1 (MQ=255)

                                                                    cGTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCacaa                       <  1:4686493/93‑1 (MQ=255)

                                                                     gTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCacaac                      >  1:5714582/1‑93 (MQ=255)

                                                                     gTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCacaac                      >  2:5684717/1‑93 (MQ=255)

                                                                     gTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCacaac                      <  2:1863021/93‑1 (MQ=255)

                                                                     gTTTGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTATGCCACAAGTAGATGCacaac                      <  1:4227462/93‑1 (MQ=255)

                                                                        tGTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAgg                   >  1:6281581/1‑93 (MQ=255)

                                                                         gTTCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGc                  <  1:3580273/93‑1 (MQ=255)

                                                                          ttCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCt                 <  2:198426/93‑1 (MQ=255)

                                                                          ttCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCt                 >  1:2752398/1‑93 (MQ=255)

                                                                          ttCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCt                 >  1:485399/1‑93 (MQ=255)

                                                                           tCCAGCCAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTa                <  1:2594484/93‑1 (MQ=255)

                                                                              aGCCAATTTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTAttt             >  2:2694351/1‑93 (MQ=255)

                                                                                ccAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTAttttt           >  1:73440/1‑93 (MQ=255)

                                                                                ccAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGGACAACAGGCTAttttt           >  1:593707/1‑93 (MQ=255)

                                                                                ccAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTAttttt           <  2:5916530/93‑1 (MQ=255)

                                                                                 cAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTa          <  2:1458867/93‑1 (MQ=255)

                                                                                 cAATGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTa          <  2:3600739/93‑1 (MQ=255)

                                                                                  aaTGTCAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTAt         >  2:874731/1‑93 (MQ=255)

                                                                                       cAGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTATCGCCg    <  2:6127027/93‑1 (MQ=255)

                                                                                        aGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTATCGCCGt   >  2:481443/1‑93 (MQ=255)

                                                                                        aGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTATCGCCGt   <  2:107169/93‑1 (MQ=255)

                                                                                        aGCACCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTATCGCCGt   <  2:905338/93‑1 (MQ=255)

                                                                                         gcgcCCTGGCGGCATGTAATGCACAAGATGCCGCGCTGTAACGGCAAAGCTGTGCCTCAAGTAGATGCACAACAGGCTATTTTTATCGCCGtt  <  2:5098254/93‑1 (MQ=255)

                                                                                         gcgcCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTATCGCCGtt  <  2:5287333/93‑1 (MQ=255)

AGCCAGAGGAAAACTCCCGCTGGCGGCTGCTTGGGAAAACGGAAACTACGCCTGCACTTGAGGTCTATCGTTTGTTCCAGCCAATGTCAGCGCCCTGGCGGCATGGAATGCACAATATGCCGCGCTGTAACGGCAAAGCTGTGCCACAAGTAGATGCACAACAGGCTATTTTTATCGCCGTT  >  NC_000913/4199877‑4200058

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