BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,736,609	G→A	21.8%	P190L (CCA→CTA)	bax ←	putative glucosaminidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,736,609	0	G	A	21.8%	253.3 / 63.4	130	P190L (CCA→CTA)	bax	putative glucosaminidase
Reads supporting (aligned to +/- strand): ref base G (53/48); new base A (17/11); total (71/59)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 134 total aligned reads displayed.

CGGAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCTGCAGCCA  >  NC_000913/3736510‑3736708

cGGAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTg                                                                                                     <  2:3380670/100‑1 (MQ=255)

 ggAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTgg                                                                                                    >  2:2564169/1‑100 (MQ=255)

  gAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGgc                                                                                                   <  1:789362/100‑1 (MQ=255)

    aCGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGacgc                                                                                                 >  1:2771558/1‑100 (MQ=255)

     cGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCa                                                                                                <  2:2116341/100‑1 (MQ=255)

     cGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTCGCGCa                                                                                                >  1:117349/1‑100 (MQ=255)

      gAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCCCCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCAt                                                                                               <  1:2364408/100‑1 (MQ=255)

      gAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCAt                                                                                               >  1:2598733/1‑100 (MQ=255)

         gagTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTgg                                                                                            >  1:1862867/1‑100 (MQ=255)

         gagTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTgg                                                                                            <  2:165122/100‑1 (MQ=255)

             aaGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACa                                                                                        >  1:2542871/1‑100 (MQ=255)

             aaGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACa                                                                                        >  1:1792760/1‑100 (MQ=255)

              aGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACaa                                                                                       >  2:2087074/1‑100 (MQ=255)

                ccGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACg                                                                                     <  1:717660/100‑1 (MQ=255)

                ccGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACg                                                                                     <  2:2705068/100‑1 (MQ=255)

                  gggTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTGACTTTACCTGGCGCATTGGTACAACGTc                                                                                   >  2:155154/1‑100 (MQ=255)

                   ggTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTcc                                                                                  >  2:1388224/1‑100 (MQ=255)

                   ggTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTcc                                                                                  <  1:598308/100‑1 (MQ=255)

                     tGCGTATTCAGGTTATTTATATAGTCGTTCACTGATTCTTTGACAGAACTAAATTGTGAGTACCTTTTCATTTTATCTGGCGTATTGGTACAACg                                                                                     <  1:208672/95‑1 (MQ=255)

                      gCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCttt                                                                               >  1:1682242/1‑100 (MQ=255)

                       ngTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTc                                                                              >  1:2386162/2‑100 (MQ=255)

                       cGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTc                                                                              <  1:3299070/100‑1 (MQ=255)

                        gTATTCAGGTTAGTGACATAGGCGCTCACCTATTCTTTGACAGAACTAATCTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCa                                                                             <  1:3163050/100‑1 (MQ=255)

                         tATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCAt                                                                            >  2:3274859/1‑100 (MQ=255)

                         tATTCAGCGTAGTGACATAGGCGCTCACCAATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCAt                                                                            <  1:287916/100‑1 (MQ=255)

                          aTTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATg                                                                           <  1:2134546/100‑1 (MQ=255)

                          aTTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATg                                                                           <  2:1526847/100‑1 (MQ=255)

                             cAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCAt                                                                        >  2:162852/1‑100 (MQ=255)

                               ggTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCAttt                                                                      >  2:302055/1‑100 (MQ=255)

                                 ttAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTCCCCTTTCACTTTACCTGGCGCATTGGTACTACGTCCTTTCATGCATTTCa                                                                    <  1:1738728/100‑1 (MQ=255)

                                 ttAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCa                                                                    <  1:2199625/100‑1 (MQ=255)

                                 ttAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCa                                                                    >  1:2984624/1‑100 (MQ=255)

                                  tAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCAt                                                                   <  1:944049/100‑1 (MQ=255)

                                    gTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGc                                                                 <  1:1610718/100‑1 (MQ=255)

                                    gTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGc                                                                 >  1:1992841/1‑100 (MQ=255)

                                    gTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGc                                                                 >  2:1200335/1‑100 (MQ=255)

                                    gTGACATAGGCGCTAACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGAACTGGTACAACGTCCTTTCATGCATTTCATGc                                                                 >  1:1187808/1‑100 (MQ=255)

                                     tGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGcc                                                                <  2:767837/100‑1 (MQ=255)

                                      gACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCg                                                               >  1:3027170/1‑100 (MQ=255)

                                        cATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGaa                                                             >  1:1784856/1‑100 (MQ=255)

                                        cATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGaa                                                             >  2:1326171/1‑100 (MQ=255)

                                          tAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACa                                                           >  1:364652/1‑100 (MQ=255)

                                          tAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACa                                                           <  2:1477800/100‑1 (MQ=255)

                                           aGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAg                                                          <  1:520814/100‑1 (MQ=255)

                                            ggCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGACCAgg                                                         >  2:2192515/1‑100 (MQ=255)

                                             gcgcTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgt                                                        <  1:2297159/100‑1 (MQ=255)

                                                cTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgt                                                     <  2:441636/100‑1 (MQ=255)

                                                  cACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                                   <  2:876801/100‑1 (MQ=255)

                                                  cACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                                   >  2:2598274/1‑100 (MQ=255)

                                                  cACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                                   <  2:2686616/100‑1 (MQ=255)

                                                  cACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                                   >  1:157685/1‑100 (MQ=255)

                                                  cACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                                   >  1:619792/1‑100 (MQ=255)

                                                  cACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                                   <  2:1859903/100‑1 (MQ=255)

                                                   aCCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgt                                                  <  2:308356/100‑1 (MQ=255)

                                                    ctgATTCTTTGACATAACTAAGCTGTGAGTACCCTTTCACTTTACCTGGCGCTTTGGTACAACGTCCTTTCATGCATTTTATGCCGAACTGgttgttgt                                                  <  1:216954/97‑1 (MQ=255)

                                                    ccGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgtt                                                 <  2:2124515/100‑1 (MQ=255)

                                                       aTTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACTGgttgttgttgtt                                              <  1:2090192/100‑1 (MQ=255)

                                                        ttCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGttgttgttgttg                                             <  1:1570018/100‑1 (MQ=255)

                                                        ttCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGttgttgttgttg                                             <  1:3349092/100‑1 (MQ=255)

                                                         ncTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgc                                            >  2:2662717/2‑100 (MQ=255)

                                                          cTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcg                                           >  1:758883/1‑100 (MQ=255)

                                                          cTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCAATTCATGCCGAACAGCttgttgttgttg                                             >  1:1635584/1‑98 (MQ=255)

                                                           tttGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTACTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgc                                          >  2:2164278/1‑100 (MQ=255)

                                                             tGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgcgc                                        >  1:2583284/1‑100 (MQ=255)

                                                              gACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgcgc                                        >  1:581115/1‑99 (MQ=255)

                                                              gACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGcc                                       <  2:3324132/100‑1 (MQ=255)

                                                              gACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGcc                                       <  2:2172927/100‑1 (MQ=255)

                                                              gACAGAACTAAACTGTGAGTACAATTTCACTCTACCTAGCGCATTCGTACAACGTCCTTTCATTCATTTCGTGCCGGACATTGTGTTGTTGTTgcg                                           >  1:1123528/1‑96 (MQ=255)

                                                               aCAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTTGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCa                                      >  1:2893214/1‑100 (MQ=255)

                                                               aCAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCa                                      >  2:799541/1‑100 (MQ=255)

                                                               aCAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAAAGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgcgc                                        >  1:1117592/1‑98 (MQ=255)

                                                                 aGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGc                                    <  2:1050853/100‑1 (MQ=255)

                                                                  gAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCt                                   >  2:2318317/1‑100 (MQ=255)

                                                                     cTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTTCGCGCCAGCTTCg                                <  2:2276482/100‑1 (MQ=255)

                                                                     cTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCg                                >  2:563997/1‑100 (MQ=255)

                                                                      tAACCTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTc                                 <  2:3367171/98‑1 (MQ=255)

                                                                      tAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCATTTCATGCATTTCCTGCCGAACAGGTTGTTGTTGTTGCCCGACAGCTTCGa                               >  1:630263/1‑100 (MQ=255)

                                                                       aaaCTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGAc                              <  1:2738554/100‑1 (MQ=255)

                                                                       aaaCTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGAc                              <  1:1939975/100‑1 (MQ=255)

                                                                         aCTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGt                            <  1:3028920/100‑1 (MQ=255)

                                                                           tgtgAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTc                          >  1:17104/1‑100 (MQ=255)

                                                                            gtgAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTcc                         <  1:143316/100‑1 (MQ=255)

                                                                              gAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTcccc                       >  2:1899992/1‑100 (MQ=255)

                                                                                gTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCaa                     >  1:3271974/1‑100 (MQ=255)

                                                                                 tACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAAc                    >  1:466437/1‑100 (MQ=255)

                                                                                 tACCCTTTCACTTTACCTGGCGCATTGGCACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAAc                    <  2:1288672/100‑1 (MQ=255)

                                                                                   cccTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCg                  <  2:2525445/100‑1 (MQ=255)

                                                                                   cccTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCg                  >  2:1811786/1‑100 (MQ=255)

                                                                                     cTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCt                >  2:1023412/1‑100 (MQ=255)

                                                                                     cTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCt                >  1:2390410/1‑100 (MQ=255)

                                                                                        tCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTc             <  2:1378558/100‑1 (MQ=255)

                                                                                         cACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTct            <  2:1634691/100‑1 (MQ=255)

                                                                                         cACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTct            >  2:3201659/1‑100 (MQ=255)

                                                                                           cTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctgc          <  1:2880075/100‑1 (MQ=255)

                                                                                             ttACCTGGCGCATTGGTACAACGTGCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctgctg        <  1:146102/100‑1 (MQ=255)

                                                                                             ttACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctgctg        >  2:1097632/1‑100 (MQ=255)

                                                                                              tACCTGGCGCATTGGTACAACGTCCTTTCATGCATCCCATGCCCAACAGGATGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCtgctgc       >  1:256536/1‑100 (MQ=255)

                                                                                              tACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCtgctgc       <  2:684941/100‑1 (MQ=255)

                                                                                                ccTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCTGCAg     >  2:2920033/1‑100 (MQ=255)

                                                                                                  nggCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCTGCAGcc   >  2:2098730/2‑100 (MQ=255)

CGGAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCTGCAGCCA  >  NC_000913/3736510‑3736708

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