BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,585,156:1	+C	5.6%	coding (1331/2280 nt)	ydeP ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,585,156	1	.	C	5.6%	511.8 / 26.6	124	coding (1331/2280 nt)	ydeP	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (53/64); new base C (4/3); total (57/67)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 131 total aligned reads displayed.

AGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGG  >  NC_000913/1585009‑1585304

aGCTTAGTGGCTACGTGTACCGCCAAATCTAAATGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑ca                                                                                                                                                    <  1:907651/150‑1 (MQ=255)

 gCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cag                                                                                                                                                   <  1:1693544/150‑1 (MQ=255)

  cTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCATAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                                                                                  <  1:411812/150‑1 (MQ=255)

  cTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                                                                                  >  1:1693359/1‑150 (MQ=255)

   ttAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagca                                                                                                                                                 <  2:433565/150‑1 (MQ=255)

      gTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgc                                                                                                                                              <  2:1413730/150‑2 (MQ=255)

         gCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCa                                                                                                                                           <  1:219914/150‑1 (MQ=255)

             cgggtACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATCTTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTg                                                                                                                                       <  1:169151/147‑1 (MQ=255)

                gtACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCAt                                                                                                                                    >  2:902459/1‑150 (MQ=255)

                 tACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATg                                                                                                                                   <  1:1041921/150‑1 (MQ=255)

                    cGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                                >  2:26781/1‑150 (MQ=255)

                    cGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                                >  1:594889/1‑150 (MQ=255)

                     gCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTACAGGTGCATGAGGt                                                                                                                               >  2:889130/1‑150 (MQ=255)

                       cAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCACTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgg                                                                                                                             <  2:890123/150‑1 (MQ=255)

                        aaaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                                                                            <  1:608385/150‑1 (MQ=255)

                        aaaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                                                                            <  1:18178/150‑1 (MQ=255)

                        aaaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATAGGCAGCGCAAAATTGCCCGCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                                                                            <  2:516664/150‑1 (MQ=255)

                         aaTCTAATTGCGTTCACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                                                                           <  2:412299/150‑1 (MQ=255)

                         aaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGTCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                                                                           <  2:953337/150‑1 (MQ=255)

                         aaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                                                                           <  1:833948/150‑1 (MQ=255)

                           ccTAAATGCGTTAACGGTAACTCGCGCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTg                                                                                                                         <  2:1170004/149‑1 (MQ=255)

                              aaTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                                      <  1:1873474/150‑1 (MQ=255)

                              aaTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATAGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGATCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGTTGAAg                                                                                                                      >  2:731774/1‑150 (MQ=255)

                              aaTTGCGATAACGGTACAGCGCTCGCTTCCCGATCTGGCATAGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                                      <  2:266980/150‑1 (MQ=255)

                                   cGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                                 >  1:1713902/1‑150 (MQ=255)

                                   cGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                                 >  2:702801/1‑150 (MQ=255)

                                   cGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGGGGGGTGACGCCATa                                                                                                                 >  1:1153115/1‑150 (MQ=255)

                                    gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAg                                                                                                                <  2:1429662/150‑1 (MQ=255)

                                     ttAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGCTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTAGGAATTG‑CAGCCTGTCGAGGTGCATGAGGTGGGGTGAAGCCATAgc                                                                                                               >  1:469102/1‑150 (MQ=255)

                                        aCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCt                                                                                                            <  1:1325202/150‑1 (MQ=255)

                                         cGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTc                                                                                                           >  2:867643/1‑150 (MQ=255)

                                          ggTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCa                                                                                                          >  1:1033971/1‑150 (MQ=255)

                                           gTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCTTGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                                                                         >  2:1005023/1‑150 (MQ=255)

                                           gTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                                                                         <  1:437467/150‑1 (MQ=255)

                                           gTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                                                                         <  2:815955/150‑1 (MQ=255)

                                            tACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCAGGCAGATCAATGAACGACCCTGCCCCGTACATATTGATTCCATGCTGGAAATTG‑AAGCATGGCCAGTGGCATGAGGTGGGGTGCCGCCATACCGCTcccc                                                                                                        >  2:1187996/1‑150 (MQ=255)

                                              cAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCATAATTGCCCCCCATGCAGCTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcccc                                                                                                      >  1:298424/1‑149 (MQ=255)

                                               aGCGCTCGCTTCCCGATCTGGCATTCCCAGCGCAAAATACCCCCCCATCCAGATCAATGCTCGAGCCTGCCCCGTACCTATTGCTTGCATGCGGGCAATTG‑CAGCATGTCCCGGGGCATGAGGGGGGGTGAAGCCATAGCGCGCACCCAg                                                                                                     <  2:291989/150‑1 (MQ=255)

                                                gcgcTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGa                                                                                                    >  2:576349/1‑150 (MQ=255)

                                                gcgcTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCAGAGCGCTCACCCAGa                                                                                                    >  1:1798771/1‑150 (MQ=255)

                                                 cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCCTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGGGCATGAGGTGGGGGGAAGCCATAGCGCTCACCCAgcc                                                                                                   >  1:1125955/1‑148 (MQ=255)

                                                 cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGAc                                                                                                   <  1:1518899/150‑1 (MQ=255)

                                                 cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGAc                                                                                                   >  1:1518777/1‑150 (MQ=255)

                                                    gcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCATGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAg                                                                                                <  1:444839/149‑1 (MQ=255)

                                                      gctTCCCGATCTGGCAATGCCAGCGCAAATTTGCCCCCCATGCAGATCAATGCTCGAGCCGGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGTGTGAAGCCATAGCGCTCACCCAGACGAGcc                                                                                              <  1:1671865/150‑1 (MQ=255)

                                                       ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGCGCCTGCCCCGTAGTTATTGCTTGCATGCTGGCAATTG‑CAGCATGTACAGGTGCATGAGGTGGGGTGAAGCCTTAGCGCTCACCCAGACGAGCCa                                                                                             >  2:1662835/1‑150 (MQ=255)

                                                       ctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                             >  1:1507624/1‑150 (MQ=255)

                                                         tCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGCTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGa                                                                                           >  1:827969/1‑150 (MQ=255)

                                                          cccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaa                                                                                          <  1:51947/150‑1 (MQ=255)

                                                          cccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCAGGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaa                                                                                          >  1:1095699/1‑150 (MQ=255)

                                                            cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGATGCCATAGCGCTCACCCAGACGAGCCAGAAAc                                                                                        <  2:1451468/150‑1 (MQ=255)

                                                                cTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGGGCATGCGGTGGGGTGAAGCGATAGCGCTCCCCCAGCCGCGCCAGACACTCTg                                                                                    >  2:57586/1‑150 (MQ=255)

                                                                cTGGCATTGCCAGCGCAAAATTGACCCCCATGCATATCAATGCGCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTg                                                                                    <  1:1172374/150‑1 (MQ=255)

                                                                  ggCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa                                                                                  >  1:1307747/1‑150 (MQ=255)

                                                                  ggCATTGCCAGCGCAAAATTG‑CCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAg                                                                                 >  1:1619333/1‑150 (MQ=255)

                                                                   gCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAg                                                                                 >  2:68711/1‑150 (MQ=255)

                                                                      ttGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATCTTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACg                                                                              <  2:1257177/150‑1 (MQ=255)

                                                                        gCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCCTAGCGCTCACCCAGACGCGCCAGAAACTCTGCAGACGGt                                                                            >  1:806326/1‑150 (MQ=255)

                                                                        gCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACggg                                                                            >  2:1413133/1‑149 (MQ=255)

                                                                         ccAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACTGtt                                                                           >  2:494443/1‑150 (MQ=255)

                                                                            gcgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtct                                                                        >  2:551718/1‑150 (MQ=255)

                                                                            gcgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtct                                                                        >  1:1803293/1‑150 (MQ=255)

                                                                            gcgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtcg                                                                        >  1:1487086/1‑149 (MQ=255)

                                                                               cAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg                                                                     <  1:1414919/150‑1 (MQ=255)

                                                                                 aaaTTGCCCCCGATGCAGAGCAATGCTCGAGCCTGCGCCGTAGATATTGCTTGCATGCGGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGAAGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg                                                                   <  2:1218838/150‑1 (MQ=255)

                                                                                 aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCTGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg                                                                   <  1:14189/150‑1 (MQ=255)

                                                                                 aaaTTGCCCCCCATGCACATCAATGCTCGAGCCTGCCGCGTACATATTGCTTGCATGCTGGCAATTG‑CTGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg                                                                   <  1:1338227/150‑1 (MQ=255)

                                                                                       ccccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCg                                                             <  1:26781/150‑1 (MQ=255)

                                                                                       ccccccATGCAGAGCAATCCGCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGGGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCg                                                             <  2:1560803/150‑1 (MQ=255)

                                                                                            cATGCAGATCAATGCACGAGCCTGCCCCGTACATATTGCGTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGGGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGt                                                        <  2:346082/150‑1 (MQ=255)

                                                                                              tGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTc                                                      <  2:1299808/150‑1 (MQ=255)

                                                                                                  gATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTc                                                  >  2:108438/1‑150 (MQ=255)

                                                                                                   aTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCg                                                 <  1:789524/150‑1 (MQ=255)

                                                                                                   aTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGGGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCg                                                 <  2:1865551/150‑1 (MQ=255)

                                                                                                   aTCAATGCGCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCg                                                 <  1:1700517/150‑1 (MQ=255)

                                                                                                       aTGCTCGAGCCTGCCCCGGACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGcccg                                             >  2:558223/1‑149 (MQ=255)

                                                                                                         gCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGCCGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGt                                           >  2:990157/1‑150 (MQ=255)

                                                                                                            cGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACa                                        >  1:1298540/1‑150 (MQ=255)

                                                                                                              aGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACAtt                                      >  1:520525/1‑150 (MQ=255)

                                                                                                               gCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCCCCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTa                                     >  2:1663020/1‑150 (MQ=255)

                                                                                                                  tGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTTTACATTagag                                  <  2:1428582/150‑1 (MQ=255)

                                                                                                                  ggCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTagag                                  <  2:648891/149‑1 (MQ=255)

                                                                                                                    ccccGTACATATTGCTTGCATGCTGTCTATTG‑CTGCATGTCAAGGTGGATGAGTTGGGGTGATGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACTGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgtg                                <  2:50004/150‑1 (MQ=255)

                                                                                                                    cccaGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGGCCAGGGGCCTGAGGTGGGGGGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgtg                                <  2:414976/150‑1 (MQ=255)

                                                                                                                       cGTACATCTTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGGGGGGTGAAGCCATAGCGCGCACCCAGACGAGCCAGAACCTCTGCAGACGGTTGCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTcc                             <  2:1061711/150‑1 (MQ=255)

                                                                                                                       cGTACATATTGCTTGCATGCTGGCATTTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGGCATAGCGCTCACCCAGACGAGCCAGAAAGTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTcc                             <  2:204612/150‑1 (MQ=255)

                                                                                                                           cATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATCAGGTGGGGGGACGCCATAGCGCACACCCAGACGAGCCAGAACCTCTGCAGCCGGTTTCTCGGTGATAACGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGt                         <  2:938735/150‑1 (MQ=255)

                                                                                                                            atatTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTa                        >  1:996192/1‑150 (MQ=255)

                                                                                                                                cctggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGCCGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTgg                    >  1:556197/7‑150 (MQ=255)

                                                                                                                                  tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGc                  >  1:1082570/5‑150 (MQ=255)

                                                                                                                                  tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGc                  >  2:464053/5‑150 (MQ=255)

                                                                                                                                   ggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCa                 >  2:50404/4‑150 (MQ=255)

                                                                                                                                    tgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAg                >  2:1339892/1‑150 (MQ=255)

                                                                                                                                         gCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAcc           >  2:1408205/1‑150 (MQ=255)

                                                                                                                                           tGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGCCGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCg         >  1:1711900/1‑150 (MQ=255)

                                                                                                                                              cAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCAc      <  2:1307747/150‑1 (MQ=255)

                                                                                                                                              cAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCAc      <  1:702801/150‑1 (MQ=255)

                                                                                                                                              cAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGGCGCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCAc      <  2:271319/150‑1 (MQ=255)

                                                                                                                                                aTTG‑CAGGATGTCCAGGTGACGGAGGTGGGGCGAAGCCATAGCGCACACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCTGTGCTACCGACGGTTCGGGCGCCCTGTACATTAGAGGGTCCACGTAGTGGGCAGATACCCGCACCa    <  2:1821239/150‑1 (MQ=255)

                                                                                                                                                 ttG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGCAGCCATAGCGCTCACCCAGACTAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGCCGGGTCGGTAGCCCTGTACATTAGAGTGTCCCCGTAGTGGGCAGATCCCCGCAcccg   >  1:1625277/1‑148 (MQ=255)

AGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAGG  >  NC_000913/1585009‑1585304

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