BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,215	T→A	8.4%	N123Y (AAC→TAC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,215	0	T	A	8.4%	278.9 / 13.6	107	N123Y (AAC→TAC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base T (49/49); new base A (4/5); total (53/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 107 total aligned reads displayed.

CAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTA  >  NC_000913/3566067‑3566362

cagTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                                                                                                                                    >  2:1658682/1‑150 (MQ=255)

cagTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                                                                                                                                    >  2:408981/1‑150 (MQ=255)

  gtagtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAcc                                                                                                                                                  <  2:1669520/150‑1 (MQ=255)

    agtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCGAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGa                                                                                                                                                >  2:1198360/1‑150 (MQ=255)

    agtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGATTCAAGGAAGCAAGCCGCCAGACGTCCCCGGCCACCGTTACCgc                                                                                                                                                >  1:102511/1‑149 (MQ=255)

     gtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                                                                                                                               <  1:614001/150‑1 (MQ=255)

     gtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                                                                                                                               >  2:1325676/1‑150 (MQ=255)

      tCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                                                                                                                              <  1:13902/150‑1 (MQ=255)

         ggCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCt                                                                                                                                           >  1:1636943/1‑150 (MQ=255)

           cGACTCTATCTGGCTACCGTTAACGATGATCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTgg                                                                                                                                         <  2:856076/150‑1 (MQ=255)

             aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATAGCGTATTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCTAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                       <  1:1279048/150‑2 (MQ=255)

             aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                       <  1:1223626/150‑2 (MQ=255)

             aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                       <  1:1709580/150‑2 (MQ=255)

             aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                       <  1:1703084/150‑2 (MQ=255)

             aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGt                                                                                                                                       <  2:528781/150‑1 (MQ=255)

              ctctTTCTGGCTACCGTTTACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGCGGCCTGGGTc                                                                                                                                      >  2:560515/1‑150 (MQ=255)

                ctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAt                                                                                                                                    >  1:653089/1‑150 (MQ=255)

                     tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttct                                                                                                                               <  1:471182/150‑1 (MQ=255)

                     tGGCTACCGTTAACGATGGTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttct                                                                                                                               <  1:698891/150‑1 (MQ=255)

                       gCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                                             >  2:1107968/1‑150 (MQ=255)

                       gCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                                             <  2:1594404/150‑1 (MQ=255)

                       gCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                                             <  1:548793/150‑1 (MQ=255)

                        cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTg                                                                                                                                          <  1:1324496/136‑1 (MQ=255)

                        cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTg                                                                                                                                          >  2:1324496/1‑136 (MQ=255)

                        cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                                            >  2:546367/1‑150 (MQ=255)

                          aCCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcatc                                                                                                                          <  2:1273036/150‑1 (MQ=255)

                           ccGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGTTACCGAGGGCTGGGTCATTTTAGtcatca                                                                                                                         >  2:1523504/1‑150 (MQ=255)

                           ccGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcatca                                                                                                                         >  1:126056/1‑150 (MQ=255)

                             gTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAACGAGGCCTGGGTCATTttcttaaaaaat                                                                                                                       >  2:1499144/1‑143 (MQ=255)

                              ttAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaatc                                                                                                                      <  1:888766/150‑1 (MQ=255)

                               tAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaatca                                                                                                                     <  1:546367/150‑1 (MQ=255)

                                 aCGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGc                                                                                                                   <  2:213204/150‑1 (MQ=255)

                                     tGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCtt                                                                                                               <  1:1325676/150‑1 (MQ=255)

                                       ttCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCg                                                                                                             <  2:935920/150‑1 (MQ=255)

                                         cTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGTCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCGGGGTCATTTTCTTCGTCGATCAGCTCTTCgag                                                                                                           >  1:664106/1‑150 (MQ=255)

                                          tGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgaga                                                                                                          <  1:105542/150‑1 (MQ=255)

                                           gCTTGAACATACCGTAGTCATAGCGGATGCCGAAACCGCGCCCCGGCAACCCTAACTTCGCCAGAGATTCAAGGAAGCAAGCCGCCAGACGTCCCAGCCCACCCTTCCCGAGGCCTGGGTCCATTTCTTCATCAAACCGCTCTTCgccaa                                                                                                         >  1:1446398/1‑146 (MQ=255)

                                              tGAACATACCGTAGTCATAGCGGAGGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTa                                                                                                      <  1:702440/150‑1 (MQ=255)

                                                aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGAAATCCGCCAGACGTCCGAGGCCTCGGGTACCGAGGGCTGGGTGATTTTATTCATAAATCAGGTCTTCGAGATTTTAc                                                                                                    >  2:1372042/1‑150 (MQ=255)

                                                 aCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAcc                                                                                                   <  2:1566810/150‑1 (MQ=255)

                                                 aCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAcc                                                                                                   <  1:1039068/150‑1 (MQ=255)

                                                  cATACCGTTGTCATTGCGGATGCCGTATCCGCGCCGCGGCATCCCTAACGTCGCAAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAccc                                                                                                  <  1:843290/150‑1 (MQ=255)

                                                  cATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAccc                                                                                                  >  2:567831/1‑150 (MQ=255)

                                                    tACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGCGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCa                                                                                                >  2:1059306/1‑150 (MQ=255)

                                                         tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCt                                                                                           <  1:650175/150‑1 (MQ=255)

                                                            tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTcc                                                                                        <  1:952350/150‑1 (MQ=255)

                                                             cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                                       <  1:1658682/150‑1 (MQ=255)

                                                             cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTATTCATCAATCAGCTCTTCGAGTTTTAACCCCATCGCTTCCa                                                                                       >  2:569220/1‑150 (MQ=255)

                                                               tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCCGGCCCCCGGTGCCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAGTTAACCCCCTCGCTTCCAgg                                                                                     >  1:599783/1‑149 (MQ=255)

                                                                aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGGCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTg                                                                                    <  2:126056/150‑1 (MQ=255)

                                                                aGCGGATGACGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTg                                                                                    <  1:1499144/150‑1 (MQ=255)

                                                                  cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTgcg                                                                                  <  2:592026/150‑1 (MQ=255)

                                                                              ccGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCt                                                                      >  1:659275/1‑150 (MQ=255)

                                                                                  gcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGt                                                                  >  2:422461/1‑150 (MQ=255)

                                                                                  gcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGt                                                                  >  2:428745/1‑150 (MQ=255)

                                                                                     ccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaaa                                                               >  2:202342/1‑150 (MQ=255)

                                                                                      cGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCACGCCACCCTTACCGCGGCCTGGGTCAATTAATACATCAATCAGCACTTCGAAAATTAACCCCATAGCGTCCACAGCCCCCTGTACATCTTCGTaaaa                                                              >  2:1081158/1‑149 (MQ=255)

                                                                                         cAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACTTCCCAGGCCACCGTTACCGAGGCCTGGGGCAGTTTCGTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCATCCAGTGCGCCCTGTACAACTTCGTAAATTcg                                                           >  2:1047661/1‑149 (MQ=255)

                                                                                            cccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg                                                        <  1:849016/150‑1 (MQ=255)

                                                                                            cccTAACGTCGCCAGAGAATCAAGGAAGAAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTACATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg                                                        >  2:1111186/1‑150 (MQ=255)

                                                                                                aaCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGac                                                    >  1:1218032/1‑150 (MQ=255)

                                                                                                aaCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGac                                                    <  1:1623750/150‑1 (MQ=255)

                                                                                                 aCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGaca                                                   <  1:941605/150‑1 (MQ=255)

                                                                                                  cGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTTCCGAGGCCTGGGTCATTTTCTTCATAAATCAGCTCGTCGAGATTTAACCCCATCGCTTACAGTGCGCCCTGTACCTCTTCGTAAATTCCTAGCGacca                                                  >  2:591480/1‑148 (MQ=255)

                                                                                                   gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaac                                                 <  2:450302/150‑1 (MQ=255)

                                                                                                   gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaac                                                 <  1:312417/150‑1 (MQ=255)

                                                                                                      gCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATg                                              <  2:653089/150‑1 (MQ=255)

                                                                                                       ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATgg                                             <  1:1143474/150‑1 (MQ=255)

                                                                                                        cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGATGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGAGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGc                                            <  2:1335343/150‑1 (MQ=255)

                                                                                                          gagaATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGt                                          <  1:1059306/150‑1 (MQ=255)

                                                                                                             aaTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTgg                                       >  1:188932/1‑150 (MQ=255)

                                                                                                                cAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCCTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgaga                                    >  2:1658828/1‑150 (MQ=255)

                                                                                                                  aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGc                                  <  2:188932/150‑1 (MQ=255)

                                                                                                                   ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGTGAGCg                                 >  1:62462/1‑150 (MQ=255)

                                                                                                                    gAAGCAAGCAGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCATCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGt                                <  1:258281/150‑1 (MQ=255)

                                                                                                                       gcaagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTATTCAAAAATCACCTCTTCGAGATTTAACCCCAACGCTTCCAGTGCGCCCAGTACATCTTCGTAAATTACTAGCGACAACATCGCGTTGGAGAGCGTACg                             >  2:1686778/1‑150 (MQ=255)

                                                                                                                         aagcCGCCAGACGTCCCAGGCCACCGTTAACGACGCCTGGGTCATTTTGTTCAACAATAAGCTATTAGAGATTTTAACCCATAGCTTACACTGAGCCCGGTTAAAATTAGTAAAATTATAGCGACAACAAGGAGGTTGGGAGCGTaaagc                           >  2:1422052/1‑146 (MQ=255)

                                                                                                                           gccgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAGTTCCCCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAGTCCTAGCGACAACATGGCGGTGGCGAGCGTACGGCCa                         >  1:17985/1‑150 (MQ=255)

                                                                                                                             cgccAGTCGTCCCAGGCCACCGTTACCGAGTCCGGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATAGCTTCCAGTGCGCCCTGTACATGTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAAt                       <  1:1081158/150‑1 (MQ=255)

                                                                                                                             cgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCATCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAAt                       >  2:728484/1‑150 (MQ=255)

                                                                                                                                aaaatgacccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaa                    >  2:783651/8‑150 (MQ=255)

                                                                                                                                  gACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaa                  >  2:1293841/1‑150 (MQ=255)

                                                                                                                                  gACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaa                  <  1:1210660/150‑1 (MQ=255)

                                                                                                                                    cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTTGAGAGCGTACGGCCAATCAAAAAc                <  1:1686778/150‑1 (MQ=255)

                                                                                                                                     gTCCCAGGCCACCGTTACCGAGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACt               >  2:525677/1‑150 (MQ=255)

                                                                                                                                      acccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTc              >  2:375299/2‑150 (MQ=255)

                                                                                                                                       cccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTcc             >  2:1417567/1‑150 (MQ=255)

                                                                                                                                        ccAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCa            <  1:567831/150‑1 (MQ=255)

                                                                                                                                        ccAGGCCACCGTTACCGAGGCCTGGGTCATTTTCATCAACAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAATCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCa            <  1:1372042/150‑1 (MQ=255)

                                                                                                                                         cAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt           >  1:1088976/1‑150 (MQ=255)

                                                                                                                                         cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCCTCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGAATTAGGAGAGCGTACGGCCAATCAAAACCTCCag           >  2:282727/1‑149 (MQ=255)

                                                                                                                                         cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt           >  2:1209644/1‑150 (MQ=255)

                                                                                                                                         cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt           >  1:904067/1‑150 (MQ=255)

                                                                                                                                          cggCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCGCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATc          <  2:1303303/149‑1 (MQ=255)

                                                                                                                                           ggCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCg         <  1:560515/150‑1 (MQ=255)

                                                                                                                                             ccACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAc       >  2:1392764/1‑150 (MQ=255)

                                                                                                                                              cACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACa      >  1:45877/1‑150 (MQ=255)

                                                                                                                                                ccGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAgg    >  1:608953/1‑150 (MQ=255)

                                                                                                                                                  gTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGgta  >  1:1442952/1‑150 (MQ=255)

                                                                                                                                                  gTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGgta  <  2:904067/150‑1 (MQ=255)

CAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTA  >  NC_000913/3566067‑3566362

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