BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	937,417	G→C	A16P (GCA→CCA)	lolA →	lipoprotein chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	937,417	0	G	C	98.6%	529.9 / ‑1.0	145	A16P (GCA→CCA)	lolA	lipoprotein chaperone
Reads supporting (aligned to +/- strand): ref base G (2/0); new base C (69/74); total (71/74)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 158 total aligned reads displayed.

TAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACG  >  NC_000913/937320‑937513

tAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTacca                                                                                                <  2:273316/100‑1 (MQ=255)

tAGAGAGGAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTacca                                                                                                <  2:2801959/100‑1 (MQ=255)

 agagagTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTtagtagatc                                                                                               >  1:1178853/1‑97 (MQ=255)

  gagagTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTaccaag                                                                                              >  2:3649860/1‑100 (MQ=255)

   agagTAGAGGGAACTCCCGATCGGTATTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTaccaagc                                                                                             <  1:947651/100‑1 (MQ=255)

      gTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAagcagc                                                                                          >  2:824252/1‑100 (MQ=255)

      gTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAagcagc                                                                                          >  1:1381799/1‑100 (MQ=255)

        agagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCAGCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGt                                                                                        >  1:3978250/1‑100 (MQ=255)

        agagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGt                                                                                        >  2:1907139/1‑100 (MQ=255)

        agacGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCAACTGTGCATTACTCTCAAGCTTATTACCAagcagc                                                                                          <  2:1814212/98‑1 (MQ=255)

         tacgtcaCTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGt                                                                                        >  1:2921399/7‑99 (MQ=255)

         gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGtt                                                                                       >  2:1247603/1‑100 (MQ=255)

         gagGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGtt                                                                                       >  1:3875418/1‑100 (MQ=255)

          agGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGttt                                                                                      <  2:437638/100‑1 (MQ=255)

          agGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACATGTGCATTACTCTCAAGCTTAGTACCACACAGCGttt                                                                                      >  2:3983767/1‑100 (MQ=255)

           gggAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTg                                                                                     <  1:386694/100‑1 (MQ=255)

           gggAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTg                                                                                     <  1:876501/100‑1 (MQ=255)

            ggAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTgg                                                                                    <  1:2871702/100‑1 (MQ=255)

             gAACTCCCGATCGGGAGTGACGTAATTTGAGGAACAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTggg                                                                                   <  2:2201758/100‑1 (MQ=255)

               aCTCCCGATCGGGAGTGACGTAATTTGCGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCt                                                                                 <  1:1475678/100‑1 (MQ=255)

                cTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTg                                                                                >  1:3195565/1‑100 (MQ=255)

                  cccGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGAt                                                                              >  1:2744174/1‑100 (MQ=255)

                  cccGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGAt                                                                              <  2:2145660/100‑1 (MQ=255)

                  cccGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGAt                                                                              <  2:4100469/100‑1 (MQ=255)

                   ccGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATg                                                                             >  1:2859893/1‑100 (MQ=255)

                    cGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGc                                                                            <  1:3992019/100‑1 (MQ=255)

                       tCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGc                                                                         <  2:2994207/100‑1 (MQ=255)

                        cGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCa                                                                        <  1:4100731/100‑1 (MQ=255)

                          ggAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAg                                                                      <  1:781590/100‑1 (MQ=255)

                             gTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGa                                                                   <  1:1658647/100‑1 (MQ=255)

                               gaagTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATc                                                                 <  2:2959490/97‑1 (MQ=255)

                               gACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATc                                                                 <  2:1972416/100‑1 (MQ=255)

                               gACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATc                                                                 >  1:2394161/1‑100 (MQ=255)

                                aCGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCt                                                                <  1:3781369/100‑1 (MQ=255)

                                aCGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCt                                                                <  1:1441603/100‑1 (MQ=255)

                                 cGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTg                                                               <  2:2331239/100‑1 (MQ=255)

                                 cGTAATTTGAGGAATAATGATGAAAAAAATCGCAATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGAAAGCGATCTg                                                               >  2:2754150/1‑100 (MQ=255)

                                   tAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGaa                                                             >  2:368608/1‑100 (MQ=255)

                                    aaTTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGaaa                                                            >  2:3023234/1‑100 (MQ=255)

                                     aTTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGaaaa                                                           >  2:3857261/1‑100 (MQ=255)

                                      tttGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGaaaa                                                           <  1:4129173/99‑1 (MQ=255)

                                      tttGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAg                                                          <  2:1527561/100‑1 (MQ=255)

                                      tttGAGGAATAATGATGAAAAAAAGTGCCATCACCGGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAg                                                          <  2:2635743/100‑1 (MQ=255)

                                        tGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAgcc                                                        >  2:1634591/1‑100 (MQ=255)

                                         gAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAgccg                                                       >  2:1259846/1‑100 (MQ=255)

                                            naataatGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCAGAAAAGCCGCCt                                                    >  1:456977/2‑100 (MQ=255)

                                            gAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCt                                                    >  1:1088338/1‑100 (MQ=255)

                                             aataatGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTg                                                   <  1:3649900/100‑1 (MQ=255)

                                             aataatGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTg                                                   <  1:3597200/100‑1 (MQ=255)

                                             aataatGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTg                                                   >  1:3164789/1‑100 (MQ=255)

                                                aataatGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGAt                                                <  2:1635276/100‑1 (MQ=255)

                                                    atgaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAg                                            <  1:2296314/100‑1 (MQ=255)

                                                     tgaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGt                                           <  2:1722483/100‑1 (MQ=255)

                                                     tgaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGt                                           >  1:1516484/1‑100 (MQ=255)

                                                      gaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTc                                          >  1:468301/1‑100 (MQ=255)

                                                      gaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTc                                          <  2:1196177/100‑1 (MQ=255)

                                                      gaAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTc                                          >  2:3013848/1‑100 (MQ=255)

                                                       aaaaaaaTTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTca                                         >  1:952320/1‑100 (MQ=255)

                                                        aaaaaaTTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTcag                                        >  2:2715974/1‑100 (MQ=255)

                                                        aaaaaaTTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTcag                                        <  1:187559/100‑1 (MQ=255)

                                                        aaaaaaTTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTcag                                        <  1:2138919/100‑1 (MQ=255)

                                                         aaaaaTTGCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTcagc                                       <  1:289248/100‑1 (MQ=255)

                                                                gCCATCACCTGTGCATTACTCTCAAGCTTGGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTcc                                >  1:3546275/1‑100 (MQ=255)

                                                                gCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGTGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTcc                                <  1:1361604/100‑1 (MQ=255)

                                                                gCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTcc                                >  1:2595369/1‑100 (MQ=255)

                                                                gCCATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCACAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTcc                                >  2:2359036/1‑100 (MQ=255)

                                                                 ccATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCa                               >  2:3906924/1‑100 (MQ=255)

                                                                 ccATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCa                               >  2:750740/1‑100 (MQ=255)

                                                                 ccATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTGTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCa                               <  2:2678948/100‑1 (MQ=255)

                                                                  cATCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCAc                              >  1:605973/1‑100 (MQ=255)

                                                                   aTCACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCACAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACg                             <  1:2359059/100‑1 (MQ=255)

                                                                     tctcTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGcc                           <  2:1755240/97‑1 (MQ=255)

                                                                     cACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGcc                           <  1:2173924/100‑1 (MQ=255)

                                                                     cACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGcc                           <  2:2198630/100‑1 (MQ=255)

                                                                     cACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGcc                           <  2:2062920/100‑1 (MQ=255)

                                                                     cACCTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGcc                           <  2:2859861/100‑1 (MQ=255)

                                                                       ccTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAg                         <  2:1381788/100‑1 (MQ=255)

                                                                       ccTGTGCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAg                         >  1:3612636/1‑100 (MQ=255)

                                                                         tgtgCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCt                       >  1:2141649/1‑100 (MQ=255)

                                                                          gtgCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCtt                      >  2:1318994/1‑100 (MQ=255)

                                                                           tgCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTc                     <  2:3400001/100‑1 (MQ=255)

                                                                            gCATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTca                    >  1:1769306/1‑100 (MQ=255)

                                                                             cATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTcac                   <  2:2921367/100‑1 (MQ=255)

                                                                             cATTACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTcac                   >  1:4175055/1‑100 (MQ=255)

                                                                               ttACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACACCAGCTTcacac                 <  1:2754182/100‑1 (MQ=255)

                                                                                tACTCTCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTcacaca                >  2:3313243/1‑100 (MQ=255)

                                                                                  ctctCAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACaaa              <  2:1769293/100‑1 (MQ=255)

                                                                                   tctcAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACaaaa             <  1:730941/100‑1 (MQ=255)

                                                                                    ctcAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACaaaaa            >  2:758038/1‑100 (MQ=255)

                                                                                    ctcAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACaaaaa            >  1:2409027/1‑100 (MQ=255)

                                                                                    ctcAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACaaaaa            >  1:3054513/1‑100 (MQ=255)

                                                                                      cAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGt          >  2:1506341/1‑100 (MQ=255)

                                                                                      cAAGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGt          >  1:683610/1‑100 (MQ=255)

                                                                                       aaGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGtg         >  2:3041243/1‑100 (MQ=255)

                                                                                        aGCTTAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGtg         >  2:1981236/1‑99 (MQ=255)

                                                                                           ttAGTACCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGtgactg     <  2:2378722/100‑1 (MQ=255)

                                                                                            tagtacCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGtgactg     >  2:1500572/1‑99 (MQ=255)

                                                                                             agtacCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGtgactgac   <  1:1609181/100‑1 (MQ=255)

                                                                                              gtacCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACg  >  2:3172730/1‑100 (MQ=255)

                                                                                              gtacCAAGCAGCGTTTGGGCTGATGCCGCAAGCCATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCTCGCCAGCTTCACACAAAAAGTGACTGACg  >  2:2517185/1‑100 (MQ=255)

TAGAGAGTAGAGGGAACTCCCGATCGGGAGTGACGTAATTTGAGGAATAATGATGAAAAAAATTGCCATCACCTGTGCATTACTCTCAAGCTTAGTAGCAAGCAGCGTTTGGGCTGATGCCGCAAGCGATCTGAAAAGCCGCCTGGATAAAGTCAGCAGCTTCCACGCCAGCTTCACACAAAAAGTGACTGACG  >  NC_000913/937320‑937513

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