BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,217	A→T	5.4%	I385I (ATT→ATA) ‡	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,217	0	A	T	5.4%	278.9 / 2.9	91	I385I (ATT→ATA) ‡	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base A (65/21); new base T (3/2); total (68/23)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTC  >  NC_000913/3728069‑3728366

cgATCTGCGCCAGCCTTGCTGCGCCCAGTGGAGGCCCCACAT‑ACCCCCCCGTACGGTAATCTAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                                                                                                                                                      <  1:1551407/150‑1 (MQ=255)

cgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                                                                                                                                                      <  1:726501/150‑1 (MQ=255)

           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                                                                                              <  1:89673/150‑1 (MQ=255)

                                                         tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTcc                                                                                              >  2:675312/1‑150 (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:1326853/1‑150 (MQ=255)

                                                                      gACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTcgcg                                                                                 >  2:1039725/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                                            >  1:1734668/1‑150 (MQ=255)

                                                                                                           aCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCaaa                                            <  2:1734668/150‑1 (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                              >  2:593204/1‑150 (MQ=255)

                                                                                                                         ccccAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCg                              >  1:1776076/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)

CGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTC  >  NC_000913/3728069‑3728366

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