BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,262,935	A→G	intergenic (‑212/‑2)	hemM ← / → hemA	HemM/HemA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,262,935	0	A	G	88.7%	563.3 / 56.0	231	intergenic (‑212/‑2)	hemM/hemA	HemM/HemA
Reads supporting (aligned to +/- strand): ref base A (15/11); new base G (111/94); total (126/105)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 238 total aligned reads displayed.

TCCAAAAGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAG  >  NC_000913/1262844‑1263026

   aaaaGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATg                                                                                         >  2:4396699/1‑93 (MQ=255)

   aaaaGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATg                                                                                         >  2:2878325/1‑93 (MQ=255)

   aaaaGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATg                                                                                         <  1:5602690/93‑1 (MQ=255)

   aaaaGGGGGCGCTCTCTTTTATTAATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATg                                                                                         <  2:5526911/93‑1 (MQ=255)

    aaaGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGa                                                                                        <  2:2866438/93‑1 (MQ=255)

      aGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGAcc                                                                                      >  2:2493717/1‑93 (MQ=255)

      aGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGAcc                                                                                      >  1:6059450/1‑93 (MQ=255)

      aGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGAcc                                                                                      >  1:5770023/1‑93 (MQ=255)

      aGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGAcc                                                                                      <  1:3836256/93‑1 (MQ=255)

       gggggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGAccc                                                                                     <  1:4104675/93‑1 (MQ=255)

         gggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCtt                                                                                   >  2:9002113/1‑93 (MQ=255)

          ggCGCTCTCTTTTATTGATCTTACGCCTCCTGTATGATCCAATCAGGCTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCttt                                                                                  <  1:2188818/93‑1 (MQ=255)

          ggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCttt                                                                                  >  1:1686978/1‑93 (MQ=255)

          ggCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCCCAGGCATGACCCttt                                                                                  >  1:6323551/1‑93 (MQ=255)

           gcgcTCCCTTTTATTGATCTTACGCATTCTGTATGATGCACGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACACGACCCtttt                                                                                 <  2:5968208/93‑1 (MQ=255)

             gcTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAg                                                                               >  2:699743/1‑93 (MQ=255)

              ctctctTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGc                                                                              >  2:500287/1‑93 (MQ=255)

               tctctTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCa                                                                             >  1:9346016/1‑93 (MQ=255)

                ctctTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCAc                                                                            >  1:9247407/1‑93 (MQ=255)

                 tctTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACt                                                                           >  2:7555044/1‑93 (MQ=255)

                 tctTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACt                                                                           <  2:6401353/93‑1 (MQ=255)

                 tctTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACt                                                                           >  1:7499897/1‑93 (MQ=255)

                  ctTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTc                                                                          <  2:2182393/93‑1 (MQ=255)

                   ttttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGCCCCTTTTAGCACTCg                                                                         <  1:6009357/93‑1 (MQ=255)

                   ttttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCg                                                                         >  1:3994581/1‑93 (MQ=255)

                    tttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCgg                                                                        >  1:714453/1‑93 (MQ=255)

                     ttATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGt                                                                       >  1:2349424/1‑93 (MQ=255)

                       aTTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTAt                                                                     >  1:2909524/1‑93 (MQ=255)

                        ttGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATc                                                                    >  1:2366571/1‑93 (MQ=255)

                        ttGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTCTAGCACTCGGTAta                                                                    >  1:4630332/1‑92 (MQ=255)

                         tGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCa                                                                   >  2:937846/1‑93 (MQ=255)

                          gATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCaa                                                                  <  1:4260373/93‑1 (MQ=255)

                           aTCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAAc                                                                 >  1:7143506/1‑93 (MQ=255)

                            tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAAcc                                                                >  1:6572017/1‑93 (MQ=255)

                            tCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAAcc                                                                >  2:6743186/1‑93 (MQ=255)

                              ttACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGCATCAACCAt                                                              >  2:419618/1‑93 (MQ=255)

                              ttACGCATCCTGTATGATGCAAGCAGACTAAACCTATCAACGATGGTCTTATTTCCCACAGACATGACCATTTTAGCACTCGGTATCAAcaca                                                              >  2:3223924/1‑90 (MQ=255)

                                 cGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATaaa                                                           <  2:8412916/93‑1 (MQ=255)

                                 cGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATaaa                                                           <  2:1774121/93‑1 (MQ=255)

                                 cGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATaaa                                                           <  2:9783494/93‑1 (MQ=255)

                                  gCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATaaaa                                                          <  2:6077539/93‑1 (MQ=255)

                                  gCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATaaaa                                                          >  1:1833751/1‑93 (MQ=255)

                                   cATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAAc                                                         >  1:2763708/1‑93 (MQ=255)

                                       cTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCa                                                     >  2:8158059/1‑93 (MQ=255)

                                        tGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCAc                                                    <  2:7141607/93‑1 (MQ=255)

                                        tGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCAc                                                    <  2:1411279/93‑1 (MQ=255)

                                        tGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCAc                                                    <  2:5940235/93‑1 (MQ=255)

                                         gTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCAcc                                                   >  1:7991514/1‑93 (MQ=255)

                                         gTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCAcc                                                   >  2:9272027/1‑93 (MQ=255)

                                         gTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCAcc                                                   <  2:3569853/93‑1 (MQ=255)

                                           atgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTTCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGTACCTg                                                 >  1:9898074/1‑93 (MQ=255)

                                           atgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTg                                                 <  1:5065018/93‑1 (MQ=255)

                                           atgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTg                                                 <  2:8516280/93‑1 (MQ=255)

                                            tgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGt                                                <  2:267419/93‑1 (MQ=255)

                                            tgatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGt                                                >  1:8747573/1‑93 (MQ=255)

                                             gatgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCTGTATCAACCATAAAACGGCACCTGTa                                               >  1:1022468/1‑93 (MQ=255)

                                              atgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTAt                                              <  2:9739067/93‑1 (MQ=255)

                                              atgCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTAt                                              >  2:8310225/1‑93 (MQ=255)

                                                 cAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGc                                           <  1:6111964/93‑1 (MQ=255)

                                                     cAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCg                                       <  2:1586881/93‑1 (MQ=255)

                                                     cAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCg                                       >  2:7560136/1‑93 (MQ=255)

                                                      agccTAACCCTATCATCGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCga                                      <  2:8923218/90‑1 (MQ=255)

                                                      aGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCga                                      >  1:1994338/1‑93 (MQ=255)

                                                      aGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCga                                      <  2:3094506/93‑1 (MQ=255)

                                                       gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgag                                     >  2:6690388/1‑93 (MQ=255)

                                                       gACTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgag                                     >  2:1788381/1‑93 (MQ=255)

                                                        aCTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCgaga                                    >  2:8696560/1‑93 (MQ=255)

                                                         cTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGaa                                   <  2:3780491/93‑1 (MQ=255)

                                                         cTAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGaa                                   <  2:4125483/93‑1 (MQ=255)

                                                          tAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAAc                                  >  2:8246335/1‑93 (MQ=255)

                                                          tAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAAc                                  >  1:2751293/1‑93 (MQ=255)

                                                          tAACCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAAc                                  >  2:579809/1‑93 (MQ=255)

                                                          tAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAAc                                  >  2:4817477/1‑93 (MQ=255)

                                                           aaCCCTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACg                                 <  2:3461411/93‑1 (MQ=255)

                                                              ccTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt                              <  2:9086927/93‑1 (MQ=255)

                                                              ccTATCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACgtgt                              >  2:6885910/1‑93 (MQ=255)

                                                                 aTCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATc                           <  2:8625167/93‑1 (MQ=255)

                                                                 aTCAACGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATc                           <  1:32418/93‑1 (MQ=255)

                                                                      cGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGtttt                      >  2:7864558/1‑93 (MQ=255)

                                                                      cGTTGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGtttt                      <  1:5334033/93‑1 (MQ=255)

                                                                        ttGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGATAACGTGTATCGTTTTCg                    >  1:8970110/1‑93 (MQ=255)

                                                                        ttGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCg                    >  2:5372338/1‑93 (MQ=255)

                                                                        ttGGTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCg                    >  1:9130442/1‑93 (MQ=255)

                                                                          ggTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGcc                  >  2:4779992/1‑93 (MQ=255)

                                                                          ggTATTATTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATAGCTGCGAGAACGTGTATCGTTTTCGcc                  >  2:7496788/1‑93 (MQ=255)

                                                                            tattatTTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCgg                <  2:2774131/93‑1 (MQ=255)

                                                                              ttattTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGAt              >  2:5598692/1‑93 (MQ=255)

                                                                              ttattTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGAt              <  1:3528219/93‑1 (MQ=255)

                                                                               tattTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATa             >  2:2235354/1‑93 (MQ=255)

                                                                                attTCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATaa            <  2:1509860/93‑1 (MQ=255)

                                                                                 tttCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAg           >  2:2102391/1‑93 (MQ=255)

                                                                                  ttCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGc          >  2:1071785/1‑93 (MQ=255)

                                                                                   tCCCGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCt         >  2:4539217/1‑93 (MQ=255)

                                                                                     ccGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTca       >  1:3629147/1‑92 (MQ=255)

                                                                                      cGCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCAGGATAAGCTcaa      <  1:8607491/93‑3 (MQ=255)

                                                                                       gCAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAAt     <  1:193691/93‑1 (MQ=255)

                                                                                        cAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc    <  2:6301041/93‑1 (MQ=255)

                                                                                        cAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc    <  2:5686766/93‑1 (MQ=255)

                                                                                        cAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc    >  1:6901090/1‑93 (MQ=255)

                                                                                        cAGGCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc    >  1:8490530/1‑93 (MQ=255)

TCCAAAAGGGGGCGCTCTCTTTTATTGATCTTACGCATCCTGTATGATGCAAGCAGACTAACCCTATCAACGTTGGTATTATTTCCCGCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAG  >  NC_000913/1262844‑1263026

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