BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,477,885	C→A	5.1%	I13I (ATC→ATA)	dsdX →	D‑serine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,477,885	0	C	A	5.1%	301.0 / 5.1	107	I13I (ATC→ATA)	dsdX	D‑serine transporter
Reads supporting (aligned to +/- strand): ref base C (53/48); new base A (4/2); total (57/50)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 108 total aligned reads displayed.

ATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAAC  >  NC_000913/2477736‑2478031

aTCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATc                                                                                                                                                    <  2:695897/150‑1 (MQ=255)

 tCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCa                                                                                                                                                   <  1:1391833/150‑1 (MQ=255)

       cATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCg                                                                                                                                             <  1:1700123/150‑1 (MQ=255)

          tGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGATGCTGATAAGCAGCgtgc                                                                                                                                          <  1:1670583/150‑2 (MQ=255)

           gATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCACATCT‑GGGTTGTGAGCACGCTGGTTATCAGCAACGTGtt                                                                                                                                         >  2:313062/1‑150 (MQ=255)

           gATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGtt                                                                                                                                         >  1:143386/1‑150 (MQ=255)

           gATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGtt                                                                                                                                         <  2:729166/150‑1 (MQ=255)

           gATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGtt                                                                                                                                         <  1:531146/150‑1 (MQ=255)

              gCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAAt                                                                                                                                      <  1:412406/150‑1 (MQ=255)

               cAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCTTTGTTTTTTTGTCCAATATCATCAGGTTAATAACAGGGGAAGGTGAGATATGACATCTCAAATCT‑GCGTTGTGAGCACGCTGCTTATCAGCATCCTGATAAtt                                                                                                                                     >  2:838583/1‑150 (MQ=255)

                aaTTTACCTCATGCGAAAGGCAAATTTTATCGTATGGCAGCCTCCCTTGTTTTTTTGTCCAATATCATCAGGGTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTGGTGAGCACGCTGCTTATCAGCATCGTGTTAATTg                                                                                                                                    <  2:1596245/150‑1 (MQ=255)

                      ccTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGa                                                                                                                              >  2:1111079/1‑150 (MQ=255)

                        tCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGAcc                                                                                                                            >  1:1628099/1‑150 (MQ=255)

                            gtgAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTAAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCg                                                                                                                        <  2:1010486/150‑1 (MQ=255)

                             tgATAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGt                                                                                                                       <  1:1443313/150‑1 (MQ=255)

                              gAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTg                                                                                                                      <  1:950387/150‑1 (MQ=255)

                              gAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTg                                                                                                                      >  2:528262/1‑150 (MQ=255)

                               aaaGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGa                                                                                                                     <  2:1592244/150‑1 (MQ=255)

                                 aGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAg                                                                                                                   >  1:497726/1‑150 (MQ=255)

                                  ggCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGt                                                                                                                  >  1:690951/1‑150 (MQ=255)

                                   gCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGtt                                                                                                                 >  1:216036/1‑150 (MQ=255)

                                   gCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCAAAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGAACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGtt                                                                                                                 >  2:6272/1‑150 (MQ=255)

                                         tttATCGTTAGTCAGCCTGCGTTGTGTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAAtt                                                                                                           <  2:1216826/150‑1 (MQ=255)

                                           tATCGTTTGTCAGCATGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTcc                                                                                                         <  2:1394512/150‑1 (MQ=255)

                                             tCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCAc                                                                                                       <  1:1338792/150‑1 (MQ=255)

                                             tCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCAc                                                                                                       <  2:143386/150‑1 (MQ=255)

                                                  tGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTTATCACAGGGGAAGGTGAGATATGGACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTTTCAGCATCGTGTTAATTGTTCTGACCATCGTGAAGTTCTAATTCCACCCGtt                                                                                                  >  1:860014/1‑150 (MQ=255)

                                                    tCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATCTGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACTATCGTGAAGTTCAAATTCCACCCGTTTc                                                                                                <  1:771769/150‑1 (MQ=255)

                                                       gCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGATGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTgg                                                                                             >  1:921517/1‑150 (MQ=255)

                                                         cTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGgcg                                                                                           <  1:998461/150‑1 (MQ=255)

                                                           gCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTTATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCt                                                                                         >  1:1433944/1‑150 (MQ=255)

                                                            cgatgttTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGAACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTg                                                                                        >  2:1150512/4‑150 (MQ=255)

                                                              ttgttTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACAATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGtt                                                                                      >  2:883928/1‑150 (MQ=255)

                                                                gttTTTTTGTCGAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGc                                                                                    <  2:1488201/150‑1 (MQ=255)

                                                                    ttttGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGc                                                                                <  1:920566/150‑1 (MQ=255)

                                                                     tttGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGcc                                                                               <  2:1601213/150‑1 (MQ=255)

                                                                         tCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGct                                                                           <  2:235291/150‑1 (MQ=255)

                                                                         tCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGct                                                                           <  2:439451/150‑1 (MQ=255)

                                                                          ccAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGctt                                                                          >  2:415895/1‑150 (MQ=255)

                                                                          ccAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGctt                                                                          >  1:787259/1‑150 (MQ=255)

                                                                          ccAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGctt                                                                          >  2:632692/1‑150 (MQ=255)

                                                                           cAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGcttc                                                                         <  2:161900/150‑1 (MQ=255)

                                                                            aaTATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGcttct                                                                        >  2:311147/1‑150 (MQ=255)

                                                                               atcatcAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCg                                                                     <  2:1357773/150‑1 (MQ=255)

                                                                               atcatcAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCg                                                                     <  2:473012/150‑1 (MQ=255)

                                                                               atcatcAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCg                                                                     <  2:954364/150‑1 (MQ=255)

                                                                               atcatcAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCg                                                                     <  1:1541926/150‑1 (MQ=255)

                                                                                 tatcaGGTTTATCACTGGGGAAGGTGAGATTTGCACTCTCAAATCT‑GGGTTGTGAGCACTCTGCTTATCAGCATCGTGTTAATTGTACTGAACATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTg                                                                   <  2:156194/149‑1 (MQ=255)

                                                                                 catcaGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAAACT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTTGCGCTGTTGCTGGCAAGATTCTTCGTg                                                                   >  2:613280/1‑150 (MQ=255)

                                                                                  atcaGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTgg                                                                  >  2:11300/1‑150 (MQ=255)

                                                                                    caGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCCTCGTGGGa                                                                >  2:743436/1‑150 (MQ=255)

                                                                                     aGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCGGGCCAGCTTCTTCGTGGGaa                                                               <  1:311524/150‑1 (MQ=255)

                                                                                         tAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgat                                                           >  1:442614/1‑150 (MQ=255)

                                                                                          aaTCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGCACgatg                                                          >  2:855624/1‑150 (MQ=255)

                                                                                            tCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgatgat                                                        <  1:604866/150‑1 (MQ=255)

                                                                                                aGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGc                                                    >  2:1028165/1‑150 (MQ=255)

                                                                                                 ggggAAGGTGAGATATGCACTCTCAAATCT‑GGGGTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCa                                                   <  2:1531815/150‑1 (MQ=255)

                                                                                                     aaGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTTCTGACCATCGTGAAGTTCTAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATggg                                               <  1:1150512/150‑1 (MQ=255)

                                                                                                      aGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATgggg                                              <  1:1245236/150‑1 (MQ=255)

                                                                                                      aGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATgggg                                              <  2:527912/150‑1 (MQ=255)

                                                                                                        gTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGcc                                            >  1:552432/1‑150 (MQ=255)

                                                                                                         tGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCa                                           >  1:853731/1‑150 (MQ=255)

                                                                                                             atatGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTgg                                       >  2:1704703/1‑150 (MQ=255)

                                                                                                             atatGCACTCTCAAATCT‑GGGTGGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTgg                                       <  2:1323768/150‑1 (MQ=255)

                                                                                                               atGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGat                                     >  2:785193/1‑150 (MQ=255)

                                                                                                               atGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGat                                     >  2:1527719/1‑150 (MQ=255)

                                                                                                                tGCACTCTCAAATCTGGGGTTGTGAGCCCGCTGCTTATCAGCATCGTGTTAATTGTACTGACCAGCGTGAAGTGCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGat                                     <  2:1176610/150‑1 (MQ=255)

                                                                                                                 gCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGATGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGatat                                   >  1:250983/1‑150 (MQ=255)

                                                                                                                  cACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATg                                  >  1:1180164/1‑150 (MQ=255)

                                                                                                                   aCTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATgg                                 <  1:1068572/150‑1 (MQ=255)

                                                                                                                    ctctCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGt                                >  2:1054565/1‑150 (MQ=255)

                                                                                                                    ctctCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGATTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGt                                >  1:692766/1‑150 (MQ=255)

                                                                                                                     tctcAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTACTTCGTGGCAACGATGATGGGCATGGGGCCACTGGATATGGTa                               >  2:804948/1‑150 (MQ=255)

                                                                                                                          aactcggggTTGTGAGCCCGCAGCTTATCAGCATCGTGTTAATTCTACTGACCCTCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCGGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATg                           <  1:854573/144‑1 (MQ=255)

                                                                                                                             cT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTAt                       >  1:1672513/1‑150 (MQ=255)

                                                                                                                             cT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGATGCTGGCCAGATTCTTAGTGGGAACGATGATTGGCATGGGGACACTGGATATGGTAAATGCTAt                       >  2:1045719/1‑150 (MQ=255)

                                                                                                                                gggTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTg                     >  2:1198558/1‑150 (MQ=255)

                                                                                                                                 ggTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGa                    >  1:1292593/1‑150 (MQ=255)

                                                                                                                                  gTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGaa                   >  1:1265591/1‑150 (MQ=255)

                                                                                                                                  gTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCAAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGaa                   >  1:1070766/1‑150 (MQ=255)

                                                                                                                                    tgtgAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAg                 <  1:528262/150‑1 (MQ=255)

                                                                                                                                    aattaacACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAg                 >  2:408803/7‑150 (MQ=255)

                                                                                                                                     gtgAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGt                <  1:1340905/150‑1 (MQ=255)

                                                                                                                                      ttaacACGATGCTGATAAGCAGCGTGGTAATTGTACTGACCATCGTGAAGTTCAAATTCCACACGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGAGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg               >  1:899378/5‑150 (MQ=255)

                                                                                                                                      tgAGCACGCTGCTTATCAGCATCGTGTTAATTGTCCTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg               <  1:360313/150‑1 (MQ=255)

                                                                                                                                      tgAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg               >  1:1462973/1‑150 (MQ=255)

                                                                                                                                        cgCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACGATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGTTAAATGCTATTGAAAGTTGa             <  1:6272/149‑1 (MQ=255)

                                                                                                                                          cACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATAAGGTAAATGCTATTGAAAGTGGAAt           >  1:932018/1‑150 (MQ=255)

                                                                                                                                          cACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg               <  1:1379290/146‑1 (MQ=255)

                                                                                                                                          cACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg               >  2:1379290/1‑146 (MQ=255)

                                                                                                                                          cACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAAt           >  1:320982/1‑150 (MQ=255)

                                                                                                                                            cGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTg         >  1:164607/1‑150 (MQ=255)

                                                                                                                                             gctgctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTgg        <  1:1011993/150‑1 (MQ=255)

                                                                                                                                              ctgctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc       <  2:818741/150‑1 (MQ=255)

                                                                                                                                              ctgctTATCAGCATCGTGTTAATGGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTCCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc       <  2:818738/150‑1 (MQ=255)

                                                                                                                                                gctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCgg     <  2:1243658/150‑1 (MQ=255)

                                                                                                                                                gctTATCAGCATCGGGGTCATTGTACTGACCATCGTGAAGTTCACATGCCACCCGTTTCTGGCGCTGGTGCTGGCCAGCTTCTTCGTGGGAACGATGAGGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCgg     <  2:830339/150‑1 (MQ=255)

                                                                                                                                                   tATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAAc  >  1:1187844/1‑150 (MQ=255)

                                                                                                                                                   tATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAAc  >  1:1530242/1‑150 (MQ=255)

                                                                                                                                                   tATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCAGGCTCTTTTGCTGGCAAGCTTCTTCGTGGTAACGATGATGTGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAAc  >  1:1070805/1‑150 (MQ=255)

ATCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCT‑GGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAAC  >  NC_000913/2477736‑2478031

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