BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,758,946	G→T	8.2%	L115L (CTC→CTA)	sufE ←	sulfur acceptor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,758,946	0	G	T	8.2%	623.8 / 32.2	170	L115L (CTC→CTA)	sufE	sulfur acceptor protein
Reads supporting (aligned to +/- strand): ref base G (75/81); new base T (7/7); total (82/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 170 total aligned reads displayed.

TGTGCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAA  >  NC_000913/1758805‑1759094

tgtgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTg                                                                                                                                              <  2:1854189/150‑1 (MQ=255)

tgtgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTg                                                                                                                                              >  1:1225670/1‑150 (MQ=255)

 gtgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTgg                                                                                                                                             <  1:262841/150‑1 (MQ=255)

 gtgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTgg                                                                                                                                             <  1:469884/150‑1 (MQ=255)

 gtgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTgg                                                                                                                                             <  1:469894/150‑1 (MQ=255)

  tgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTggg                                                                                                                                            <  2:2146432/150‑1 (MQ=255)

  tgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTggg                                                                                                                                            <  2:1016560/150‑1 (MQ=255)

  tgCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTggg                                                                                                                                            <  2:278118/150‑1 (MQ=255)

        tCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTCAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGAGGTTGGGTGAgcg                                                                                                                                      >  1:558283/1‑150 (MQ=255)

        tCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAgcg                                                                                                                                      >  1:1464840/1‑150 (MQ=255)

               aaaCATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTCTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCAtttt                                                                                                                               >  2:1733777/1‑150 (MQ=255)

                aaCATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCAttttt                                                                                                                              >  1:74508/1‑150 (MQ=255)

                 aCATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCAtttttt                                                                                                                             <  2:1516889/150‑1 (MQ=255)

                   aTTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCa                                                                                                                           >  2:1255598/1‑150 (MQ=255)

                    ttGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCaa                                                                                                                          <  2:1321580/150‑1 (MQ=255)

                      gACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAAc                                                                                                                        >  1:1382397/1‑150 (MQ=255)

                      gACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAAc                                                                                                                        >  1:1019780/1‑150 (MQ=255)

                      gACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAAc                                                                                                                        <  1:1581911/150‑1 (MQ=255)

                      gACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAAc                                                                                                                        >  1:295749/1‑150 (MQ=255)

                        ccATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCa                                                                                                                      <  2:965025/150‑1 (MQ=255)

                         cATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCAc                                                                                                                     <  1:1817750/150‑1 (MQ=255)

                          aTCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTTAGTGCATCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGTACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACg                                                                                                                    <  2:872947/150‑1 (MQ=255)

                          aTCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACg                                                                                                                    <  2:873373/150‑1 (MQ=255)

                           tCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACgg                                                                                                                   >  1:2301688/1‑150 (MQ=255)

                             cGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGAc                                                                                                                 >  1:2141275/1‑150 (MQ=255)

                                    gtgAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATc                                                                                                          <  1:2319868/150‑1 (MQ=255)

                                       aGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGaa                                                                                                       >  1:33895/1‑150 (MQ=255)

                                        gCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGaaa                                                                                                      >  2:727414/1‑150 (MQ=255)

                                         ccAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAAt                                                                                                     >  2:1363093/1‑150 (MQ=255)

                                         ccAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAAt                                                                                                     >  1:177830/1‑150 (MQ=255)

                                         ccAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAAt                                                                                                     <  1:2253453/150‑1 (MQ=255)

                                          cAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAAtt                                                                                                    >  2:1739529/1‑150 (MQ=255)

                                           aaCCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTg                                                                                                   <  1:2007464/150‑1 (MQ=255)

                                           aaCCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTg                                                                                                   <  1:113902/150‑1 (MQ=255)

                                              cGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACa                                                                                                >  2:1368248/1‑150 (MQ=255)

                                              cGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACa                                                                                                >  2:458404/1‑150 (MQ=255)

                                              cGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACa                                                                                                >  2:1210477/1‑150 (MQ=255)

                                                     aaGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCt                                                                                         <  1:199710/150‑1 (MQ=255)

                                                      aGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTg                                                                                        >  2:1947699/1‑150 (MQ=255)

                                                           tCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTgcggcg                                                                                   >  1:1998738/1‑150 (MQ=255)

                                                           tCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTgcggcg                                                                                   >  1:1198149/1‑150 (MQ=255)

                                                           tCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTgcggcg                                                                                   >  1:1587969/1‑150 (MQ=255)

                                                           tCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTgcggcg                                                                                   >  1:1662650/1‑150 (MQ=255)

                                                            ccTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGt                                                                                  <  2:1594186/150‑1 (MQ=255)

                                                               tttagtttagCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCAt                                                                               >  2:1586521/1‑150 (MQ=255)

                                                               tttagtttagCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCAt                                                                               <  1:819910/150‑1 (MQ=255)

                                                                ttagtttagCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATc                                                                              <  1:934356/150‑1 (MQ=255)

                                                                   gtttagCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGa                                                                           >  1:1305194/1‑150 (MQ=255)

                                                                     ttagCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATc                                                                         <  1:74581/150‑1 (MQ=255)

                                                                     ttagCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATc                                                                         <  1:1501629/150‑1 (MQ=255)

                                                                     ttagCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATc                                                                         <  1:1607822/150‑1 (MQ=255)

                                                                       agCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGt                                                                       >  1:703189/1‑150 (MQ=255)

                                                                       agCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGt                                                                       >  1:2227112/1‑150 (MQ=255)

                                                                        gCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTa                                                                      <  1:382925/150‑1 (MQ=255)

                                                                        gCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTa                                                                      <  1:382952/150‑1 (MQ=255)

                                                                         cTAAGTGCAGCGTCTTTGGCGCTAATTGCGCTAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCTAAATTGACAATATCCTGCGGCGTCATCTGATCGTag                                                                     <  2:1492440/150‑1 (MQ=255)

                                                                         cTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTag                                                                     >  1:1777257/1‑150 (MQ=255)

                                                                         cTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTag                                                                     <  1:1310650/150‑1 (MQ=255)

                                                                          tAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTaga                                                                    <  1:1733479/150‑1 (MQ=255)

                                                                           aaGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTagag                                                                   >  2:2024637/1‑150 (MQ=255)

                                                                                 aGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCATGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATaaa                                                             <  1:1329192/150‑1 (MQ=255)

                                                                                 aGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATaaa                                                             >  2:2256764/1‑150 (MQ=255)

                                                                                 aGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATaaa                                                             <  1:628808/150‑1 (MQ=255)

                                                                                 aGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATaaa                                                             <  1:161129/150‑1 (MQ=255)

                                                                                  gCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAg                                                            <  2:2301688/150‑1 (MQ=255)

                                                                                  gCGGCTTTGGCGCGAAGTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAg                                                            <  2:168680/150‑1 (MQ=255)

                                                                                       tttGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAgacgac                                                       >  2:1730572/1‑150 (MQ=255)

                                                                                           gcgcGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCa                                                   >  2:922140/1‑150 (MQ=255)

                                                                                                      gcgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCttt                                        >  2:1516144/1‑150 (MQ=255)

                                                                                                      gcgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCttt                                        >  2:1797637/1‑150 (MQ=255)

                                                                                                      gcgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCttt                                        >  2:395557/1‑150 (MQ=255)

                                                                                                      gcgAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCttt                                        <  2:1653516/150‑1 (MQ=255)

                                                                                                          atcatcGCTTCCACACCTTGAGAACGAGATGGGGTGAGATGTTGGGTGAGCGCAAATTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTcaca                                    <  2:1139104/150‑1 (MQ=255)

                                                                                                              tcGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcg                                >  2:289395/1‑150 (MQ=255)

                                                                                                               cGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgc                               >  2:1648751/1‑150 (MQ=255)

                                                                                                               cGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgc                               >  2:2077391/1‑150 (MQ=255)

                                                                                                               cGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgc                               >  2:513890/1‑150 (MQ=255)

                                                                                                                gCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATcgcc                              >  1:98907/1‑150 (MQ=255)

                                                                                                                        aCCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCt                      >  2:563824/1‑150 (MQ=255)

                                                                                                                          cTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGt                    <  2:2385900/150‑1 (MQ=255)

                                                                                                                           ttGTGAACGAGTTGGGGTGAGATGTTGGGTGAGCGCAATTTTTACAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCTACATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTc                   <  1:618797/150‑1 (MQ=255)

                                                                                                                           ttGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTc                   >  1:998712/1‑150 (MQ=255)

                                                                                                                               taaCGAGATGTGTTGATATGTTTGGTGAGATAAATTTTTTCAAACCACGGACGAACATCGAAATTGAAAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGccc               <  2:1115231/149‑1 (MQ=255)

                                                                                                                               gAACGAGATGGGGTGAGATGTTGGGTGAGCCCAATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGccc               <  2:1115230/150‑1 (MQ=255)

                                                                                                                                aaCGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCt              <  1:1368248/150‑1 (MQ=255)

                                                                                                                                aaCGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCt              <  1:458404/150‑1 (MQ=255)

                                                                                                                                aaCGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCt              >  1:1425134/1‑150 (MQ=255)

                                                                                                                                aaCGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCt              >  2:575066/1‑150 (MQ=255)

                                                                                                                                aaCGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCt              >  1:1733399/1‑150 (MQ=255)

                                                                                                                                         gggTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaa     <  2:1998738/150‑1 (MQ=255)

                                                                                                                                         gggTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaa     <  2:1587969/150‑1 (MQ=255)

                                                                                                                                         gggTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaa     <  2:1587963/150‑1 (MQ=255)

                                                                                                                                         gggTTAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaa     <  2:465225/150‑1 (MQ=255)

                                                                                                                                          ggTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaat    >  2:424154/1‑150 (MQ=255)

                                                                                                                                           gTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaata   <  1:176635/150‑1 (MQ=255)

                                                                                                                                           gTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaata   >  2:1833693/1‑150 (MQ=255)

                                                                                                                                           gTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaata   <  2:1305194/150‑1 (MQ=255)

                                                                                                                                            tGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaataa  >  1:1699786/1‑150 (MQ=255)

                                                                                                                                            tGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCaataa  >  1:1080209/1‑150 (MQ=255)

                                                                                                                                            tGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGACGCCCTGTAATTCaataa  >  1:1263927/1‑150 (MQ=255)

TGTGCAGGTCAGTGTAAACATTGACCATCCGGCAATGTGAGCCAACCGGATGAAAGCTGTCCTTTTAGTTTAGCTAAGTGCAGCGGCTTTGGCGCGAATTGCGCGAATCATCGCTTCCAGACCTTGTGAACGAGATGGGGTGAGATGTTGGGTGAGCGCCATTTTTTCAAACCACGGACGCACATCGAAATTGACAATATCCTGCGGCGTCATCTGATCGTAGAGAATAAAGACGACCGCAATAAGCCCTTTCACAATCGCCGCATCGCTGTCGCCCTGTAATTCAATAA  >  NC_000913/1758805‑1759094

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