BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,499,885	C→G	5.0%	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.0%	385.3 / 9.2	120	pseudogene (650/942 nt)	yjgX	pseudogene fragment
Reads supporting (aligned to +/- strand): ref base C (54/58); new base G (2/4); total (58/62)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

Only 100 of 120 total aligned reads displayed.

CATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGG  >  NC_000913/4499737‑4500031

cATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCt                                                                                                                                                   >  2:1692501/1‑150 (MQ=255)

    cACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGAcc                                                                                                                                               <  1:1720059/150‑1 (MQ=255)

    cACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGAcc                                                                                                                                               <  2:1236765/150‑1 (MQ=255)

    cACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGAcc                                                                                                                                               >  2:439009/1‑150 (MQ=255)

     aCTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCg                                                                                                                                              >  1:540167/1‑150 (MQ=255)

     aCTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTAGCTGACCg                                                                                                                                              >  1:108617/1‑150 (MQ=255)

      cTGAAAAATACCAGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGCCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGa                                                                                                                                             >  1:747180/1‑150 (MQ=255)

        gAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACa                                                                                                                                           >  1:1575970/1‑150 (MQ=255)

        gAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACa                                                                                                                                           >  2:1596/1‑150 (MQ=255)

            aaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtgc                                                                                                                                       <  2:1498201/150‑4 (MQ=255)

             aTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtgcc                                                                                                                                      <  2:1576202/150‑5 (MQ=255)

              tACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtgcca                                                                                                                                     <  1:1322508/150‑6 (MQ=255)

              tACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtgcca                                                                                                                                     <  1:898602/150‑6 (MQ=255)

              tACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCaa                                                                                                                                     >  1:1008575/1‑150 (MQ=255)

                ccTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAgg                                                                                                                                   >  1:1370269/1‑150 (MQ=255)

                ccTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAgg                                                                                                                                   <  2:1262337/150‑1 (MQ=255)

                ccTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAgg                                                                                                                                   <  2:511266/150‑1 (MQ=255)

                  tGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGcc                                                                                                                                 <  1:1702404/150‑1 (MQ=255)

                   gCGAAATGACGGGCACTGATTTTTCTACTGACTGGCTGATGAGACGTTATGTCACTGGCGGTAATAATAAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCa                                                                                                                                <  2:1298350/150‑1 (MQ=255)

                   gCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCa                                                                                                                                >  1:1016206/1‑150 (MQ=255)

                    cGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAg                                                                                                                               >  2:792349/1‑150 (MQ=255)

                    cGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCAAAATGATGAACAGGATTGGCACTGTCGATGACCGACAGCCCAAGGCCAg                                                                                                                               >  2:189856/1‑150 (MQ=255)

                         tGACGGGCAGTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAg                                                                                                                          <  2:656207/150‑1 (MQ=255)

                         tGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAg                                                                                                                          >  2:1688666/1‑150 (MQ=255)

                         tGAAGGGCACTGATTATTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACACGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAg                                                                                                                          <  1:23624/150‑1 (MQ=255)

                           actggCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTa                                                                                                                        <  1:411935/147‑1 (MQ=255)

                           aCGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTa                                                                                                                        >  1:1261647/1‑150 (MQ=255)

                            cGGGCACTGATTATTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTaa                                                                                                                       <  2:318169/150‑1 (MQ=255)

                                   tGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATgg                                                                                                                <  1:1278954/150‑1 (MQ=255)

                                   tGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATgg                                                                                                                <  1:875881/150‑1 (MQ=255)

                                    gATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAACGCATGGc                                                                                                               >  1:137480/1‑150 (MQ=255)

                                    gATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGc                                                                                                               >  2:733431/1‑150 (MQ=255)

                                    gATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGc                                                                                                               >  1:1569341/1‑150 (MQ=255)

                                    gATTTTTCTACTGACGGGATGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGAACGACAGCCCAATGCAAGAAAAGTAAAGCATGGc                                                                                                               >  1:244435/1‑150 (MQ=255)

                                    gATTTTTCTACTGACCGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGc                                                                                                               <  2:453865/150‑1 (MQ=255)

                                       ttttCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGaa                                                                                                            >  1:946227/1‑150 (MQ=255)

                                          tCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATg                                                                                                         >  2:1298805/1‑150 (MQ=255)

                                           cTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGa                                                                                                        <  1:1512348/150‑1 (MQ=255)

                                               tGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTg                                                                                                    <  1:553968/150‑1 (MQ=255)

                                                 aCGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTTGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACACGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAg                                                                                                  <  1:1543641/150‑1 (MQ=255)

                                                 aCGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAg                                                                                                  <  1:1606704/150‑1 (MQ=255)

                                                   gggCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTa                                                                                                >  2:676559/1‑150 (MQ=255)

                                                   gggCTGATGAGACGTTATCTCACTGGCGGTAATAATCAGTGGCACACTGTAGACTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTa                                                                                                <  1:75130/150‑1 (MQ=255)

                                                      cTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCg                                                                                             >  2:392564/1‑150 (MQ=255)

                                                       tgatgaGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCgg                                                                                            <  2:1168581/150‑1 (MQ=255)

                                                        gatgaGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGc                                                                                           >  2:1648750/1‑150 (MQ=255)

                                                        gatgaGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGc                                                                                           >  1:1014286/1‑150 (MQ=255)

                                                            agaCGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAg                                                                                       <  2:1016206/150‑1 (MQ=255)

                                                               cGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGAt                                                                                    <  2:1575970/150‑1 (MQ=255)

                                                               cGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGAt                                                                                    >  2:1139240/1‑150 (MQ=255)

                                                                 ttATGTCACTGG‑GGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATAcc                                                                                 <  1:1258101/150‑1 (MQ=255)

                                                                   aTGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCa                                                                                <  1:1692501/150‑1 (MQ=255)

                                                                    tGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCAt                                                                               >  2:527220/1‑150 (MQ=255)

                                                                    tGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCAt                                                                               <  1:313528/150‑1 (MQ=255)

                                                                      tCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCa                                                                             <  1:1660100/150‑1 (MQ=255)

                                                                         cTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGc                                                                          >  2:1568928/1‑150 (MQ=255)

                                                                            gCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTgg                                                                       <  2:1005170/150‑1 (MQ=255)

                                                                                 aataatCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAg                                                                  <  2:640132/150‑1 (MQ=255)

                                                                                  ataatCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGa                                                                 >  1:1665263/1‑150 (MQ=255)

                                                                                   taatCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGaa                                                                >  1:719508/1‑150 (MQ=255)

                                                                                      tCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCGAGCTGGCCAAGAAAGt                                                             <  1:1262747/150‑1 (MQ=255)

                                                                                        aGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTa                                                           <  2:540167/150‑1 (MQ=255)

                                                                                        aGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTa                                                           <  2:626164/150‑1 (MQ=255)

                                                                                          tGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGACGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTGATc                                                         >  2:1120314/1‑150 (MQ=255)

                                                                                           ggCACACTGTAGCCTCCCTGCACATGACCATAATGATGTACAGGATTGGCAGTGTCGCTGACCGAGTGCCCAAGGCCTGAAAAGTATAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATcc                                                        <  2:1511916/150‑1 (MQ=255)

                                                                                            gCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGt                                                             >  2:50031/1‑144 (MQ=255)

                                                                                            gCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGt                                                             <  1:50031/144‑1 (MQ=255)

                                                                                            gCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCa                                                       <  2:719508/150‑1 (MQ=255)

                                                                                               cacTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAATTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGtt                                                    >  2:136048/1‑150 (MQ=255)

                                                                                               cacTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGtt                                                    <  2:622638/150‑1 (MQ=255)

                                                                                                 cTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTa                                                  >  1:981548/1‑150 (MQ=255)

                                                                                                   gTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACt                                                <  1:1139240/150‑1 (MQ=255)

                                                                                                     aGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTATAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCCGTGTACTGa                                              >  2:1245284/1‑150 (MQ=255)

                                                                                                     aGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGa                                              <  2:946227/150‑1 (MQ=255)

                                                                                                      gCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGATAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGTAGGATACCATCAAGCTGGCCTAGAAAGTTATCCAGTTTACTGAt                                             <  2:747180/150‑1 (MQ=255)

                                                                                                      gCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGAt                                             <  1:1355981/150‑1 (MQ=255)

                                                                                                      gCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGAt                                             >  2:640110/1‑150 (MQ=255)

                                                                                                          cccTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTg                                         >  2:759465/1‑150 (MQ=255)

                                                                                                          cccTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTg                                         >  2:1284272/1‑150 (MQ=255)

                                                                                                          cccTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGATGACCGACAGCCCAACGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAACGTTATCCAGTTTACTGATGCTg                                         >  2:120548/1‑150 (MQ=255)

                                                                                                            cTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGc                                       >  2:91277/1‑150 (MQ=255)

                                                                                                              tcacaTTACCATAATGATGAACAGGTTTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAATAGTAAAGCATGGCGAAATGACGTGAGTACCGGAGAAGGATACCATCAAGCTGGCCAAGAATGTTATCCAGTTTACTGATGCTGGCGa                                     <  1:189856/149‑1 (MQ=255)

                                                                                                               cacaTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAg                                    <  2:1370269/150‑1 (MQ=255)

                                                                                                                acaTGACCATAATGATGATCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAATTTATCCAGTTTACTGATGCTGGCGAgg                                   >  2:230530/1‑150 (MQ=255)

                                                                                                                acaTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGCgg                                   >  2:235568/1‑150 (MQ=255)

                                                                                                                acaTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAgg                                   >  1:782983/1‑150 (MQ=255)

                                                                                                                  aTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTa                                 <  1:1688666/150‑1 (MQ=255)

                                                                                                                    gACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTAATCAAGATTACTGATGCTGGCGAGGTaaa                               >  2:1382711/1‑149 (MQ=255)

                                                                                                                            tgatgaACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGGTTACTGCTGCGGGCGAGGTACCAGGCAAtt                       >  1:209401/1‑150 (MQ=255)

                                                                                                                                  aCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGACAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCgg                 <  2:1059790/150‑1 (MQ=255)

                                                                                                                                  aCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTAAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCgg                 <  1:143649/150‑1 (MQ=255)

                                                                                                                                       ccttgggcTGTCGGTCAGCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATa            >  2:1633086/8‑150 (MQ=255)

                                                                                                                                        ttGGCACTGTCGCTGCCCGACAGCCCACGGCCAGAAAAGTAACGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCAGCACGGGGGCCTAGAAATTTATCCAGTTGACGGATGGTGGCGCGGTACCCGGCATTTTGGCggggggcg           >  2:848094/1‑146 (MQ=255)

                                                                                                                                        ttGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATAc           >  1:1352714/1‑150 (MQ=255)

                                                                                                                                          ggCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTg         <  2:1014286/150‑1 (MQ=255)

                                                                                                                                          ggCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTg         >  2:1694457/1‑150 (MQ=255)

                                                                                                                                          ggCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTg         >  2:1275466/1‑150 (MQ=255)

                                                                                                                                           gCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGc        <  1:733431/150‑1 (MQ=255)

                                                                                                                                           gCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGc        <  2:1155635/150‑1 (MQ=255)

                                                                                                                                                 tCGCTGACCGACATCCCAAGGCCAGAAAAGTAAAGCATGGCGAACTGACGTGAGTAAGGGCGAAGGATACCATCAAGCTGGCCAAGAAATTTATCCAGTTTACTGATGCTGGCCAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAgg  <  2:244435/150‑1 (MQ=255)

CATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGG  >  NC_000913/4499737‑4500031

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