BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,877,536	G→C	A87G (GCT→GGT)	b2755 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,877,536	0	G	C	85.9%	929.8 / 126.6	312	A87G (GCT→GGT)	b2755	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (24/20); new base C (140/128); total (164/148)
Fisher's exact test for biased strand distribution p-value = 8.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 312 total aligned reads displayed.

GTACGACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAACCCGTGTACC  >  NC_000913/2877447‑2877627

gtacgACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCAt                                                                                          >  2:4867215/1‑93 (MQ=255)

   cgACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATaaa                                                                                       <  1:2361417/93‑1 (MQ=255)

   cgACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATaaa                                                                                       >  1:934270/1‑93 (MQ=255)

   cgACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATaaa                                                                                       <  1:6315800/93‑1 (MQ=255)

    gACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAac                                                                                      <  1:1248355/93‑1 (MQ=255)

    gACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAac                                                                                      >  2:5985147/1‑93 (MQ=255)

    gACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAac                                                                                      >  1:3785059/1‑93 (MQ=255)

      ccTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAacac                                                                                    >  2:2518410/1‑93 (MQ=255)

      ccTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAacac                                                                                    <  2:1070804/93‑1 (MQ=255)

      ccTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAacac                                                                                    <  1:1572151/93‑1 (MQ=255)

       cTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACg                                                                                   >  1:427452/1‑93 (MQ=255)

        ttCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGa                                                                                  >  1:4768171/1‑93 (MQ=255)

         tCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGaa                                                                                 <  2:4991107/93‑1 (MQ=255)

         tCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGaa                                                                                 <  2:5405894/93‑1 (MQ=255)

           aGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACGAGAACCATAAACACGAACg                                                                               >  1:5722610/1‑93 (MQ=255)

                cAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGc                                                                          >  1:2422497/1‑93 (MQ=255)

                cAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGc                                                                          >  1:4604356/1‑93 (MQ=255)

                 aaaTCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCt                                                                         >  2:4559163/1‑93 (MQ=255)

                  aaTCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCtt                                                                        >  2:4225759/1‑93 (MQ=255)

                  aaTCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCtt                                                                        >  2:444626/1‑93 (MQ=255)

                   aTCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTc                                                                       >  2:4299416/1‑93 (MQ=255)

                    tCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTcc                                                                      >  2:3350537/1‑93 (MQ=255)

                    tCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTcc                                                                      >  2:3702120/1‑93 (MQ=255)

                     cTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTccc                                                                     <  1:1987733/93‑1 (MQ=255)

                     cTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTccc                                                                     >  2:662907/1‑93 (MQ=255)

                     cTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTccc                                                                     <  2:236447/93‑1 (MQ=255)

                     cTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTccc                                                                     >  1:5605025/1‑93 (MQ=255)

                     cTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTccc                                                                     <  1:4528243/93‑1 (MQ=255)

                      ttCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTcccc                                                                    <  2:4749036/93‑1 (MQ=255)

                       tCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCCCCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTccccc                                                                   <  1:7467847/93‑1 (MQ=255)

                       tCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTccccc                                                                   <  2:5614126/93‑1 (MQ=255)

                           ccAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCcccaccc                                                               <  1:3721100/93‑1 (MQ=255)

                            cAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccaccca                                                              <  1:1951117/93‑1 (MQ=255)

                            cAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCcccaccca                                                              <  2:1364796/93‑1 (MQ=255)

                            cAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCcccaccca                                                              <  2:5832312/93‑1 (MQ=255)

                             agagCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCAt                                                             >  1:172392/1‑93 (MQ=255)

                               agCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATAc                                                           <  2:3621705/93‑1 (MQ=255)

                                gCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATAcc                                                          >  1:1012986/1‑93 (MQ=255)

                                gCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATAcc                                                          >  1:3614244/1‑93 (MQ=255)

                                gCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATAcc                                                          >  1:7168080/1‑93 (MQ=255)

                                gCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATAcc                                                          >  1:6528138/1‑93 (MQ=255)

                                      tttGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAAt                                                    <  1:1299617/93‑1 (MQ=255)

                                         gCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCGGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGtt                                                 >  2:7254205/1‑93 (MQ=255)

                                         gCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGtt                                                 <  1:4319988/93‑1 (MQ=255)

                                          ccTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTc                                                <  2:690848/93‑1 (MQ=255)

                                          ccTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTc                                                <  1:2208159/93‑1 (MQ=255)

                                           cTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTcc                                               <  1:3424405/93‑1 (MQ=255)

                                           cTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTcc                                               <  1:5813789/93‑1 (MQ=255)

                                           cTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTcc                                               <  1:4717415/93‑1 (MQ=255)

                                            tGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCa                                              >  1:5193354/1‑93 (MQ=255)

                                            tGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCa                                              <  2:5091612/93‑1 (MQ=255)

                                            tGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCa                                              >  2:6854455/1‑93 (MQ=255)

                                             gATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCaa                                             <  1:3892319/93‑1 (MQ=255)

                                             gATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCaa                                             <  2:621864/93‑1 (MQ=255)

                                              aTAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAAc                                            <  1:2835743/93‑1 (MQ=255)

                                              aTAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAAc                                            <  2:2587802/93‑1 (MQ=255)

                                               tAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACt                                           <  2:3825145/93‑1 (MQ=255)

                                               tAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACt                                           >  1:590787/1‑93 (MQ=255)

                                                 gagCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTg                                         >  1:3261470/1‑93 (MQ=255)

                                                  agcagcTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTgc                                        <  1:4870906/93‑1 (MQ=255)

                                                     agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCa                                     <  2:3099265/93‑1 (MQ=255)

                                                     agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCa                                     <  2:741250/93‑1 (MQ=255)

                                                       cTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGc                                   >  2:1882564/1‑93 (MQ=255)

                                                       cTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGc                                   <  1:616403/93‑1 (MQ=255)

                                                       cTTATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGc                                   <  1:1819674/93‑1 (MQ=255)

                                                        ttATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGcc                                  >  2:3455595/1‑93 (MQ=255)

                                                        ttATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGcc                                  <  1:145656/93‑1 (MQ=255)

                                                         tATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCa                                 <  2:1722449/93‑1 (MQ=255)

                                                         tATCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCa                                 >  2:3549821/1‑93 (MQ=255)

                                                          aTCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAg                                >  1:7526677/1‑93 (MQ=255)

                                                           tCTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAgg                               <  2:6304815/93‑1 (MQ=255)

                                                            cTGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGc                              >  2:5703737/1‑93 (MQ=255)

                                                             tGAACGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCg                             >  1:6687531/1‑93 (MQ=255)

                                                               aaCGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTa                           >  2:2006829/1‑93 (MQ=255)

                                                               aaCGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTa                           <  1:7040483/93‑1 (MQ=255)

                                                                aCGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTAc                          >  1:4088847/1‑93 (MQ=255)

                                                                aCGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTAc                          <  2:5131648/93‑1 (MQ=255)

                                                                aCGCGCACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTAc                          <  1:5408623/93‑1 (MQ=255)

                                                                  gcgcACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAg                        >  2:6457535/1‑93 (MQ=255)

                                                                   cgcACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGc                       <  2:3132631/93‑1 (MQ=255)

                                                                    gcACCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCt                      >  1:268325/1‑93 (MQ=255)

                                                                      aCCTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGc                    >  2:2758524/1‑93 (MQ=255)

                                                                       ccTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCCGCCCGGCGTACAGCTGCa                   >  1:3006555/1‑93 (MQ=255)

                                                                       ccTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCa                   >  1:2612531/1‑93 (MQ=255)

                                                                       ccTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCa                   <  1:5561069/93‑1 (MQ=255)

                                                                        cTCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCAt                  >  2:1489152/1‑93 (MQ=255)

                                                                         tCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATg                 >  1:3474021/1‑93 (MQ=255)

                                                                         tCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATg                 >  1:3332313/1‑93 (MQ=255)

                                                                         tCCAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATg                 >  2:2029672/1‑93 (MQ=255)

                                                                          ccAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGc                <  1:3947035/93‑1 (MQ=255)

                                                                          ccAGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGc                <  2:2736382/93‑1 (MQ=255)

                                                                            aGGCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGa              >  1:7355930/1‑93 (MQ=255)

                                                                             ggCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCTCCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGaa             <  2:2459143/93‑1 (MQ=255)

                                                                              gCTGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGaaa            >  1:4515510/1‑93 (MQ=255)

                                                                               cTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAAc           >  1:7182428/1‑93 (MQ=255)

                                                                                tGACCAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAAcc          <  2:5545631/93‑1 (MQ=255)

                                                                                    cAGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAACCCgtg      >  2:5172259/1‑93 (MQ=255)

                                                                                     aGAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAACCCgtgt     >  1:5029414/1‑93 (MQ=255)

                                                                                      gAACCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGTGTACAGCTGCATGCGAAACCCGTGTa    <  2:1593161/93‑1 (MQ=255)

                                                                                        aGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAACCCGTGTAcc  <  2:7322892/93‑1 (MQ=255)

GTACGACCTTCAGACGCAAATCTTCATCCAGAGCAAGTTTTGCCTGATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAACCCGTGTACC  >  NC_000913/2877447‑2877627

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