BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,406,093	A→G	6.7%	I39V (ATC→GTC)	smrA →	DNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,406,093	0	A	G	6.7%	345.6 / 11.9	119	I39V (ATC→GTC)	smrA	DNA endonuclease
Reads supporting (aligned to +/- strand): ref base A (51/59); new base G (4/4); total (56/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 120 total aligned reads displayed.

ATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTG  >  NC_000913/1405947‑1406242

 tatCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGa                                                                                                                                                   >  2:32886/1‑150 (MQ=255)

  atCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                                                                  >  1:1854169/1‑150 (MQ=255)

  atCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCAGCTGATACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGac                                                                                                                                                  >  2:735277/1‑150 (MQ=255)

    cgcgCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGacac                                                                                                                                                <  1:1063508/150‑1 (MQ=255)

      cgcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTAGACGCCAAGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTCTCGACACGc                                                                                                                                              >  2:1861248/1‑150 (MQ=255)

       gcTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTACTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGGGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCt                                                                                                                                             <  2:333318/150‑1 (MQ=255)

          ctTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCa                                                                                                                                          >  1:1323373/1‑150 (MQ=255)

           ttAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                                                                         <  1:530202/150‑3 (MQ=255)

           ttAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCAGCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGcgg                                                                                                                                         <  1:1275822/150‑3 (MQ=255)

            tacttGCGGGAGGAGGTACCATGACACTTGACGACAAATCGCTGTTTCTAGACGCCATGGACGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                                                                        <  1:353812/147‑1 (MQ=255)

            tAATTGCGGGAGGATGTAACATGAACCTTGACGACAATTCGCTTTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTGTCGATACGCTGCGGc                                                                                                                                        <  1:1175217/150‑1 (MQ=255)

            tAATTGCGGGAGGAGGTAACATGAACGTTGACGACAAATCGCTGTTTCTTGACGCCATGGCAGATGTCCAGCCGCTGAAACGTGCTACCGATGACCACTGGCTTCCATCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGc                                                                                                                                        <  1:1281226/150‑1 (MQ=255)

             aaTAGCGGGAGGATGTAACATGAACCTAGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCt                                                                                                                                       <  1:1562078/150‑1 (MQ=255)

               ttGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTg                                                                                                                                     <  1:1829538/150‑1 (MQ=255)

                tGCTGGAGGATGTAACATGAAGCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCATCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                                                                    <  1:502716/150‑1 (MQ=255)

                tGCGGGAGGATGTAACATGAACGTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGGGCTACCGATGTCCACTGGCATCCAACGCGTAACCGACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGa                                                                                                                                    <  2:190921/150‑1 (MQ=255)

                    ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAt                                                                                                                                >  2:1691237/1‑150 (MQ=255)

                    ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAAt                                                                                                                                >  1:1416217/1‑150 (MQ=255)

                    ggaggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACACCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGAAGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGCCACGCGGCAGCTTGATAAt                                                                                                                                >  1:221826/1‑150 (MQ=255)

                       ggaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTc                                                                                                                             <  1:1624341/150‑1 (MQ=255)

                        gaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCCGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                            >  1:158133/1‑150 (MQ=255)

                        gaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                            <  2:1699018/150‑1 (MQ=255)

                        gaTGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGAGGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTcc                                                                                                                            <  2:907562/150‑1 (MQ=255)

                               aaTGAACCTTGAAGTCAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCAACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGTCACGCTGCAGCTTGATAATTTCCTcacctc                                                                                                                     <  2:45078/149‑3 (MQ=255)

                                    aCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAt                                                                                                                >  1:885958/1‑150 (MQ=255)

                                    aCCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGATGCAGCTTGATAATTTCCTCACAACCGGAt                                                                                                                >  2:1320166/1‑150 (MQ=255)

                                     ccTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGAtt                                                                                                               >  1:658277/1‑150 (MQ=255)

                                        tGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTtct                                                                                                            <  1:808414/150‑1 (MQ=255)

                                         gacgacAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTATCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCGTGATAATTTCCTCACCACCGGATTtctc                                                                                                           >  2:1518991/1‑150 (MQ=255)

                                          acgacAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCg                                                                                                          >  2:1753685/1‑150 (MQ=255)

                                           cgacAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGCCACGCTGCAGCTTGATAATTTCCTCACCACCGGTTTTCTCGa                                                                                                         >  2:562891/1‑150 (MQ=255)

                                              cAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcat                                                                                                      <  2:906608/150‑1 (MQ=255)

                                              cAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTAGATAATTTCCTCACCACCGGATTTCTCGAcat                                                                                                      >  1:1344931/1‑150 (MQ=255)

                                               aaaTCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatc                                                                                                     <  2:1797869/150‑1 (MQ=255)

                                                aaTCGCTGTTTCTTGACGCCTTGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCCCCACCGGATTTCTCGAcatca                                                                                                    >  2:1422393/1‑150 (MQ=255)

                                                 aTCGCTGTTTCTTGACGCCATGGAAGATGACCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGAcatcat                                                                                                   <  2:1112604/150‑1 (MQ=255)

                                                    gCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATccc                                                                                                <  2:324920/150‑1 (MQ=255)

                                                      tGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                              <  1:1651728/150‑1 (MQ=255)

                                                      tGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCAc                                                                                              <  2:1191490/150‑1 (MQ=255)

                                                        tttCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTa                                                                                            <  1:424101/150‑1 (MQ=255)

                                                            ttGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTc                                                                                        >  2:1137610/1‑150 (MQ=255)

                                                             tGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                                                                       <  2:1598891/150‑1 (MQ=255)

                                                             tGACGCCATGGAAGATGTACAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCa                                                                                       <  2:1591343/150‑1 (MQ=255)

                                                              gACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAg                                                                                      <  2:1058062/150‑1 (MQ=255)

                                                                cGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGcc                                                                                    >  1:799380/1‑150 (MQ=255)

                                                                  ccATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGc                                                                                  >  1:618727/1‑150 (MQ=255)

                                                                   cATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCt                                                                                 >  1:1414066/1‑150 (MQ=255)

                                                                       gAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCAACGTATTTCTCGCCATCATCCCACTAAGTCAGCAGCTGGAg                                                                             >  1:917525/1‑150 (MQ=255)

                                                                         aGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGtt                                                                           <  2:1854169/150‑1 (MQ=255)

                                                                            tGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTcg                                                                        <  2:1399686/150‑1 (MQ=255)

                                                                                cAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTACCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCCGCCGCTGGAGTTTcggcgg                                                                    >  2:546307/1‑150 (MQ=255)

                                                                                 aGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCggg                                                                   <  1:1215880/150‑1 (MQ=255)

                                                                                  gCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGa                                                                  <  2:94432/150‑1 (MQ=255)

                                                                                     gcggataCGTGCTACCGATGTGCACTGGCATCCAACGCGTAACCAACCTGCGCCGCAGCGTATCGACCCGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                                                               <  2:126213/144‑1 (MQ=255)

                                                                                     gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                                                               >  2:95002/1‑150 (MQ=255)

                                                                                     gCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAgg                                                               >  1:1465430/1‑150 (MQ=255)

                                                                                       tGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGt                                                             >  1:1800648/1‑150 (MQ=255)

                                                                                        gAAACGTGCTACCGATGTCCACTGGCATCCACCGCGTAACCAACGTGCGCCGCAGCGTATCCCCCCCCTGCAGCTTGCTAATTTCCTCACCCCCGGAGGTCGGGACATCTTCCCACTTCGTCAGACGCGGGAGTTTCGGCAGGAAGGGtt                                                            >  1:302101/1‑150 (MQ=255)

                                                                                        gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt                                                            >  2:1056908/1‑150 (MQ=255)

                                                                                        gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGtt                                                            <  2:1323373/150‑1 (MQ=255)

                                                                                        gAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTGCGGCGGGAAGGGtt                                                            >  2:135933/1‑150 (MQ=255)

                                                                                          aaCGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGc                                                          <  1:1166981/150‑1 (MQ=255)

                                                                                                 tACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgg                                                   <  2:409384/150‑1 (MQ=255)

                                                                                                  aCCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATggg                                                  >  1:1195288/1‑150 (MQ=255)

                                                                                                   ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg                                                 >  1:1214211/1‑150 (MQ=255)

                                                                                                   ccGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATgggg                                                 <  2:1651662/150‑1 (MQ=255)

                                                                                                    cGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGt                                                <  2:1344931/150‑1 (MQ=255)

                                                                                                     gATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACCTCATCCCACTAAGTAAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACAGgggggg                                               >  2:487248/1‑148 (MQ=255)

                                                                                                       tGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTATGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCt                                             <  1:238441/150‑1 (MQ=255)

                                                                                                          ccACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGa                                          >  2:450051/1‑150 (MQ=255)

                                                                                                           cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                         <  2:618635/150‑1 (MQ=255)

                                                                                                           cACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGAt                                         >  2:1163752/1‑150 (MQ=255)

                                                                                                            aCTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGTGTGCTGGATa                                        <  1:735277/150‑1 (MQ=255)

                                                                                                               ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                     >  2:390952/1‑150 (MQ=255)

                                                                                                               ggCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGc                                     <  1:257605/150‑1 (MQ=255)

                                                                                                                    ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                >  1:1719375/1‑150 (MQ=255)

                                                                                                                    ccAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTgcgc                                <  2:1017634/150‑1 (MQ=255)

                                                                                                                      aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAg                              >  2:1520995/1‑150 (MQ=255)

                                                                                                                      aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAg                              <  2:1171926/150‑1 (MQ=255)

                                                                                                                      aaCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAg                              <  1:1398316/150‑1 (MQ=255)

                                                                                                                       aCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGt                             >  1:587261/1‑150 (MQ=255)

                                                                                                                       aCGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGGTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGt                             <  1:1247255/150‑1 (MQ=255)

                                                                                                                         gcgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTgg                           <  2:1404837/150‑1 (MQ=255)

                                                                                                                          cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                          <  2:522155/150‑1 (MQ=255)

                                                                                                                          cgTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGt                          >  2:692417/1‑150 (MQ=255)

                                                                                                                           gTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTa                         >  2:1321876/1‑150 (MQ=255)

                                                                                                                            tAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTaa                        <  1:1605591/150‑1 (MQ=255)

                                                                                                                              aCCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATATCAACACCACCGGATTTCTCGACATCCTCCCACTAAGACAGCCGCTGGAGTTTCGGAGGCAAGGGATCCAACATGGGGTGAAGGATAAGCTGCGCAGTggaatat                      >  1:238241/1‑145 (MQ=255)

                                                                                                                                  aCGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATcc                  <  1:1188620/150‑1 (MQ=255)

                                                                                                                                   cGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCg                 <  1:896857/150‑1 (MQ=255)

                                                                                                                                    gTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGc                <  1:1691237/150‑1 (MQ=255)

                                                                                                                                      gcgcCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGCCATCATCCCCCTAAGTCAGCCTCTGGCGTTTCGGCGGTAATGGTGGGACCATGGGGTGCTGGATAAGCAGCGCAGTGCGAAATAGCCGCaa              >  1:382492/1‑150 (MQ=255)

                                                                                                                                       agctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTACGGCGGGAAGGGTTGCAAGATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAAc             >  2:276527/5‑150 (MQ=255)

                                                                                                                                         ctgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAg           >  1:1754934/3‑150 (MQ=255)

                                                                                                                                          tgcAGCGTGTCGATACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAAAATGGGGTGGTGGATAAGCTGCGCTGTGGTAAATATCCGCAACAgg          >  2:920567/2‑150 (MQ=255)

                                                                                                                                            cAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCg        <  1:1058903/150‑1 (MQ=255)

                                                                                                                                              gCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAg      <  1:1137610/150‑1 (MQ=255)

                                                                                                                                                 tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCt   <  1:1875670/150‑1 (MQ=255)

                                                                                                                                                 tATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAACTCAGCCGCTGGAGGTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCt   <  1:1154378/150‑1 (MQ=255)

                                                                                                                                                  aTCGACACGCTGCAGCTTGATAATTTCCTCACCACCCGATTTCTCGACATCATCCCACTATGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCTACAGGCGAGCCTg  <  2:1233449/150‑1 (MQ=255)

ATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTGTTTCTTGACGCCATGGAAGATGTCCAGCCGCTGAAACGTGCTACCGATGTCCACTGGCATCCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTG  >  NC_000913/1405947‑1406242

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