BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,499,885	C→G	5.1%	pseudogene (650/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,885	0	C	G	5.1%	257.8 / 6.9	98	pseudogene (650/942 nt)	yjgX	pseudogene fragment
Reads supporting (aligned to +/- strand): ref base C (50/42); new base G (2/3); total (53/45)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

Only 100 of 102 total aligned reads displayed.

CATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTA  >  NC_000913/4499737‑4500033

cATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGc                                                                                                                                                      <  2:506645/149‑1 (MQ=255)

aaTCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGGGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGc                                                                                                                                                      <  2:498158/148‑1 (MQ=255)

 atacACTTAAAAAAACTTGCGAAATGCCGGGCTCTGATTTTTATAATTACGGGCTGAAGCGACGTTATTTCACTGTCGGTAATAATCAGTGGCACACTGTAGCCTCCCTTCACATGACCATAATGATGAACAGGATTGGCACTGTCGCt                                                                                                                                                     <  2:1029851/146‑1 (MQ=255)

 aTCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCt                                                                                                                                                     <  1:820118/149‑1 (MQ=255)

   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)

            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)

                                                       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                     >  2:516600/1‑149 (MQ=255)

                                                                                                                                 aaCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTc                     >  1:277866/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)

CATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTA  >  NC_000913/4499737‑4500033

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