BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,544	C→A	100%	S13I (AGC→ATC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,544	0	C	A	97.3%	547.3 / 2.0	163	S13I (AGC→ATC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (3/2); new base A (82/75); total (86/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

CGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCG  >  NC_000913/3566447‑3566640

  ttGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAg                                                                                              >  2:2228027/1‑100 (MQ=255)

  ttGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAg                                                                                              <  1:3498089/100‑1 (MQ=255)

     cGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGt                                                                                           >  2:1097695/1‑100 (MQ=255)

     cGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGt                                                                                           >  2:796739/1‑100 (MQ=255)

      gTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTg                                                                                          >  1:1036649/1‑100 (MQ=255)

       ttCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTgg                                                                                         >  2:3887352/1‑100 (MQ=255)

         cAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGc                                                                                       >  2:1623138/1‑100 (MQ=255)

           gCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGa                                                                                     >  2:544883/1‑100 (MQ=255)

           gCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGa                                                                                     <  2:2837299/100‑1 (MQ=255)

           gCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGa                                                                                     >  1:1799147/1‑100 (MQ=255)

            ccATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGAt                                                                                    <  1:3525448/100‑1 (MQ=255)

            ccATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGAt                                                                                    <  2:3835582/100‑1 (MQ=255)

            ccATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGAt                                                                                    >  1:1814995/1‑100 (MQ=255)

             cATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGATAAGCGTGGGCGATg                                                                                   <  2:3856667/100‑1 (MQ=255)

             cATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATg                                                                                   >  1:3463921/1‑100 (MQ=255)

             cATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATg                                                                                   <  1:3338092/100‑1 (MQ=255)

             cATTCATGTTTATTGGCGACGACCGGGTCCGTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATg                                                                                   <  2:1791952/100‑1 (MQ=255)

              aTTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGa                                                                                  <  2:1425433/100‑1 (MQ=255)

                 cATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAata                                                                               >  2:3774492/1‑100 (MQ=255)

                   tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                             <  2:1613303/100‑1 (MQ=255)

                   tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                             <  2:844612/100‑1 (MQ=255)

                   tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                             <  1:3703125/100‑1 (MQ=255)

                   tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                             >  1:1976363/1‑100 (MQ=255)

                   tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                             >  2:3032695/1‑100 (MQ=255)

                     tttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTa                                                                           <  2:3377297/100‑1 (MQ=255)

                       tATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTaaa                                                                         <  1:228597/100‑1 (MQ=255)

                        aTTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAAc                                                                        <  1:1559618/100‑1 (MQ=255)

                         ttGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACg                                                                       >  2:3153767/1‑100 (MQ=255)

                         ttGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACg                                                                       <  2:1243332/100‑1 (MQ=255)

                          tGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACgn                                                                      <  2:3220094/100‑2 (MQ=255)

                          tGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACgg                                                                      >  1:98465/1‑100 (MQ=255)

                           ggCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTATGATAAGCGTGGGCGATGAATATGTAAACGGa                                                                     >  2:217242/1‑100 (MQ=255)

                           ggCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGa                                                                     >  2:2641836/1‑100 (MQ=255)

                            gCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGa                                                                     >  2:3316961/1‑99 (MQ=255)

                             cgacgaCCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGc                                                                   >  1:1607248/1‑100 (MQ=255)

                             cgacgaCCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGc                                                                   <  1:2122262/100‑1 (MQ=255)

                                    cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCata                                                            <  1:3316983/100‑1 (MQ=255)

                                    cGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCata                                                            <  1:1691545/100‑1 (MQ=255)

                                     gggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCAATGAATATGTAAACGGAGCATTCatat                                                           >  1:1286002/1‑100 (MQ=255)

                                        tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAgg                                                        >  1:1711471/1‑100 (MQ=255)

                                         ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGCGCGATGCATCTGTACACGGAGCATTCATATAGGc                                                       >  1:374591/1‑100 (MQ=255)

                                         ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGc                                                       <  1:627173/100‑1 (MQ=255)

                                         ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGc                                                       >  2:1977149/1‑100 (MQ=255)

                                             tCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGttt                                                   <  1:816182/100‑1 (MQ=255)

                                              ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTc                                                  <  2:518531/100‑1 (MQ=255)

                                               cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGGTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGACGTTTcc                                                 >  1:770476/1‑100 (MQ=255)

                                               cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                                 >  1:54298/1‑100 (MQ=255)

                                               cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                                 <  2:20290/100‑1 (MQ=255)

                                               cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                                 <  1:276784/100‑1 (MQ=255)

                                                 aTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTg                                               <  2:2946697/100‑1 (MQ=255)

                                                  tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                              <  2:1976353/100‑1 (MQ=255)

                                                   cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                             >  2:1860163/1‑100 (MQ=255)

                                                   cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                             <  1:1042040/100‑1 (MQ=255)

                                                   cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                             <  1:2931589/100‑1 (MQ=255)

                                                   cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                             >  1:3247478/1‑100 (MQ=255)

                                                   cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                              <  1:1795043/99‑1 (MQ=255)

                                                    gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaa                                            >  2:261020/1‑100 (MQ=255)

                                                     tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaaa                                           <  1:3113418/100‑1 (MQ=255)

                                                     tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaaa                                           >  2:2648964/1‑100 (MQ=255)

                                                      aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                                          <  2:1814154/100‑1 (MQ=255)

                                                      aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                                          >  2:2685121/1‑100 (MQ=255)

                                                         cATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                                       <  2:2060590/100‑1 (MQ=255)

                                                           tCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                                     >  2:2716002/1‑100 (MQ=255)

                                                            cAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                                    <  2:3524985/100‑1 (MQ=255)

                                                             aGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa                                   <  1:2497985/100‑1 (MQ=255)

                                                             aGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCa                                   >  1:1355448/1‑100 (MQ=255)

                                                              gCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                                  <  2:2875064/100‑1 (MQ=255)

                                                                ttGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc                                <  2:1535950/100‑1 (MQ=255)

                                                                ttGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc                                <  1:1446830/100‑1 (MQ=255)

                                                                 tGTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCg                               >  1:2811645/1‑100 (MQ=255)

                                                                  gTAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGa                              >  1:2571246/1‑100 (MQ=255)

                                                                   tAAGCTATAGTGTGCTTAAGAGCTTCTACGATAAGCGTCGGCGCTGAATATGTAAACGGAGCATTCATATAGGTGTTTCCTGAAAACTATTTCAAGCGAt                             <  2:3193367/100‑1 (MQ=255)

                                                                   tAAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt                             <  2:140285/100‑1 (MQ=255)

                                                                     aGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAg                           >  1:2153015/1‑100 (MQ=255)

                                                                        gATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaa                        >  1:3845217/1‑100 (MQ=255)

                                                                          tAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg                      <  2:1492864/100‑1 (MQ=255)

                                                                            gagTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCt                    >  2:527766/1‑100 (MQ=255)

                                                                              gTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                  >  1:1483724/1‑100 (MQ=255)

                                                                              gTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                  >  1:69812/1‑100 (MQ=255)

                                                                               tGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCAc                 <  2:1644959/100‑1 (MQ=255)

                                                                                gCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg                >  1:3579981/1‑100 (MQ=255)

                                                                                gCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg                <  2:1576221/100‑1 (MQ=255)

                                                                                gCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg                >  2:3517273/1‑100 (MQ=255)

                                                                                 cTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg               >  2:966295/1‑100 (MQ=255)

                                                                                 cTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg               >  2:3017169/1‑100 (MQ=255)

                                                                                 cTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg               <  2:2927969/100‑1 (MQ=255)

                                                                                 cTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg               >  1:78698/1‑100 (MQ=255)

                                                                                  ttAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGt              >  1:249348/1‑100 (MQ=255)

                                                                                   tAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTa             <  1:1977159/100‑1 (MQ=255)

                                                                                   tAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTa             >  1:1535847/1‑100 (MQ=255)

                                                                                   tAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTa             <  1:8773/100‑1 (MQ=255)

                                                                                     agagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg           >  1:1596461/1‑100 (MQ=255)

                                                                                       agCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGAc         <  1:3904337/100‑1 (MQ=255)

                                                                                        gCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACt        >  2:2514489/1‑100 (MQ=255)

                                                                                        gCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACt        <  2:1174508/100‑1 (MQ=255)

                                                                                          ttCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTc      <  2:3703720/100‑1 (MQ=255)

                                                                                            cTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgc    >  1:438764/1‑100 (MQ=255)

                                                                                              aCGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgccg  >  2:3471363/1‑100 (MQ=255)

                                                                                              aCGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgccg  >  2:581765/1‑100 (MQ=255)

                                                                                              aCGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgccg  <  1:1213316/100‑1 (MQ=255)

CGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCG  >  NC_000913/3566447‑3566640

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: `‑`