BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	934,731	A→G	Q762R (CAG→CGG)	ftsK →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	934,731	0	A	G	99.4%	512.7 / ‑4.5	161	Q762R (CAG→CGG)	ftsK	cell division protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (100/59); total (102/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 222 total aligned reads displayed.

TTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATC  >  NC_000913/934641‑934820

ttACGCCAATTGTTGACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggc                                                                                         <  1:3552375/93‑1 (MQ=255)

ttACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggc                                                                                         <  2:1454934/93‑1 (MQ=255)

 tACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggcc                                                                                        <  1:2507289/93‑1 (MQ=255)

  aCGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggccg                                                                                       >  1:4217705/1‑93 (MQ=255)

   cGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggccgc                                                                                      <  2:6208390/93‑1 (MQ=255)

   cGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCgccgcggccgc                                                                                      <  2:7546859/93‑1 (MQ=255)

       aaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTa                                                                                  >  1:3971575/1‑93 (MQ=255)

       aaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTa                                                                                  <  2:1934465/93‑1 (MQ=255)

       aaTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTa                                                                                  <  1:2803764/93‑1 (MQ=255)

        aTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTAt                                                                                 <  1:5160309/93‑1 (MQ=255)

        aTTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTAt                                                                                 <  1:8169755/93‑1 (MQ=255)

         ttgttgACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATc                                                                                <  2:4417673/93‑1 (MQ=255)

         ttgttgACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATc                                                                                <  2:3152035/93‑1 (MQ=255)

         ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATc                                                                                >  1:8688954/1‑93 (MQ=255)

         ttgttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATc                                                                                >  2:2725999/1‑93 (MQ=255)

           tgtgACCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcag                                                                              <  2:4058379/91‑1 (MQ=255)

           gttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcag                                                                              <  2:4214792/93‑1 (MQ=255)

           gttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcag                                                                              <  1:1495738/93‑1 (MQ=255)

           gttgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcag                                                                              >  1:8322482/1‑93 (MQ=255)

            ttgAGCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcagc                                                                             <  1:3505021/93‑1 (MQ=255)

              gAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATcagcac                                                                           >  1:7335618/1‑92 (MQ=255)

                aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGcc                                                                         >  1:3255354/1‑93 (MQ=255)

                aCCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGcc                                                                         >  1:5609338/1‑93 (MQ=255)

                aCCTGTACAGCAGCCGCAACAACCGGTAGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGGATCaccagct                                                                         >  2:1149150/1‑92 (MQ=18)

                   tGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCcaca                                                                      >  1:7759945/1‑93 (MQ=255)

                    gTACAGCAGCCGCAACACCCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCacaa                                                                     <  1:2358326/93‑1 (MQ=255)

                    gTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCacaa                                                                     <  1:1124209/93‑1 (MQ=255)

                      aaggcaccCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCacaaca                                                                   <  1:7931741/90‑1 (MQ=17)

                       cagcagCCGCAACAACCGGTTGCACCGCAGCAGCACTATCAGCAGCCGCAACAACCAGTTCCACCGCAGCCGCAGTATCAGCAGCCACAACAg                                                                  >  2:1496256/1‑93 (MQ=255)

                        agcggcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGc                                                                 <  1:4264993/93‑1 (MQ=255)

                        agcagcCGCAACACCCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGc                                                                 <  1:8113116/93‑1 (MQ=255)

                        agcaccCGAAACACCCGGTTCCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGc                                                                 <  2:7715124/93‑1 (MQ=17)

                         gcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGGATCAGCAGCTACAACAGcc                                                                >  2:6331046/1‑93 (MQ=255)

                             ccGCAACAACGGGTGGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGtt                                                            <  2:5451084/93‑1 (MQ=255)

                               gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                                                          >  2:2277112/1‑93 (MQ=255)

                               gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                                                          >  2:4085219/1‑93 (MQ=255)

                               gCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgc                                                          >  1:509683/1‑93 (MQ=255)

                                 aacaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgcgc                                                        >  2:6176254/1‑93 (MQ=255)

                                 aacaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTgcgc                                                        >  1:1203661/1‑93 (MQ=255)

                                  acaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgcc                                                       >  2:6505176/1‑93 (MQ=255)

                                  acaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgcc                                                       >  2:6809119/1‑93 (MQ=255)

                                   aaacCGGTTGCACCGCAGCAGCAATACCACCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccg                                                      <  1:1330825/92‑1 (MQ=17)

                                    aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccgc                                                     >  2:8680181/1‑93 (MQ=255)

                                    aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgccgc                                                     <  2:6465704/93‑1 (MQ=255)

                                    aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGagcagc                                                     >  1:4963848/1‑93 (MQ=255)

                                          ttGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCcac                                               >  1:6462181/1‑93 (MQ=255)

                                            gCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa                                             >  2:2120776/1‑93 (MQ=255)

                                            cccccGCAGCACCAATATCAGCAGCCGCAACAACCAGTCCCGCCGCGCCCGCATTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACaa                                             <  2:8278967/90‑1 (MQ=18)

                                             ctccGCAGCAGCAATACCAGCAGCCGCACCAGCCAGTTCCGCGTCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAat                                            <  2:5244009/91‑1 (MQ=18)

                                             cACCGCAGCAGAAATACCACCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAat                                            <  2:5011057/93‑1 (MQ=21)

                                               ccGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAatat                                          >  1:3921738/1‑93 (MQ=255)

                                                cGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATc                                         >  2:1202691/1‑93 (MQ=255)

                                                cGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATc                                         >  2:2597915/1‑93 (MQ=255)

                                                cGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCCGTTGCGCCGCAGCCACAATATc                                         >  2:1554655/1‑93 (MQ=255)

                                                 gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATca                                        >  1:1767760/1‑93 (MQ=255)

                                                  cagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcag                                       <  2:1872824/93‑1 (MQ=255)

                                                   agcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagc                                      <  2:475268/93‑1 (MQ=255)

                                                   agcagcaACATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATcagc                                      <  1:654423/93‑1 (MQ=255)

                                                      agcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCagcagc                                   >  1:7231882/1‑93 (MQ=255)

                                                       gcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGcc                                  <  2:727700/93‑1 (MQ=255)

                                                         aaTATCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGc                                >  1:8068432/1‑93 (MQ=255)

                                                         aaTATCAGCAGCCGCAACAACCAGTTCCGCCGCGCCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGc                                <  1:3056372/93‑1 (MQ=255)

                                                          atatCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGca                               >  2:39090/1‑93 (MQ=255)

                                                           tatCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaa                              >  2:6041405/1‑93 (MQ=255)

                                                           tacaagcagcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaa                              <  1:4765137/89‑1 (MQ=255)

                                                           aaccagcggCCGCAACACCCAGTCCCGCCGCGGCCGCAGTACCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaa                              <  2:3746660/90‑1 (MQ=21)

                                                            atCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaac                             >  1:8769349/1‑93 (MQ=255)

                                                             tCAGCAGCCGCAACAACCAGTTCCGCCGCGGCCTCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaaca                            <  1:8519443/93‑1 (MQ=255)

                                                               gccggcCGAACCAACCAGTCCCGCCGGGCCTGAGGTCTCCGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCaacaac                          <  2:60768/89‑1 (MQ=12)

                                                               cgcagcCCCACCACCCAGTCCCGCCGCGGCCGCAGTATCACCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCaacaac                          <  2:3156161/92‑1 (MQ=21)

                                                               agcagcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCaacaac                          <  1:5447568/93‑1 (MQ=255)

                                                               agcagcCGCAACAACCAGTCCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCaacaac                          <  2:6101963/93‑1 (MQ=255)

                                                                  agcCGCAACACCCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCgg                       <  1:8389790/93‑1 (MQ=255)

                                                                   gcCGCAACACCCAGTCCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                      <  2:1860031/93‑1 (MQ=255)

                                                                   gcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCGGCAACAACCGGt                      >  1:2223309/1‑93 (MQ=255)

                                                                   gcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                      >  2:7248319/1‑93 (MQ=255)

                                                                   gcCGCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                      >  2:8232053/1‑93 (MQ=255)

                                                                      gCAACAACCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgc                   >  1:7029638/1‑93 (MQ=255)

                                                                       caacaaCCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTggc                  >  1:6836065/1‑91 (MQ=255)

                                                                       caacaaCCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgcg                  >  1:5381598/1‑93 (MQ=255)

                                                                       caacaaCCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgcg                  >  2:1089310/1‑93 (MQ=255)

                                                                        aacaacCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTgcgc                 >  1:1812598/1‑93 (MQ=255)

                                                                         acaccCAGTCCCCCGGCGGCAGCAGTATCACCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGcc                <  2:5613141/93‑1 (MQ=18)

                                                                         acaacCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGcc                >  2:1959761/1‑93 (MQ=255)

                                                                         acaacCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGcc                >  1:4479675/1‑93 (MQ=255)

                                                                          caacCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCca               >  1:7549844/1‑93 (MQ=255)

                                                                           accgaGTCCCGCCGCGGCCAAAGAATCAGCAGCAACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcac              <  2:7955511/89‑1 (MQ=17)

                                                                           accaaGTCCCGCCGCGGCCGCATTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcac              <  1:8012998/89‑1 (MQ=18)

                                                                           aacCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcac              <  1:7451080/93‑1 (MQ=255)

                                                                           aacCAGTTCCCCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcac              <  1:6279606/93‑1 (MQ=255)

                                                                            acCAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca             >  1:7179166/1‑93 (MQ=255)

                                                                             ccAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGCTGCGCCGCAGCCACAATATCAGCAGCCGCAACAGCCGGTTGCGCCAcag            >  2:385524/1‑93 (MQ=255)

                                                                              cAGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagc           >  1:6209464/1‑93 (MQ=255)

                                                                               cgccccgccgCGGCCGCAATAGCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagca          <  1:2109989/89‑1 (MQ=17)

                                                                               aGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagca          <  1:2591126/93‑1 (MQ=255)

                                                                               aGTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagca          >  2:8566838/1‑93 (MQ=255)

                                                                                gTTCCGCCGCGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagcag         >  2:4134789/1‑93 (MQ=255)

                                                                                   ccgccgGGGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCaa      <  2:2433976/93‑1 (MQ=255)

                                                                                     gccgcggccgcaGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAata    <  2:4527644/93‑1 (MQ=255)

                                                                                       cgcggccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATc  >  2:5537380/1‑93 (MQ=255)

TTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATC  >  NC_000913/934641‑934820

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