BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,219	T→A	5.5%	I385F (ATT→TTT) ‡	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,219	0	T	A	5.5%	279.8 / 2.6	93	I385F (ATT→TTT) ‡	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (66/21); new base A (4/2); total (70/23)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGC  >  NC_000913/3728080‑3728368

aGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCAGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa                                                                                                                                             <  1:865690/150‑1 (MQ=255)

aGCCTTGCGGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGAAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCAAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa                                                                                                                                             <  1:149397/150‑1 (MQ=255)

 gCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCTCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGGATTAAACCGCaaagtg                                                                                                                                            <  2:1387614/150‑6 (MQ=255)

 gCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaagtg                                                                                                                                            <  2:1647525/150‑6 (MQ=255)

  ccTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTTATATCCTCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATg                                                                                                                                           <  1:1350887/150‑1 (MQ=255)

                ttGCTGTCCCCACAT‑CCCCCCCCTTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                                                                                                                             <  2:148998/149‑1 (MQ=255)

                gTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAGGc                                                                                                                             <  1:1608271/150‑1 (MQ=255)

                      gCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                                                                                                                       <  1:693522/150‑1 (MQ=255)

                      gCCCCACAA‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                                                                                                                       <  1:1065349/150‑1 (MQ=255)

                       ccccACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                                                                                                                       <  2:793499/150‑1 (MQ=255)

                            caT‑CCCACCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCTGGTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg                                                                                                                 <  1:761463/150‑1 (MQ=255)

                                 cccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCCTCTGCCAGCGCTTa                                                                                                             >  1:645076/1‑150 (MQ=255)

                                     cccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCACCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAccc                                                                                                         <  1:119558/150‑1 (MQ=255)

                                      ccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa                                                                                                        <  1:50718/150‑1 (MQ=255)

                                           cGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGTGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCt                                                                                                   >  2:741038/1‑150 (MQ=255)

                                           cGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATACCCCCCGcaa                                                                                                   >  1:695858/1‑148 (MQ=255)

                                           cGGTAATCGAGCTGCTGACCGCTGATATACGCCAGCATCTGACGCCAGTACTCCCTACGCGCCCCGCCCCCAATCAACGTAACCCTTTGCGGTTTATACCCGCAGGCATGCCCGACATCCATGCCATCTGCAAGAGCATTTCCAACGcgg                                                                                                   >  2:1188349/1‑148 (MQ=255)

                                             gTAATCGAGCTGCTGCCCGCTGATATCCGCCAGCATATGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATACCCCCCGCCTTc                                                                                                 >  2:1272347/1‑150 (MQ=255)

                                              tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCCCGCCTTcc                                                                                                >  1:86509/1‑150 (MQ=255)

                                              tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTcc                                                                                                >  2:675312/1‑150 (MQ=255)

                                              tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTcc                                                                                                <  1:89673/150‑1 (MQ=255)

                                                   gAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCTCATAACCCCCGCCTTCCAGCAc                                                                                           >  1:453823/1‑150 (MQ=255)

                                                    aGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGATTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCAcc                                                                                          <  1:214988/150‑2 (MQ=255)

                                                     gctgctGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACCTCCATGCCATCTGCCATCGCATAACCCACGCCTTCCAGCACTg                                                                                         <  1:364497/150‑1 (MQ=255)

                                                      ctgctgACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAGCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGc                                                                                        >  1:1178248/1‑150 (MQ=255)

                                                      ctgctgACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGCCATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGc                                                                                        >  2:1663544/1‑150 (MQ=255)

                                                         ctgACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCCTGCCATCTGCCGGCGCCTAACCCCCGCCTTCCCGCCCTGcgcg                                                                                     >  2:1231873/1‑150 (MQ=255)

                                                           gACCGCTGATATCCGCCAGCATGTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCAGTTTAATACCGCAGGCATGCACGACTTCCATGCCAGCTGCCAGCGCATAACCCACGCCTTCCTGCACTGCTcgcg                                                                                   <  2:1724694/150‑1 (MQ=255)

                                                           gACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTcgcg                                                                                   >  2:1039725/1‑150 (MQ=255)

                                                           gACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTcgcg                                                                                   >  2:1326853/1‑150 (MQ=255)

                                                            aCCGCTGATATCCGCCAGCATCTGACGCCTGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCgcgc                                                                                  >  2:1236239/1‑150 (MQ=255)

                                                                cTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATAATTCCGCAGGCATGCACGACCTCCATGCCATCTGCCAGCGCATAACCCACTCCTGCCAGCACTGCTCGCTcccgg                                                                              >  1:263444/1‑147 (MQ=255)

                                                                  gATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTc                                                                            <  2:646350/150‑1 (MQ=255)

                                                                  gATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTc                                                                            <  1:1577615/150‑1 (MQ=255)

                                                                    tatCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAt                                                                          >  2:1243456/1‑150 (MQ=255)

                                                                     atCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCttt                                                                         >  2:310220/1‑150 (MQ=255)

                                                                     atCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAtt                                                                         >  2:471198/1‑150 (MQ=255)

                                                                      tCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCCTCTGCCAGCGCATACCCAACGCCTTCCAGCACTGCTCGCGCCAGTTCATTg                                                                        >  1:1213094/1‑150 (MQ=255)

                                                                      tCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCCCGCCTTCCAGCACTGCTCGCGCCAGTTCATTg                                                                        >  2:66514/1‑150 (MQ=255)

                                                                       ccGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAACAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGACAGCGCATAACCCACGCCATCCAGCACTGCTCGCGCCAGTTCATTgg                                                                       >  1:614940/1‑150 (MQ=255)

                                                                        cGCCAGCATCTGACGCCAGTACTCCCTACGCGCCCCGCCCCCAATAAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCCTGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTggg                                                                      >  2:811062/1‑150 (MQ=255)

                                                                        cGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACCTCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTggg                                                                      >  1:1817420/1‑150 (MQ=255)

                                                                         gCCTGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTgggg                                                                     <  1:1683548/150‑1 (MQ=255)

                                                                            aGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCa                                                                  >  2:1687713/1‑150 (MQ=255)

                                                                                 cTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtg                                                             >  1:787236/1‑150 (MQ=255)

                                                                                 cTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGGtg                                                             >  2:1439133/1‑150 (MQ=255)

                                                                                  tGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCCCGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCAGGTtga                                                            >  2:690088/1‑150 (MQ=255)

                                                                                     cGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATTCCGCAGGCATGCACGACATCCCTGCCATCTGCCAGCGCATAACCAACGCCTTCCAGCACTGCTCGCGCCAGATCATTGGGGCCATGTtgatga                                                         >  1:609788/1‑150 (MQ=255)

                                                                                     cGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtgatga                                                         >  2:755525/1‑150 (MQ=255)

                                                                                     cGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtgatga                                                         >  2:954740/1‑150 (MQ=255)

                                                                                     cGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTtgatga                                                         >  1:947514/1‑150 (MQ=255)

                                                                                      gCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAg                                                        >  1:741379/1‑150 (MQ=255)

                                                                                       ccAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACCCTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACACGCGCCGTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGt                                                       >  1:912905/1‑150 (MQ=255)

                                                                                       ccAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGt                                                       >  1:1525779/1‑150 (MQ=255)

                                                                                             cTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGCCATCCATGCCATCTGCCAGCGCATAACCCCCGCCTTCCCGCCCTGCTCGCGCCAGTTCATTGGGGCCATGGGGATGAGTCAAAcc                                                 >  1:749067/1‑150 (MQ=255)

                                                                                              tCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACAGCCATGCCATCTGCCAGCGCATAACCCACGCCGTCCAGCACTGCTCGCGCGAGTTCATTGGGGCCAGGTTGATGGGTCAAAccc                                                >  2:4717/1‑149 (MQ=255)

                                                                                               cACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATCCCGCAGGCATGGCCGGCGGCCCTGCCATCTGCCGGCGCCTACCCCCCGCCTGCCAGCACTGCTCGCGCCAGTTCCGTCGGGCCAGGGTGCTGAGTCAGCCCaa                                               >  1:1111790/1‑150 (MQ=255)

                                                                                               cACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCCCGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaa                                               >  2:1343523/1‑150 (MQ=255)

                                                                                                aCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaaa                                              <  2:1734668/150‑1 (MQ=255)

                                                                                                aCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaaa                                              >  1:1734668/1‑150 (MQ=255)

                                                                                                 cTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAg                                             >  2:1805838/1‑150 (MQ=255)

                                                                                                    cgcgCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCAGGTTGATGAGTCAAACCAAAGaaa                                          >  1:1557476/1‑150 (MQ=255)

                                                                                                    cgcgCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGGTCATTGGGGCCATTTTGAGGCGTCAAACcaaacaaa                                          >  2:325152/1‑150 (MQ=255)

                                                                                                        ccccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAAGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAGTGGGGCCATGTTGCTGAGTCAAACCAATGAAAAccc                                      >  1:308428/1‑150 (MQ=255)

                                                                                                         cccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAAGCCATCAGCCAGCGCATAACCCACGCCTTCACGCACAGCTCGCGCCAGTTCATAGGGGCCAAGTTGATGAGTCAAAGAAAATAAAGcccc                                     >  1:872875/1‑150 (MQ=255)

                                                                                                         cccgccccCAATCAACGTAACACTTTGCGGATTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCGAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAcccc                                     >  2:1069435/1‑150 (MQ=255)

                                                                                                          ccgccccCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGCCATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAAccccc                                    >  2:234508/1‑150 (MQ=255)

                                                                                                             cccccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATACCCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTc                                 >  1:558519/1‑150 (MQ=255)

                                                                                                             cccccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCAT‑CCCCCCGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCg                                >  2:1564215/1‑150 (MQ=255)

                                                                                                              ccccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCg                                >  1:1776076/1‑150 (MQ=255)

                                                                                                              ccccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCg                                >  2:593204/1‑150 (MQ=255)

                                                                                                                  aaTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCCCGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTg                            >  1:1614756/1‑150 (MQ=255)

                                                                                                                    tCAACGTCACACTTTGCGGTTTAATACCGCAGGCATGCACGACAGCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTggg                          <  2:1664655/150‑1 (MQ=255)

                                                                                                                      aaCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGa                        >  1:1137272/1‑150 (MQ=255)

                                                                                                                        cGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAACACCCCCTTCGCCTGGGGatt                      >  1:1807060/1‑150 (MQ=255)

                                                                                                                        cGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCAGCCAGCGCATAACCCACGCCTTCCAGCACTGCTAGCGCCAGTACATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTACGCCTGGGGatt                      >  1:652610/1‑150 (MQ=255)

                                                                                                                          tAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGattat                    >  2:1614295/1‑150 (MQ=255)

                                                                                                                          tAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCAAGTTGAAGAGTCACACCAAAGAAAACCCCCTTAGCCTGGGGattat                    >  2:1582809/1‑150 (MQ=255)

                                                                                                                            catgcctGCGGTATTAAACCGCAGGCATGCACGACATCCATGCCATCTCCCAGAGCATAACGCACGCCTTCCAGCACTGCGAGCGCCAGTTCAGTGGGGCCATGTTGATGAGTCAACCCAAAGAAACCCGCCTTCGCCGGGGGATGATtg                  >  1:280049/7‑150 (MQ=255)

                                                                                                                            acacTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtg                  >  1:638586/1‑150 (MQ=255)

                                                                                                                             cacTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgt                 >  2:175893/1‑150 (MQ=255)

                                                                                                                               gcctGCGGTATTAAACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATtgtgt               >  2:818217/4‑150 (MQ=255)

                                                                                                                               gcctGCGGTATTAAACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCGGGGGATTATtgtgt               >  2:825220/4‑150 (MQ=255)

                                                                                                                                 ttGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTgg             <  1:1039725/150‑1 (MQ=255)

                                                                                                                                   gCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCg           >  2:281984/1‑150 (MQ=255)

                                                                                                                                    cGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGt          >  1:1365252/1‑150 (MQ=255)

                                                                                                                                    cGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGt          >  2:1373904/1‑150 (MQ=255)

                                                                                                                                      gTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGTGATTATTGTGTGGCGTAc        >  2:251229/1‑150 (MQ=255)

                                                                                                                                        ttAATACCGCAGGCATGCACGACCTCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCCCTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCGTCGCCTGGGGAGTATTGTGTGGCGTACGc      >  2:1225733/1‑150 (MQ=255)

                                                                                                                                        ttAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCTGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAACCCCCCGTCGCCGGGGGATTATTGTGTGGCGTACGc      >  2:1219656/1‑150 (MQ=255)

                                                                                                                                        ttAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAACCCCCCTTCGCCTGGGGATTTTTGTGTGGCGTACGc      >  1:1146372/1‑150 (MQ=255)

                                                                                                                                          aaTACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCCAACCAAAGAAAACCCCCTTCGCCAGGGGATTATTGGGTGGCGTACGCTc    >  1:1136854/1‑150 (MQ=255)

                                                                                                                                          aaTACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTc    >  2:1400458/1‑150 (MQ=255)

                                                                                                                                            tACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCCGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGCTTATTGTGTGGCGTACGCTCGc  >  1:980233/1‑150 (MQ=255)

AGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGC  >  NC_000913/3728080‑3728368

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