BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	157,879	C→A	5.1%	L416L (CTG→CTT)	pcnB ←	poly(A) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	157,879	0	C	A	5.1%	346.6 / 4.8	107	L416L (CTG→CTT)	pcnB	poly(A) polymerase
Reads supporting (aligned to +/- strand): ref base C (40/61); new base A (4/2); total (44/63)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTAATGTCGT  >  NC_000913/157732‑158028

tGCGGTACCCTCAGGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAg                                                                                                                                                     <  2:1018167/150‑1 (MQ=255)

 gCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGt                                                                                                                                                    <  2:1542390/150‑1 (MQ=255)

   gtaaCCCTGACGACGTGGTGCTCGTTTGCGTGGACGACGAGTACGACGACTCGGTGAGGGTTCTTCATCCAGCTCGTTGATCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTc                                                                                                                                                  <  1:302074/147‑1 (MQ=255)

   ggTCCCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGGGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTc                                                                                                                                                  <  1:229065/150‑1 (MQ=255)

   ggTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTc                                                                                                                                                  >  1:554386/1‑150 (MQ=255)

    gTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCa                                                                                                                                                 <  1:1540474/150‑1 (MQ=255)

      aCCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGCAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGc                                                                                                                                               >  1:1449230/1‑150 (MQ=255)

        ccTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCCCCATTTCACCAGACGCTGCAGTTCAGCGt                                                                                                                                             >  2:139992/1‑150 (MQ=255)

         cTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGtt                                                                                                                                            >  1:1473156/1‑150 (MQ=255)

           cACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTAc                                                                                                                                          >  2:1215720/1‑150 (MQ=255)

           cACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGAACTGCAGCGtctg                                                                                                                                          <  1:561518/150‑4 (MQ=255)

               acgTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTc                                                                                                                                      <  2:163213/150‑1 (MQ=255)

               acgTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTc                                                                                                                                      <  2:466053/150‑1 (MQ=255)

                  tggtgCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAAc                                                                                                                                   <  2:799891/150‑1 (MQ=255)

                   ggtgCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACt                                                                                                                                  >  2:1474862/1‑150 (MQ=255)

                      gcgcGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCa                                                                                                                               <  2:1197796/150‑1 (MQ=255)

                      gcgcGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCa                                                                                                                               <  1:1679533/150‑1 (MQ=255)

                       cgcgTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAg                                                                                                                              <  1:690139/150‑1 (MQ=255)

                       cgcgTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAg                                                                                                                              <  1:794630/150‑1 (MQ=255)

                       cgcgTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAg                                                                                                                              <  1:947660/150‑1 (MQ=255)

                         cgTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCt                                                                                                                            >  2:18114/1‑150 (MQ=255)

                         cgTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCGCAACTTCAGcg                                                                                                                            >  1:1250352/1‑149 (MQ=255)

                          gtgtGCGTGGACGACGAGTACGAGGACGCGGGGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCA‑CCCACCATTTCTCCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCg                                                                                                                          <  1:251119/147‑1 (MQ=255)

                          gTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTc                                                                                                                           <  2:860370/150‑1 (MQ=255)

                          gTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCATTTTTGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTTCAGTTCAGCGTTACTCTCAACTTCAGCTc                                                                                                                           <  1:667330/150‑1 (MQ=255)

                          gTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTGCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCA‑CCCACCTTTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCtgg                                                                                                                          <  1:251120/150‑3 (MQ=255)

                            ttGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGc                                                                                                                         <  1:793002/150‑1 (MQ=255)

                              gCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGACCTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCaa                                                                                                                       >  2:100719/1‑150 (MQ=255)

                                gTGGACGACGAGTACGACGCCGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAgg                                                                                                                     >  2:1475340/1‑150 (MQ=255)

                                  ggACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGcc                                                                                                                   <  2:1610131/150‑1 (MQ=255)

                                  ggACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGcc                                                                                                                   <  1:404284/150‑1 (MQ=255)

                                   gacgacGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCa                                                                                                                  >  2:996294/1‑150 (MQ=255)

                                           gTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGACCTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCACCAGGGCa                                                                                                          >  1:841209/1‑150 (MQ=255)

                                           gTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCa                                                                                                          <  1:741017/150‑1 (MQ=255)

                                           gTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCa                                                                                                          <  2:554386/150‑1 (MQ=255)

                                           gTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCa                                                                                                          >  2:725441/1‑150 (MQ=255)

                                                acgacgCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGc                                                                                                     <  1:1215720/150‑1 (MQ=255)

                                                  gacgCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCg                                                                                                   >  1:918441/1‑150 (MQ=255)

                                                  gacgCGGGGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCg                                                                                                   <  1:1182646/150‑1 (MQ=255)

                                                   acgCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGccac                                                                                                  <  1:657641/150‑3 (MQ=255)

                                                     gcgGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCAc                                                                                                <  1:32838/150‑1 (MQ=255)

                                                       ggTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACgg                                                                                              <  1:1654029/150‑1 (MQ=255)

                                                            cGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACtt                                                                                         >  2:523676/1‑150 (MQ=255)

                                                            cGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACtt                                                                                         >  1:1350323/1‑150 (MQ=255)

                                                             ggTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTa                                                                                        <  2:1589907/150‑1 (MQ=255)

                                                              gTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCACCAGGTCATAAGCCGCACGGAACTTAg                                                                                       >  1:1143130/1‑150 (MQ=255)

                                                              gTTCTTCATCCAGCGCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAg                                                                                       <  2:448502/150‑1 (MQ=255)

                                                                 cttcATCCAGCTCGTTGCGCATCCGTTTTTGGTCTGGTGGCGCGGAGACCTGGACCTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACGGGTCATAAGCCGCACGGAACTTAGGAt                                                                                    >  2:271021/1‑150 (MQ=255)

                                                                 cttcATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGAt                                                                                    <  1:750552/150‑1 (MQ=255)

                                                                   tcATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGc                                                                                  <  1:725441/150‑1 (MQ=255)

                                                                    cATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACACCACCATTTCACCAGACGCTGCAGTTCAGCGTTAAGCTCAACTTCAGCTCGCAAGGCCAACAGGTCAAAAGCAGAACCGAACTTAGGATGCt                                                                                 >  2:609224/1‑150 (MQ=255)

                                                                      tCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTcc                                                                               >  1:1347109/1‑150 (MQ=255)

                                                                      tCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTcc                                                                               <  2:736061/150‑1 (MQ=255)

                                                                       ccAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCCACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCCCGGAACTTAGGATGCTccc                                                                              >  1:1039205/1‑149 (MQ=255)

                                                                       ccAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCca                                                                              <  2:435460/150‑1 (MQ=255)

                                                                       ccAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCca                                                                              <  2:1508788/150‑1 (MQ=255)

                                                                        cAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCCCCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCCACTTCCGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCcag                                                                             >  1:1191728/1‑150 (MQ=255)

                                                                        cAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGCACTTAGGATGCTCcag                                                                             >  1:1480753/1‑150 (MQ=255)

                                                                        cAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGc                                                                                  <  1:731275/145‑1 (MQ=255)

                                                                        cAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGc                                                                                  >  2:731275/1‑145 (MQ=255)

                                                                        cAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCcag                                                                             >  2:1469566/1‑150 (MQ=255)

                                                                        cAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATATCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACATGTCATAAGCCGCACGGAACTTAGGATGCTCcag                                                                             >  1:1616170/1‑150 (MQ=255)

                                                                         aGCTCGTTGAGCATCGCTTTTTGGTCTGGTGTCGCGGAAACCTGGACCTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCGCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCcagc                                                                            <  2:86154/150‑1 (MQ=255)

                                                                         aGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGATCTCACCCCACCATTTCACCAGACCCTGCAGTTCAGCGTTACGCTCAACTTCCGCTCGCAAGGCCACCAGGTCATAAGCCGCACGGAACTTAGGATGCTCcagc                                                                            >  2:470404/1‑150 (MQ=255)

                                                                                       gccTTTTTGGTCTGGTGGCGCGGATACCTTGTAGTGACCGCACCATTTCACGAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGTTGATCCAGCAGTTTCAATGcgcg                                                              <  2:931281/149‑1 (MQ=255)

                                                                                          tttttGGTCTGGTGGAGCGGAAACCTGGAACTCAACCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCTACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGttt                                                           <  2:102105/150‑1 (MQ=255)

                                                                                              tGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTAccc                                                       <  2:1347109/150‑1 (MQ=255)

                                                                                              tGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTAccc                                                       <  2:1473156/150‑1 (MQ=255)

                                                                                              tGGTCAGGTGGCGCGGAAACCTGGAACTCACCCCACCAATTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTACCATGCGCGTTAAccc                                                       >  2:234644/1‑150 (MQ=255)

                                                                                               ggTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCt                                                      <  2:1540037/150‑1 (MQ=255)

                                                                                               ggTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCt                                                      <  1:580360/150‑1 (MQ=255)

                                                                                               ggTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCt                                                      >  1:1574762/1‑150 (MQ=255)

                                                                                                 tCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTga                                                    <  1:1475340/150‑1 (MQ=255)

                                                                                                           cgGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTTCCCTGACGACGGGGCa                                          >  2:1284931/1‑150 (MQ=255)

                                                                                                           cgGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACa                                          <  1:100719/150‑1 (MQ=255)

                                                                                                             gAAACCTGGAACTCACCCCACCCTTTCACACGACGCTGCAGGTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATa                                        <  2:340371/150‑1 (MQ=255)

                                                                                                             gAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATa                                        <  1:744115/150‑1 (MQ=255)

                                                                                                               aaCCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACg                                      <  2:1449230/150‑1 (MQ=255)

                                                                                                                 ccTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgca                                    <  2:638660/150‑1 (MQ=255)

                                                                                                                 ccTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgca                                    >  2:776175/1‑150 (MQ=255)

                                                                                                                   tGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgcaac                                  <  1:395221/150‑1 (MQ=255)

                                                                                                                       aCTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgcaactgca                              <  1:18114/150‑1 (MQ=255)

                                                                                                                       aCTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgcaactgca                              <  2:691278/150‑1 (MQ=255)

                                                                                                                          cACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgcaactgcaact                           >  1:116664/1‑150 (MQ=255)

                                                                                                                               caccaTTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAg                      >  2:1188029/1‑150 (MQ=255)

                                                                                                                                 ccaTTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGat                    <  2:1616170/150‑1 (MQ=255)

                                                                                                                                  caTTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGata                   >  2:1132535/1‑150 (MQ=255)

                                                                                                                                     ttCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATcg                <  1:1469566/150‑1 (MQ=255)

                                                                                                                                         cgtaaCGCTGAACTGCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCggg            >  1:989585/5‑150 (MQ=255)

                                                                                                                                          cAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCGCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCACTTTCCCTGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGt           >  1:1222341/1‑150 (MQ=255)

                                                                                                                                            gACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTa         <  1:1469905/150‑1 (MQ=255)

                                                                                                                                            cgagCTGAAGTTGAGCGTAACGCTGAACTGCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgcaact                                 >  2:1061937/4‑126 (MQ=38)

                                                                                                                                            cgagCTGAAGTTGAGCGTAACGCTGAACTGCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACgcaact                                 <  1:1061937/123‑1 (MQ=38)

                                                                                                                                             aCGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTaa        <  1:1534192/150‑1 (MQ=255)

                                                                                                                                              cGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTAAt       >  2:1309060/1‑150 (MQ=255)

                                                                                                                                              agCTGAAGTTGAGCGTAACGCTGAACTGCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTAAt       >  1:665974/2‑150 (MQ=255)

                                                                                                                                               gCTGAAGTTGAGCGTAACGCTGAACTGCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTGAGGATGCTCCAGCAGTATCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGCTATCGCGGGTTAATg      >  1:263374/1‑150 (MQ=255)

                                                                                                                                                cTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTAATgt     <  1:1613418/150‑1 (MQ=255)

                                                                                                                                                cTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTAATgt     <  1:1132535/150‑1 (MQ=255)

                                                                                                                                                   cAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTAATgtcgt  >  2:197056/1‑150 (MQ=255)

TGCGGTACCCTCACGACGTGGTGCGCGTTTGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCTGGTGGCGCGGAAACCTGGAACTCACCCCACCATTTCACCAGACGCTGCAGTTCAGCGTTACGCTCAACTTCAGCTCGCAAGGCCAACAGGTCATAAGCCGCACGGAACTTAGGATGCTCCAGCAGTTTCCATGCGCGTTTACCCTGACGACGGGACATACGCAACTGCAACTGCCAGATATCGCGGGTTAATGTCGT  >  NC_000913/157732‑158028

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