BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,499,887	G→C	5.0%	pseudogene (648/942 nt)	yjgX ←	pseudogene fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,499,887	0	G	C	5.0%	274.5 / 6.8	100	pseudogene (648/942 nt)	yjgX	pseudogene fragment
Reads supporting (aligned to +/- strand): ref base G (53/40); new base C (2/3); total (56/44)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

CCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAA  >  NC_000913/4499740‑4500034

ccACTGAAAAATACCTGCGAAATGACGGGCTCTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATCAACAGGATTGGCACTGTCGCTGa                                                                                                                                                    <  1:288725/149‑1 (MQ=255)

ccACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGa                                                                                                                                                    <  1:734902/149‑1 (MQ=255)

   cTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCg                                                                                                                                                 >  1:1229596/1‑149 (MQ=255)

      aaaaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACa                                                                                                                                              <  1:843410/149‑1 (MQ=255)

        aaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgt                                                                                                                                            <  1:999903/149‑2 (MQ=255)

         aaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTTGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCAAAATGATGACCAGCATTGGCACTGTCGCTGACCGACAGcc                                                                                                                                           >  2:37882/1‑149 (MQ=255)

         aaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtg                                                                                                                                           <  2:1406960/149‑3 (MQ=255)

          aTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtgc                                                                                                                                          <  2:812199/149‑4 (MQ=255)

                 cGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATACTCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGCTGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCa                                                                                                                                   <  1:259663/149‑1 (MQ=255)

                  gAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGATCAGGATTGGCACTGTCGCTGCCCGACAGCCCAAGGCCAg                                                                                                                                  >  1:216574/1‑149 (MQ=255)

                          gtgCACTGATTTTTCTACAGACGGGCTGATGAGACGTAATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTaa                                                                                                                          <  2:1317500/147‑1 (MQ=255)

                          gggCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGATCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGCAAAGTaa                                                                                                                          >  2:879145/1‑149 (MQ=255)

                                tGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATg                                                                                                                    >  1:733892/1‑149 (MQ=255)

                                tGATTTATCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATg                                                                                                                    <  2:360345/149‑1 (MQ=255)

                                    ttttCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGa                                                                                                                >  2:1453652/1‑149 (MQ=255)

                                    ttttCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGa                                                                                                                >  1:542280/1‑149 (MQ=255)

                                       tCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGaaag                                                                                                             >  1:1237860/1‑148 (MQ=255)

                                       tCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAAt                                                                                                             <  1:1115320/149‑1 (MQ=255)

                                         tACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGa                                                                                                           >  2:1516982/1‑149 (MQ=255)

                                            tGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGt                                                                                                        <  1:82194/149‑1 (MQ=255)

                                              actggCGGATGAGACGTTATGTCAATGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGa                                                                                                      <  2:1056703/146‑1 (MQ=255)

                                                  gCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACTTGAGTAc                                                                                                  <  2:968038/149‑1 (MQ=255)

                                                    tgatgaGACGTTATGTCCCTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCg                                                                                                >  1:819543/1‑149 (MQ=255)

                                                    tgaggaGACGTTATGTCCCTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATCGCGAAATGACGTGAGTACCg                                                                                                >  1:819542/1‑149 (MQ=255)

                                                      atgaGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGc                                                                                              >  1:1344532/1‑149 (MQ=255)

                                                        gagaCGTTTTGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGATAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGa                                                                                            >  2:483484/1‑149 (MQ=255)

                                                        gagaCGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGa                                                                                            >  2:541601/1‑149 (MQ=255)

                                                         tgaCGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGaa                                                                                           <  1:234623/148‑1 (MQ=255)

                                                         acgcGTTTTGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCTTACTGATGAACCGGATTGGCACTGTCGCTGACCGACAGCCGAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGaa                                                                                           <  2:1537460/146‑1 (MQ=255)

                                                            cGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGa                                                                                        <  2:164661/149‑1 (MQ=255)

                                                             gTTATGTCACTGGCGGTAATAATCAGTGGCCCACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTCTCGCTTACCGGCAGCCCTAGGCCAGAAAAGTATGGCAGTGCTAAATTACGTCGGTACCGGCGAAgggt                                                                                       >  2:1237509/1‑147 (MQ=255)

                                                             gTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATTAACAGGATTGGCACTGTCGATGACCGCCAGCCCTAGGACTGAAAAGGAAAGCATGGCGAAATGAGGTGAGTACCGGGGAAGGAt                                                                                       >  2:362255/1‑149 (MQ=255)

                                                               tATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATAc                                                                                     <  2:462245/149‑1 (MQ=255)

                                                                atctCACTAGCGGTAATAATCAGTGGCAAACTGTAGCCTCCCTGCACATGACAATAAAGATGAAAAGGATTGGCACTGTCGCTGACCGACAGCCCAATGCCAGAAAAGAAAAGCATGGCGAAATGACGTGAGTACCGGCGAATGATAcc                                                                                    >  2:801366/4‑149 (MQ=255)

                                                                  gTCACTGGCGGCAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGAGGAACAGGATTGGCACTGGCGCTGACCGACAGCCCAAGGCCAGAAAAGAAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCAt                                                                                  <  2:1275319/149‑1 (MQ=255)

                                                                        ggCGGTAATAATCAGAGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCt                                                                            >  1:1349206/1‑149 (MQ=255)

                                                                         gCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTg                                                                           >  1:552488/1‑149 (MQ=255)

                                                                          cGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTgg                                                                          <  1:541601/149‑1 (MQ=255)

                                                                           ggTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGCCAGCCCAAGGCCAGAAATGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGc                                                                         >  2:1144549/1‑149 (MQ=255)

                                                                                  atCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGCCCGCACAAGGCCAGAAAAGTAATGCATTGCGAAATGACGTGAGTGCCGGCGAAGGATCCCATCAAGATCGCCAAGaca                                                                  >  2:422741/1‑147 (MQ=255)

                                                                                  atCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGaaa                                                                  <  2:1540086/149‑1 (MQ=255)

                                                                                  atCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGaaa                                                                  <  1:816740/149‑1 (MQ=255)

                                                                                  atCAGTGGCAAACTGTAGCCTCCCTGCACATGACCATAATGAGTAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAACCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGaaa                                                                  <  1:1385680/149‑1 (MQ=255)

                                                                                   tCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAg                                                                 <  2:610773/149‑1 (MQ=255)

                                                                                     aGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGtt                                                               <  2:1140008/149‑1 (MQ=255)

                                                                                       tGGCACACTGAACCCTCCCTGCACATGACCATAATGATCCACAGGATTGGCCCAGTCGCGTACCGACAGCCCAAGGCCAGAAAACTAAAGCAAGGCGAAATGAATTGAGTACCGGCGAAGGATACCATCAAGATGGCCAAGAAAGTTAt                                                             <  2:857898/149‑1 (MQ=255)

                                                                                         gCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATcc                                                           >  1:1296570/1‑149 (MQ=255)

                                                                                          cacaCTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAACGCATGGCGAAATGACGGGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCa                                                          <  1:502654/149‑1 (MQ=255)

                                                                                          cacaCTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCa                                                          <  2:188429/149‑1 (MQ=255)

                                                                                               tGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTa                                                     >  2:720229/1‑149 (MQ=255)

                                                                                                gTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTAc                                                    <  1:1251357/149‑1 (MQ=255)

                                                                                                gTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTAc                                                    <  2:1080857/149‑1 (MQ=255)

                                                                                                   gCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGa                                                 >  2:169626/1‑149 (MQ=255)

                                                                                                    ccTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGTGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGAt                                                >  1:405763/1‑149 (MQ=255)

                                                                                                      tCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGc                                              <  1:1242321/149‑1 (MQ=255)

                                                                                                        ccTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTg                                            >  2:1077733/1‑149 (MQ=255)

                                                                                                             tcaTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGGTGGCCAAGAAAGTTATCCAGTTTCCTGATGCTGGCGAg                                       <  1:1204808/148‑1 (MQ=255)

                                                                                                             acaTGACCATAATGATGAACATGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAg                                       <  1:879145/149‑1 (MQ=255)

                                                                                                             acaTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAg                                       <  1:1453652/149‑1 (MQ=255)

                                                                                                               attaCCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGt                                     <  1:1303954/146‑1 (MQ=255)

                                                                                                               aTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGCCGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGt                                     >  1:1412992/1‑149 (MQ=255)

                                                                                                               aTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACAGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGGTATCCATTTTACTGATGCTGGCGAGGt                                     >  2:859417/1‑149 (MQ=255)

                                                                                                                 gACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCAAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATCCCATCAAGCTGGCAAAGAAAGTTATGCAGTTTTCTCATGCTGGAGAGGTaa                                   >  2:1371804/1‑149 (MQ=255)

                                                                                                                 gACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACACCCCAATGCCAGAAAAGTAAAGCATGGCGAAATGAAGTGAGTACCGGCGAAGGATAACATCCAGCTGGCCAAGAAAGTTCTCCCGTTTACGGATGCTGGCGAGGTaa                                   >  1:1023017/1‑149 (MQ=255)

                                                                                                                    cATAATGATGAACAGGATTGGCACTGTCGCTGACCGCCAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAg                                >  2:497304/1‑149 (MQ=255)

                                                                                                                    cATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAg                                >  2:74249/1‑149 (MQ=255)

                                                                                                                    cATAATGATGAACACGATTGGCACTGTCGCTGACCGACAGCGCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAg                                <  1:1063098/149‑1 (MQ=255)

                                                                                                                     aTAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAgg                               <  1:228812/149‑1 (MQ=255)

                                                                                                                      tAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCTAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCAAGTTTACTGATGCTGGCGAGATAACAGGc                              >  1:1505368/1‑149 (MQ=255)

                                                                                                                      tAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGc                              >  2:1541445/1‑149 (MQ=255)

                                                                                                                       aaTGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCa                             >  2:1068644/1‑149 (MQ=255)

                                                                                                                         tgatgaACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAAt                           >  1:685415/1‑149 (MQ=255)

                                                                                                                          gatgaACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACAGGCGAAGGATACCATCAAGATGGCCAAGAAAGTTATCCAGTTTACTGATGATGGCGAGGTAACAGGCTAtt                          >  2:343244/1‑149 (MQ=255)

                                                                                                                           atgaACAGGATTGGCACTGTCGCTGACCGACAGCACAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGAAACCAGCAAGCTGGCCAAGAAAGTTATCAAGTTTACTGATGCTGGCGAGGTAAAAGGCAAttt                         >  1:1127660/1‑149 (MQ=255)

                                                                                                                            tgaACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAAtttt                        >  1:748059/1‑149 (MQ=255)

                                                                                                                              aaCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTCACAGGCATATTTTc                      >  2:200501/1‑149 (MQ=255)

                                                                                                                              aaCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTc                      >  1:277866/1‑149 (MQ=255)

                                                                                                                              aaCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTc                      >  2:516600/1‑149 (MQ=255)

                                                                                                                               accggtatGGCACTTTGCCTGTCCGACTTCCCAAGGCCATAAAAGTAAAGCCTTGGGAAATGCCGTGTGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCg                     <  2:266653/142‑1 (MQ=255)

                                                                                                                               aCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCg                     <  2:364006/149‑1 (MQ=255)

                                                                                                                                cAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACAGAAGCTGGCGAGGTCACAGGCAATTTTTCgg                    >  2:1530019/1‑149 (MQ=255)

                                                                                                                                cAGCATTGGCACTGTCGCTGACCGACAGCCACAGGCCAGAAAAGTAAAGCATGGCGAATTGCCGTGAGTACCGGCGAAGGAGACCATGACGATAGCCAAGATCGTTATCCAGTTTACTGATCCTGGCGAGGTAACAGGCCATTTTTctg                    >  2:846792/1‑147 (MQ=255)

                                                                                                                                 aGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCggg                   <  1:74249/149‑1 (MQ=255)

                                                                                                                                 aGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCggg                   <  2:1296570/149‑1 (MQ=255)

                                                                                                                                   gATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGAGGCTGGCGAGGTAACCGGCAATTTTTCGGGGa                 >  1:1158870/1‑149 (MQ=255)

                                                                                                                                     ttGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAACGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACAATCAAGCTGGACAAGAAAGTTATACAGTTTACTGATGCTGTCGAGGTAACAGGCAATTTTTCGTGGATa               >  2:1547739/1‑149 (MQ=255)

                                                                                                                                     gtGGCACAGTCTCTGACCGACAGCCCAAGGCCAGAAAACTAAAGCATCGCGAAATGACGTGAGTACCGGCGAAGGAAACCATCAAGCTGGCCAAGAAAGTTATCCAGTTAACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATa               <  1:37882/148‑1 (MQ=255)

                                                                                                                                     cttgggcTGTCGGTCAGCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATa               >  1:1269519/7‑149 (MQ=255)

                                                                                                                                      ttgggcTGTCGGTCAGCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTTCCGGCGAAGGTTACCAGCAAGGAGGCCAACAAAGTTATACAGTGGACTGATGCTGGCGAGGTAACTGGCATTTTTTTGGGGATAc              >  2:1540725/6‑149 (MQ=255)

                                                                                                                                      tGGCACTGTCGCGGCCCGACAGCCCAAGGCCAGAAAAGTAACGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATAc              <  2:1237860/149‑1 (MQ=255)

                                                                                                                                       ggCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACt             >  1:643519/1‑149 (MQ=255)

                                                                                                                                          aatGTCGCTGACCGACAGCCGAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCTAAGGATATCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGATGTAAGAGGCAATTTATCGGGGATACTGCt          <  1:801366/147‑1 (MQ=255)

                                                                                                                                          aCTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGCTACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGCTGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCt          >  1:70322/1‑149 (MQ=255)

                                                                                                                                          aCTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGcg          >  1:100456/1‑148 (MQ=255)

                                                                                                                                            tGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTcc        >  2:1253318/1‑149 (MQ=255)

                                                                                                                                                gCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGt    >  1:69871/1‑149 (MQ=255)

                                                                                                                                                 cTGCCCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCAACAAGCTGGCCAACAAATTTATCCAGTATACAGAAGCTGGCGAGCTAACAGTCAATTTATCGGGCATACTGCTCCAggaa   >  1:627680/1‑147 (MQ=255)

                                                                                                                                                 cTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTa   >  1:597695/1‑149 (MQ=255)

                                                                                                                                                  tGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAATTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAGCAGGCAATTTTTCGGGGATACTGCTCCAGGGaa  >  2:1087009/1‑149 (MQ=255)

                                                                                                                                                  tGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTaa  <  1:1516982/149‑1 (MQ=255)

CCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAA  >  NC_000913/4499740‑4500034

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