BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,369,263	A→G	V316A (GTG→GCG)	b3224 ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,369,263	0	A	G	97.2%	631.1 / 5.7	211	V316A (GTG→GCG)	b3224	sialic acid transporter
Reads supporting (aligned to +/- strand): ref base A (3/3); new base G (93/112); total (96/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 216 total aligned reads displayed.

CAAACGTACGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCACATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGAAAGCA  >  NC_000913/3369173‑3369352

 aaaCGTACGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCAt                                                                                        >  1:1981189/1‑93 (MQ=255)

   aCGTACGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTg                                                                                      <  2:2024143/93‑1 (MQ=255)

   aCGTACGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTg                                                                                      >  1:3985458/1‑93 (MQ=255)

      tACGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCt                                                                                   <  1:3881899/93‑1 (MQ=255)

       aCGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTa                                                                                  <  2:814184/93‑1 (MQ=255)

        cGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTAc                                                                                 <  2:51713/93‑1 (MQ=255)

         gCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTACAGAACAGCGCATTGGCTACa                                                                                >  1:4632720/1‑93 (MQ=255)

          cTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAg                                                                               >  1:877881/1‑93 (MQ=255)

          cTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAg                                                                               <  1:3250425/93‑1 (MQ=255)

            ttGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTa                                                                             <  2:1639012/93‑1 (MQ=255)

            ttGCGCGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTa                                                                             <  2:2328502/93‑1 (MQ=255)

              gCGGTTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATg                                                                           <  1:114883/93‑1 (MQ=255)

              gCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATg                                                                           <  1:961269/93‑1 (MQ=255)

                 ggTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCgg                                                                        >  2:2544387/1‑93 (MQ=255)

                 ggTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCgg                                                                        <  1:48150/93‑1 (MQ=255)

                  gTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCggg                                                                       >  1:1613564/1‑93 (MQ=255)

                  gTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCggg                                                                       <  1:2760542/93‑1 (MQ=255)

                     cccAGCCAGTCACCGAGGACGCCACCTCCGCAGCATCCCACCGCCGCGCCACAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTta                                                                    <  2:1725615/93‑1 (MQ=255)

                     cccAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTta                                                                    >  2:2417054/1‑93 (MQ=255)

                      ccagccagTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCACATTGGCTACAGTATGCGGGTtat                                                                   >  2:4150153/1‑93 (MQ=255)

                       cagccagTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTtata                                                                  <  2:3719680/93‑1 (MQ=255)

                        atccagTCACCGACGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATaa                                                                 <  2:2289577/91‑1 (MQ=255)

                        agccagTCCCCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATaa                                                                 <  1:1019776/93‑1 (MQ=255)

                        agccagTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATaa                                                                 <  1:2114637/93‑1 (MQ=255)

                         gccagTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCACATTGGCTACAGTATGCGGGTTATAAg                                                                >  1:3165080/1‑93 (MQ=255)

                          ccagTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGc                                                               >  2:3777703/1‑93 (MQ=255)

                           cagTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGcc                                                              >  1:2023974/1‑93 (MQ=255)

                             gTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAg                                                            >  1:665645/1‑93 (MQ=255)

                             gTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAg                                                            <  2:2124164/93‑1 (MQ=255)

                              tCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGa                                                           >  1:2881048/1‑93 (MQ=255)

                                aCCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATc                                                         >  1:4510357/1‑93 (MQ=255)

                                 ccGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCg                                                        >  1:3024254/1‑93 (MQ=255)

                                  cGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCgg                                                       >  1:3270132/1‑93 (MQ=255)

                                   gAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGt                                                      <  1:3788689/93‑1 (MQ=255)

                                    aGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGtt                                                     <  2:3416062/93‑1 (MQ=255)

                                    aGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGtt                                                     >  1:1685234/1‑93 (MQ=255)

                                     ggAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGttt                                                    <  2:2911039/93‑1 (MQ=255)

                                           cACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGAt                                              <  2:4628113/93‑1 (MQ=255)

                                            aCCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATa                                             <  1:4365685/93‑1 (MQ=255)

                                            aCCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATa                                             >  1:4189991/1‑93 (MQ=255)

                                            aCCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATa                                             <  1:2642563/93‑1 (MQ=255)

                                             ccTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATAc                                            <  1:3602067/93‑1 (MQ=255)

                                              cTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACg                                           <  1:3030985/93‑1 (MQ=255)

                                               tACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGt                                          <  1:1851339/93‑1 (MQ=255)

                                                aCGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGtt                                         <  1:4336122/93‑1 (MQ=255)

                                                aCGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGtt                                         <  1:1123006/93‑1 (MQ=255)

                                                aCGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGtt                                         <  2:4363923/93‑1 (MQ=255)

                                                aCGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGtt                                         <  1:815759/93‑1 (MQ=255)

                                                 cGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTg                                        <  1:749271/93‑1 (MQ=255)

                                                 cGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTg                                        >  2:1903418/1‑93 (MQ=255)

                                                  gcagcaTCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTgg                                       >  2:4051398/1‑93 (MQ=255)

                                                   cagcaTCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGgc                                      <  1:3687387/93‑1 (MQ=255)

                                                   cagcaTCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGgc                                      <  1:1527657/93‑1 (MQ=255)

                                                     gcaTCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGgcag                                    >  2:3177150/1‑93 (MQ=255)

                                                      caTCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGgcagc                                   <  2:1241361/93‑1 (MQ=255)

                                                       aTCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGgcagca                                  >  2:3279380/1‑93 (MQ=255)

                                                       aTCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGgcagca                                  >  1:2499579/1‑93 (MQ=255)

                                                        tCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGcagcag                                 >  2:3020110/1‑93 (MQ=255)

                                                         cccACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCagcagc                                <  1:3820578/93‑1 (MQ=255)

                                                          ccACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAgcg                               >  1:886215/1‑93 (MQ=255)

                                                           cACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAgcgc                              <  2:4401251/93‑1 (MQ=255)

                                                           cACCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAgcgc                              <  2:2234686/93‑1 (MQ=255)

                                                            aCCGCCGCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGcc                             <  1:1708910/93‑1 (MQ=255)

                                                             ccgccgCGCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCt                            >  2:3571230/1‑93 (MQ=255)

                                                                 cgcgCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAAt                        <  1:2537657/93‑1 (MQ=255)

                                                                 cgcgCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAAt                        <  1:604283/93‑1 (MQ=255)

                                                                 cgcgCCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAAt                        <  2:2484278/93‑1 (MQ=255)

                                                                   cgcCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCg                      <  2:116591/93‑1 (MQ=255)

                                                                    gcCAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCgg                     >  1:3203537/1‑93 (MQ=255)

                                                                    gcCAAAGCCACTAAAGAACAGCACATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCgg                     >  2:3651913/1‑93 (MQ=255)

                                                                     ccAAAGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGc                    >  2:2161118/1‑93 (MQ=255)

                                                                       aaaGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCa                  >  2:1484058/1‑93 (MQ=255)

                                                                       aaaGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCa                  <  2:2455933/93‑1 (MQ=255)

                                                                        aaGCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCAt                 >  1:4473267/1‑93 (MQ=255)

                                                                          gCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGa               <  2:4235866/93‑1 (MQ=255)

                                                                          gCCACTAAAGAACAGCGCATTGGCTACAGTATGCGGGCTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGa               >  2:435355/1‑93 (MQ=255)

                                                                           ccACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAg              >  1:1072497/1‑93 (MQ=255)

                                                                           ccACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAg              <  2:2700939/93‑1 (MQ=255)

                                                                            cACTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGt             <  2:331381/93‑1 (MQ=255)

                                                                             aCTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTa            >  1:2186226/1‑93 (MQ=255)

                                                                              cTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTGAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTag           >  2:1625393/1‑93 (MQ=255)

                                                                              cTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTag           >  2:860710/1‑93 (MQ=255)

                                                                              cTAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAACCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTag           <  1:999753/93‑1 (MQ=255)

                                                                               tAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTaga          >  2:2192043/1‑93 (MQ=255)

                                                                               tAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTaga          <  1:1178721/93‑1 (MQ=255)

                                                                               tAAAGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTaga          >  1:2274189/1‑93 (MQ=255)

                                                                                aaaGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTagag         <  1:2788486/93‑1 (MQ=255)

                                                                                 aaGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAgg        >  1:1401742/1‑93 (MQ=255)

                                                                                  aGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGa       <  1:2324368/93‑1 (MQ=255)

                                                                                  aGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGa       >  2:3841621/1‑93 (MQ=255)

                                                                                  aGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGa       <  2:582060/93‑1 (MQ=255)

                                                                                  aGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGa       >  2:707989/1‑93 (MQ=255)

                                                                                  aGAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGa       >  2:1438660/1‑93 (MQ=255)

                                                                                   gAACAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGaa      <  2:902504/93‑1 (MQ=255)

                                                                                    aaCAGCGCATTGGCTACAGTATGCGTGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGaaa     <  1:439277/93‑1 (MQ=255)

                                                                                    aaCAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGaaa     >  2:3958559/1‑93 (MQ=255)

                                                                                    aaCAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGaaa     >  1:4263631/1‑93 (MQ=255)

                                                                                     aCAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGAAAg    >  1:1418427/1‑93 (MQ=255)

                                                                                      cAGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGAAAGc   <  1:1389354/93‑1 (MQ=255)

                                                                                       aGCGCATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGAAAGCa  <  1:3097094/93‑1 (MQ=255)

CAAACGTACGCTTTGCGGGTTCCCAGCCAGTCACCGAGGAAGCCACCTACGCAGCATCCCACCGCCGCGCCAAAGCCACTAAAGAACAGCACATTGGCTACAGTATGCGGGTTATAAGCCAGATCGGTTTTCAGATACGTTGGCAGCAGCGCCTGAATCGGCCATGAGTAGAGGAAAGCA  >  NC_000913/3369173‑3369352

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