BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,415	A→G	T230A (ACG→GCG)	b3614 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,415	0	A	G	99.6%	874.9 / ‑5.2	270	T230A (ACG→GCG)	b3614	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (144/125); total (145/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 280 total aligned reads displayed.

CGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAT  >  NC_000913/3786325‑3786504

cGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCg                                                                                         <  1:2149751/93‑1 (MQ=255)

 gCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACTCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTgcgc                                                                                        >  1:6041805/1‑93 (MQ=255)

 gCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTgcgc                                                                                        <  2:2717250/93‑1 (MQ=255)

   aaCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTg                                                                                      >  2:5681553/1‑93 (MQ=255)

   aaCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTg                                                                                      >  2:4446563/1‑93 (MQ=255)

    aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGCCATTGCGCTgg                                                                                     >  2:3548438/1‑93 (MQ=255)

    aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                     >  2:745727/1‑93 (MQ=255)

    aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                     >  1:1243702/1‑93 (MQ=255)

    aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                     <  2:5992689/93‑1 (MQ=255)

     cGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGt                                                                                    <  2:1752577/93‑1 (MQ=255)

     cGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGt                                                                                    <  1:3437627/93‑1 (MQ=255)

      ggTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTg                                                                                   >  1:4322441/1‑93 (MQ=255)

      ggTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTg                                                                                   <  1:3972528/93‑1 (MQ=255)

      ggTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTg                                                                                   >  2:5061478/1‑93 (MQ=255)

          tggacattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAg                                                                               <  1:3390967/89‑1 (MQ=255)

           gACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGc                                                                              >  2:1023567/1‑93 (MQ=255)

           gACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGc                                                                              >  2:1058670/1‑93 (MQ=255)

            aCCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGcc                                                                             <  2:48337/93‑1 (MQ=255)

            aCCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGcc                                                                             >  2:5047176/1‑93 (MQ=255)

             ccattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCca                                                                            <  1:3919683/93‑1 (MQ=255)

              cattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcag                                                                           <  1:879831/93‑1 (MQ=255)

               attgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagc                                                                          >  1:5505437/1‑93 (MQ=255)

               attgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagc                                                                          <  2:4199038/93‑1 (MQ=255)

               attgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagc                                                                          >  2:5323102/1‑93 (MQ=255)

                ttgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagca                                                                         >  1:2371438/1‑93 (MQ=255)

                  gccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCagcagc                                                                       <  2:2933218/93‑1 (MQ=255)

                   ccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCt                                                                      <  1:2744583/93‑1 (MQ=255)

                    cattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCtt                                                                     >  2:1394960/1‑93 (MQ=255)

                    cattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCtt                                                                     >  2:5942630/1‑93 (MQ=255)

                     attgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTg                                                                    >  2:3148400/1‑93 (MQ=255)

                      ttgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGc                                                                   <  1:5147125/93‑1 (MQ=255)

                          gCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaa                                                               <  1:4180363/93‑1 (MQ=255)

                           cATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaat                                                              <  2:3667539/93‑1 (MQ=255)

                            aTCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                             <  2:2426164/93‑1 (MQ=255)

                             tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                            <  1:2015797/93‑1 (MQ=255)

                             tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                            >  1:6146455/1‑93 (MQ=255)

                               cTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAc                                                          <  2:874318/93‑1 (MQ=255)

                                tCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAcc                                                         >  1:1107768/1‑93 (MQ=255)

                                tCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAcc                                                         >  2:1489429/1‑93 (MQ=255)

                                 cACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCg                                                        >  1:5211545/1‑93 (MQ=255)

                                  aCCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGc                                                       <  2:3484975/93‑1 (MQ=255)

                                   cccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCa                                                      <  1:4275389/93‑1 (MQ=255)

                                    ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     >  1:4708577/1‑93 (MQ=255)

                                    ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     <  2:3235862/93‑1 (MQ=255)

                                      tGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGt                                                   >  1:3597000/1‑93 (MQ=255)

                                       gCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGtt                                                  <  1:4862465/93‑1 (MQ=255)

                                         aaCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGa                                                >  2:5432263/1‑93 (MQ=255)

                                           cGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATa                                              >  1:978918/1‑93 (MQ=255)

                                             gTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACt                                            <  1:3477728/93‑1 (MQ=255)

                                              tACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACtt                                           >  2:2166805/1‑93 (MQ=255)

                                              tACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACtt                                           <  1:5858628/93‑1 (MQ=255)

                                               aCGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTc                                          >  1:3253016/1‑93 (MQ=255)

                                                 gcgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTa                                        >  1:3025477/1‑93 (MQ=255)

                                                 gcgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTa                                        <  1:5232786/93‑1 (MQ=255)

                                                  cgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTac                                       <  1:1310094/93‑1 (MQ=255)

                                                   gtgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTaca                                      >  1:3872730/1‑93 (MQ=255)

                                                   gtgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTaca                                      >  1:3227517/1‑93 (MQ=255)

                                                    tgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTacac                                     <  2:5593550/93‑1 (MQ=255)

                                                    tgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTacac                                     >  2:3042103/1‑93 (MQ=255)

                                                    tgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTacac                                     <  1:761689/93‑1 (MQ=255)

                                                     gtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAcc                                    <  1:3172479/93‑1 (MQ=255)

                                                      ttGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccg                                   >  2:6039580/1‑93 (MQ=255)

                                                       tGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgc                                  <  2:3346310/93‑1 (MQ=255)

                                                        gCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgcc                                 <  2:456482/93‑1 (MQ=255)

                                                         cAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgccg                                >  1:3817182/1‑93 (MQ=255)

                                                          aaCAGATGGTTTATAACCTGCCGCCAGACTTTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGa                               >  2:3806550/1‑93 (MQ=255)

                                                            cAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaa                             >  1:1762686/1‑93 (MQ=255)

                                                             aGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaa                            >  2:2315058/1‑93 (MQ=255)

                                                             aGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaa                            >  1:1977499/1‑93 (MQ=255)

                                                              gATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaaa                           <  1:2890984/93‑1 (MQ=255)

                                                              gATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaaa                           >  1:5942983/1‑93 (MQ=255)

                                                              gATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaaa                           >  2:5535357/1‑93 (MQ=255)

                                                              gATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaaa                           <  2:1659189/93‑1 (MQ=255)

                                                               aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          <  2:2052908/93‑1 (MQ=255)

                                                                 ggTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGc                        <  2:1279345/93‑1 (MQ=255)

                                                                 ggTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGc                        <  1:2580328/93‑1 (MQ=255)

                                                                    ttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                     >  1:439278/1‑93 (MQ=255)

                                                                    ttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                     <  2:2661850/93‑1 (MQ=255)

                                                                    ttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                     <  2:4645995/93‑1 (MQ=255)

                                                                       taACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGcc                  <  1:5683206/93‑1 (MQ=255)

                                                                        aaCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCt                 <  2:2987279/93‑1 (MQ=255)

                                                                         aCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTg                <  1:589482/93‑1 (MQ=255)

                                                                          ccTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGa               >  1:5760576/1‑93 (MQ=255)

                                                                          ccTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGa               <  1:1069156/93‑1 (MQ=255)

                                                                          ccTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGa               <  1:3519215/93‑1 (MQ=255)

                                                                           cTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAACAACGCTGTGCCTGAt              >  2:4930185/1‑93 (MQ=255)

                                                                           cTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGAt              >  1:1812573/1‑93 (MQ=255)

                                                                              ccgccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCa           <  2:4143964/93‑1 (MQ=255)

                                                                               cgccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAc          <  2:3932930/93‑1 (MQ=255)

                                                                               cgccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAc          >  2:211925/1‑93 (MQ=255)

                                                                                gccAGTCATTGCGCTGGGGAAAGCCAGTAGCTTGCTGAATGAACCGCAGGCTGATACGTCTACACTGCCGAAAAACGCTGTGCCTGACGCAcc         >  2:6304284/1‑93 (MQ=255)

                                                                                gccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAcc         >  1:2311857/1‑93 (MQ=255)

                                                                                 ccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATCAACCGCAGGTTGAGACTTCTACACCGCCGAAAAAAGCTGTGCCTGATGCACcg        >  1:374841/1‑93 (MQ=255)

                                                                                  cAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgc       <  1:3671925/93‑1 (MQ=255)

                                                                                   aGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgcg      >  2:1925536/1‑93 (MQ=255)

                                                                                    gACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGt     >  2:5079645/1‑93 (MQ=255)

                                                                                    gACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGt     >  2:5179323/1‑93 (MQ=255)

                                                                                    gACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGt     >  1:6235284/1‑93 (MQ=255)

                                                                                     aCATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTa    >  2:5835969/1‑93 (MQ=255)

                                                                                     aCATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTa    >  2:2872625/1‑93 (MQ=255)

CGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAT  >  NC_000913/3786325‑3786504

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