BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,736,609	G→A	43.3%	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	43.3%	43.2 / 112.0	100	P190L (CCA→CTA)	bax	putative glucosaminidase
Reads supporting (aligned to +/- strand): ref base G (30/26); new base A (17/26); total (47/53)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

Only 100 of 103 total aligned reads displayed.

GGAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCTGCAGCCA  >  NC_000913/3736511‑3736708

ggAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTgg                                                                                                    <  2:945360/100‑1 (MQ=255)

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

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

  aaCGAAGAGTAAGCCGGGTGCGTATTCAGGTTATTGACATAGGCTCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTTgcg                                                                                                  <  2:1925229/100‑1 (MQ=255)

   aCGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGgcgc                                                                                                 >  1:1360124/1‑100 (MQ=255)

   aCGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGgcgc                                                                                                 >  1:1473406/1‑100 (MQ=255)

   aCGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAgcgc                                                                                                 >  2:2027452/1‑100 (MQ=255)

    cGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCa                                                                                                >  2:1772332/1‑100 (MQ=255)

        gagTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTATCTAGCGCATTgg                                                                                            <  1:1573135/100‑1 (MQ=255)

         agTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGt                                                                                           <  2:219499/100‑1 (MQ=255)

          gTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTa                                                                                          <  2:2126388/100‑1 (MQ=255)

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

             aGCCGGGTGCGTATTCAGGATAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACa                                                                                        >  2:1141420/1‑99 (MQ=255)

              gCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAAc                                                                                      >  1:2112081/1‑100 (MQ=255)

              gCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAAc                                                                                      <  1:2111268/100‑1 (MQ=255)

               ccGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACg                                                                                     >  1:1446891/1‑100 (MQ=255)

               ccGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACg                                                                                     <  2:1032583/100‑1 (MQ=255)

                  ggTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACATAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTcc                                                                                  <  2:1200544/100‑1 (MQ=255)

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

                      cGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTc                                                                              >  2:463403/1‑100 (MQ=255)

                      cGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTc                                                                              <  1:1104956/100‑1 (MQ=255)

                       gTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCa                                                                             <  2:1553540/100‑1 (MQ=255)

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

                         aTTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACATTACCTGGCGCATTGGTACAACGTCCTTTCATg                                                                           <  2:1432562/100‑1 (MQ=255)

                          ttCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCACTGGTACAACGTCCTTTCATGc                                                                          >  1:1672667/1‑100 (MQ=255)

                           tCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCa                                                                         <  1:56919/100‑1 (MQ=255)

                            cAGGTTAGTGTCATAGGCGCTCACCGATTCTTTGACAGAACTAATCTGTGAGTACCCTTTCACTTTGCCTGGCGCATTGGTACAACGTCCTTTCATGCAt                                                                        <  2:1578465/100‑1 (MQ=255)

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

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

                              ggTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCAttt                                                                      <  1:1873800/100‑1 (MQ=255)

                                 tAGTGCCATAGGCGCTCCTCGTTCCTTTGATAGAACTAAACTTTGGGTAACCTTTCACTTCACCTGGCGCATTGGTACAACGTCCTTAAATGCa                                                                         <  2:2262429/94‑1 (MQ=255)

                                  aGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATg                                                                  >  1:115924/1‑100 (MQ=255)

                                  aGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATg                                                                  <  1:1540141/100‑1 (MQ=255)

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

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

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

                                    tGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGcc                                                                <  2:1723060/100‑1 (MQ=255)

                                     gACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCg                                                               <  2:845987/100‑1 (MQ=255)

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

                                      aCATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGa                                                              >  2:2133727/1‑100 (MQ=255)

                                       cATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGaa                                                             >  2:1692608/1‑100 (MQ=255)

                                        aTAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAAc                                                            >  2:191711/1‑100 (MQ=255)

                                        aTAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAAc                                                            <  1:433422/100‑1 (MQ=255)

                                         tAGGCGCCCACTGCTTCTTTAACAGAACTAAACGGTGAGTACCCTATCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACa                                                           >  2:1585618/1‑100 (MQ=255)

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

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

                                            gcgcTCACCGATTCTTTGACAGAACTAAACTGCGACTACCCTTTCACTTTACCTGGCGCATTGGTAAAACGTCATTTCATGCATTTCATGCCGAACAGgt                                                        <  2:1395491/100‑1 (MQ=255)

                                             cgcTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgtt                                                       >  1:1753286/1‑100 (MQ=255)

                                             cgcTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgt                                                        >  2:421338/1‑99 (MQ=255)

                                             cgcTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAgg                                                         >  1:56561/1‑98 (MQ=255)

                                               cTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgt                                                     >  2:340345/1‑100 (MQ=255)

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

                                               cTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCAAACAGgttgt                                                     >  1:120944/1‑100 (MQ=255)

                                                tCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgtt                                                    >  2:364566/1‑100 (MQ=255)

                                                  aCCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgt                                                  <  2:1309828/100‑1 (MQ=255)

                                                  aCCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgt                                                  <  2:1377757/100‑1 (MQ=255)

                                                  aCCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                                   <  1:2065549/99‑1 (MQ=255)

                                                   ccGCTTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgtt                                                 <  2:1105287/100‑1 (MQ=255)

                                                   ccGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgtt                                                 >  2:68367/1‑100 (MQ=255)

                                                   ccGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgtt                                                 >  1:773905/1‑100 (MQ=255)

                                                    cGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgttg                                                <  2:1139213/100‑1 (MQ=255)

                                                     gATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgttgt                                               <  2:138284/100‑1 (MQ=255)

                                                     gATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgttgt                                               <  2:1497329/100‑1 (MQ=255)

                                                      aTTCTTTGACAGAACTAAACGGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgttgtt                                              >  1:2259790/1‑100 (MQ=255)

                                                       ttCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgttgtt                                              >  2:1738752/1‑99 (MQ=255)

                                                       ttCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgttgt                                               <  1:2132424/98‑1 (MQ=255)

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

                                                         cTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTATTgcg                                           >  1:433174/1‑100 (MQ=255)

                                                          tttGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgc                                          >  1:1117727/1‑100 (MQ=255)

                                                           ttGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgcg                                         <  1:1921557/100‑1 (MQ=255)

                                                            tGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgcgc                                        <  1:1009615/100‑1 (MQ=255)

                                                               cAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAg                                     <  2:873037/100‑1 (MQ=255)

                                                                   aCTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCt                                   >  2:2059843/1‑98 (MQ=255)

                                                                   aCTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTc                                 >  1:2193456/1‑100 (MQ=255)

                                                                    cTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCg                                >  1:2267340/1‑100 (MQ=255)

                                                                    cTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCg                                <  2:172035/100‑1 (MQ=255)

                                                                    cTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCg                                <  1:1047108/100‑1 (MQ=255)

                                                                     tAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCg                                <  1:268827/99‑1 (MQ=255)

                                                                     tAAACTGTGAGTACCCTTTCACTTTAACTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGa                               >  2:789465/1‑100 (MQ=255)

                                                                      aaaCTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGAc                              >  2:1028238/1‑100 (MQ=255)

                                                                      aaaCTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGAc                              >  1:974365/1‑100 (MQ=255)

                                                                      aaaCTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGAc                              >  2:49141/1‑100 (MQ=255)

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

                                                                           gtgAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTcc                         >  2:192475/1‑100 (MQ=255)

                                                                           gtgAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGCACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTcc                         <  2:1283948/100‑1 (MQ=255)

                                                                            ctaGTACCCTTTCACTTTACCTGGCGCATTGGTACATCGTCCTTTCATGCATTTCATGCCGAACATGTTGTTGTTGTTCCGCGCCAGCTTCGACGTTcc                         <  2:1457325/97‑1 (MQ=255)

                                                                               gTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCaa                     >  1:1883061/1‑100 (MQ=255)

                                                                                tACCCTTTCACTTTACCTGGCGCAGTGGTCCAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAAc                    >  1:14769/1‑100 (MQ=255)

                                                                                tACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAAc                    <  1:1874094/100‑1 (MQ=255)

                                                                                tACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAAc                    <  1:910056/100‑1 (MQ=255)

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

                                                                                      ttCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCttt              <  2:1242879/100‑1 (MQ=255)

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

                                                                                        cACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACTTTCCCCAACCGCTTTct            <  2:374022/100‑1 (MQ=255)

                                                                                        cACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTct            >  1:530564/1‑100 (MQ=255)

                                                                                         ttccgatcTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTGGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctg           <  2:1473816/93‑1 (MQ=255)

                                                                                         aCTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctg           >  1:603987/1‑100 (MQ=255)

                                                                                          cattACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctgc          <  2:1446875/98‑1 (MQ=255)

                                                                                               ccTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGAGCGCCAGCTTCCACGTTCCCCAACCGCTTTCTGCTGCAg     >  2:1183683/1‑100 (MQ=255)

                                                                                                  ggCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctgagatcggaa  >  1:1473832/1‑91 (MQ=255)

GGAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCTGCAGCCA  >  NC_000913/3736511‑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: `‑`