BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,211,779	T→C	R364R (CGT→CGC)	rpoD →	RNA polymerase sigma(70)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,211,779	0	T	C	97.3%	686.1 / 2.3	219	R364R (CGT→CGC)	rpoD	RNA polymerase sigma(70)
Reads supporting (aligned to +/- strand): ref base T (2/4); new base C (96/117); total (98/121)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

TGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAAT  >  NC_000913/3211689‑3211867

tGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCAt                                                                                        >  2:3958025/1‑93 (MQ=255)

  aagaagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGt                                                                                      >  1:2886972/1‑93 (MQ=255)

  aagaagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGt                                                                                      >  2:2660210/1‑93 (MQ=255)

   agaagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAATAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTc                                                                                     <  2:3413446/93‑1 (MQ=255)

   agaagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAAACGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTc                                                                                     >  2:1497173/1‑93 (MQ=255)

    gaagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTcc                                                                                    >  2:3816550/1‑93 (MQ=255)

     aagTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCa                                                                                   <  2:2051185/93‑1 (MQ=255)

      agTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCAt                                                                                  >  2:2590484/1‑93 (MQ=255)

       gTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATc                                                                                 <  1:713154/93‑1 (MQ=255)

       gTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATc                                                                                 >  2:1065127/1‑93 (MQ=255)

       gTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATc                                                                                 >  2:1674035/1‑93 (MQ=255)

       gTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATc                                                                                 <  1:3579650/93‑1 (MQ=255)

        tGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCg                                                                                <  2:2743927/93‑1 (MQ=255)

        tGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCg                                                                                >  2:859899/1‑93 (MQ=255)

         gCATCGCGCCCTGCAAAAACTGCAGCAGATAGTAGAAGAAACCGGCCATACCATCGAGCAGGATAAAGATATCAATCGTCGTATGTCCATCgt                                                                               >  2:4331378/1‑92 (MQ=255)

          cATCGCGCCTTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGt                                                                              <  1:422306/93‑1 (MQ=255)

          cATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGt                                                                              <  2:3320659/93‑1 (MQ=255)

            tCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGa                                                                            >  1:1757509/1‑93 (MQ=255)

            tCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGa                                                                            >  1:129726/1‑93 (MQ=255)

            tCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGa                                                                            <  2:1050465/93‑1 (MQ=255)

             cgcgCCNTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGaa                                                                           <  1:2587273/93‑1 (MQ=255)

                gcCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCg                                                                        >  1:3262173/1‑93 (MQ=255)

                    tGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAg                                                                    <  2:1820121/93‑1 (MQ=255)

                     gCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGc                                                                   <  1:612011/93‑1 (MQ=255)

                       aaaaaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGccc                                                                 <  1:1614569/93‑1 (MQ=255)

                       aaaaaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGccc                                                                 <  1:3424199/93‑1 (MQ=255)

                        aaaaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGcccc                                                                >  2:3509859/1‑92 (MQ=255)

                         aaaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgc                                                               <  1:1921108/93‑1 (MQ=255)

                          aaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgcg                                                              >  1:3986290/1‑92 (MQ=255)

                          aaCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgcc                                                              >  1:3624310/1‑93 (MQ=255)

                           aCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgccg                                                             >  2:757486/1‑93 (MQ=255)

                           aCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgccg                                                             >  1:493290/1‑93 (MQ=255)

                           aCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgccg                                                             <  1:376161/93‑1 (MQ=255)

                           aCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgccg                                                             >  2:1914995/1‑93 (MQ=255)

                           aCTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCccgccg                                                             >  2:2267245/1‑93 (MQ=255)

                            cTGCAGCAGATTGAAGAAGAAACCGTCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGt                                                            <  2:2398218/93‑1 (MQ=255)

                            cTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGt                                                            >  1:4619007/1‑93 (MQ=255)

                                 gcagATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCgaa                                                       <  1:529681/93‑1 (MQ=255)

                                 gcagATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGGTAAAGATATCAACCGTCGCATGGCCATCGGGGAAGCTAAAGCCCACTGTGCgag                                                       >  1:2684702/1‑92 (MQ=255)

                                    gATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCgaagaa                                                    >  1:2397496/1‑93 (MQ=255)

                                    gATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCgaagaa                                                    <  1:2004900/93‑1 (MQ=255)

                                     aTTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGaagaaa                                                   <  2:4047247/93‑1 (MQ=255)

                                     aTTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGaagaaa                                                   <  2:1615077/93‑1 (MQ=255)

                                       tGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGaagaaaga                                                 <  1:636873/93‑1 (MQ=255)

                                        gaagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAagag                                                >  1:1495783/1‑93 (MQ=255)

                                        gaagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAagag                                                >  2:4436602/1‑93 (MQ=255)

                                         aagaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAgaga                                               <  2:3800369/93‑1 (MQ=255)

                                          agaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGTCCGCCGTGCGAAGAAAGAGAt                                              <  2:4671581/93‑1 (MQ=255)

                                          agaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGAt                                              <  2:3043515/93‑1 (MQ=255)

                                          agaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGAt                                              <  1:821249/93‑1 (MQ=255)

                                          agaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGAt                                              <  2:1630653/93‑1 (MQ=255)

                                           gaagaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATg                                             >  2:284517/1‑93 (MQ=255)

                                              gaaCCCGGCCTGACCATCGAGCAGTTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGtt                                          <  2:1986126/93‑1 (MQ=255)

                                              gaaACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGtt                                          >  2:132259/1‑93 (MQ=255)

                                                 aCCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGaa                                       <  1:2923221/93‑1 (MQ=255)

                                                 aCCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGaa                                       <  1:1200588/93‑1 (MQ=255)

                                                  ccGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAg                                      >  2:1199167/1‑93 (MQ=255)

                                                     gCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGa                                   <  2:3690500/93‑1 (MQ=255)

                                                      ccTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGGTGGTTGAAGCGaa                                  >  2:4526431/1‑93 (MQ=255)

                                                      ccTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGaa                                  >  1:429936/1‑93 (MQ=255)

                                                       cTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTTCGAAGAAAGAGATGGTTGAAGCGAAc                                 >  1:1593439/1‑93 (MQ=255)

                                                       cTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGaat                                 >  2:804924/1‑92 (MQ=255)

                                                       cTGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAAc                                 <  2:603814/93‑1 (MQ=255)

                                                        tGACCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACt                                >  1:188213/1‑93 (MQ=255)

                                                        tGACCATCGAGCAGGTTAAAGATATCAACCGCCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACt                                >  2:2099332/1‑93 (MQ=255)

                                                          aCCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTa                              <  1:1314997/93‑1 (MQ=255)

                                                          aCCATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTa                              <  1:407620/93‑1 (MQ=255)

                                                            cATCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACg                            >  1:3903097/1‑93 (MQ=255)

                                                             aTCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGt                           >  2:844548/1‑93 (MQ=255)

                                                             aTCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGt                           <  2:2526375/93‑1 (MQ=255)

                                                              tCGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTc                          <  1:2819283/93‑1 (MQ=255)

                                                               cGAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCt                         <  2:3101160/93‑1 (MQ=255)

                                                                gAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTg                        >  2:3775080/1‑93 (MQ=255)

                                                                gAGCAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTg                        >  1:1238612/1‑93 (MQ=255)

                                                                  gCAGGTTAAAGATATCAACCGTCGTATGTCGATCGGTAAAGCAAAAGCCCGCGGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGt                      <  1:4014288/93‑1 (MQ=255)

                                                                   cAGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGtt                     <  1:1802254/93‑1 (MQ=255)

                                                                    aGGTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGttt                    >  1:1030323/1‑92 (MQ=255)

                                                                     ggTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTAt                   <  2:3005236/93‑1 (MQ=255)

                                                                     ggTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTAt                   <  1:1794021/93‑1 (MQ=255)

                                                                     ggTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTAt                   >  2:565074/1‑93 (MQ=255)

                                                                     ggTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTAt                   >  1:4369772/1‑93 (MQ=255)

                                                                     ggTTAAAGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTAt                   >  1:4492667/1‑93 (MQ=255)

                                                                         aaaGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCt               <  1:2566362/93‑1 (MQ=255)

                                                                         aaaGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCt               <  1:3070931/93‑1 (MQ=255)

                                                                          aaGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTa              <  1:2070055/93‑1 (MQ=255)

                                                                          aaGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTa              <  2:3729970/93‑1 (MQ=255)

                                                                           aGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTAt             >  1:1113159/1‑93 (MQ=255)

                                                                           aGATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTAt             <  2:3216622/93‑1 (MQ=255)

                                                                            gATATCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATc            <  2:1375948/93‑1 (MQ=255)

                                                                             atatCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCg           <  2:1796786/93‑1 (MQ=255)

                                                                              tatCAACCTTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGc          <  2:2358656/93‑1 (MQ=255)

                                                                              tatCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCga          >  2:1121801/1‑92 (MQ=255)

                                                                              tatCAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGc          >  2:2919744/1‑93 (MQ=255)

                                                                                 cAACCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTaa       >  2:1889392/1‑93 (MQ=255)

                                                                                  aaCCGTCGCATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAg      >  1:3381941/1‑93 (MQ=255)

                                                                                     cgtcgcATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGaaa   <  1:1137801/93‑1 (MQ=255)

                                                                                      gtcgcATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAAt  <  2:2701957/93‑1 (MQ=255)

                                                                                      gtcgcATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAAt  >  2:2440749/1‑93 (MQ=255)

                                                                                      gtcgcATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAAt  >  1:3510454/1‑93 (MQ=255)

                                                                                      gtcgcATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAAt  >  1:219622/1‑93 (MQ=255)

TGAAGAAGTGCATCGCGCCCTGCAAAAACTGCAGCAGATTGAAGAAGAAACCGGCCTGACCATCGAGCAGGTTAAAGATATCAACCGTCGTATGTCCATCGGTGAAGCGAAAGCCCGCCGTGCGAAGAAAGAGATGGTTGAAGCGAACTTACGTCTGGTTATTTCTATCGCTAAGAAAT  >  NC_000913/3211689‑3211867

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