BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,986,123	T→C	Y93H (TAC→CAC)	yqeG →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,986,123	0	T	C	97.5%	435.3 / 1.0	118	Y93H (TAC→CAC)	yqeG	putative transporter
Reads supporting (aligned to +/- strand): ref base T (3/0); new base C (43/72); total (46/72)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 118 total aligned reads displayed.

GCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGTACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGAATCTCATTT  >  NC_000913/2985976‑2986270

gCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATGAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACt                                                                                                                                                  <  1:485605/151‑1 (MQ=255)

gCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACt                                                                                                                                                  <  1:1273685/151‑1 (MQ=255)

    tGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCAt                                                                                                                                              >  1:1065585/1‑151 (MQ=255)

     ggTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTc                                                                                                                                                <  1:1631605/148‑1 (MQ=255)

     ggTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTc                                                                                                                                                <  1:142154/148‑1 (MQ=255)

     ggTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCacc                                                                                                                                             <  2:1039443/151‑3 (MQ=255)

     ggTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCa                                                                                                                                               <  2:1364961/149‑1 (MQ=255)

     ggTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATc                                                                                                                                             <  2:662671/151‑1 (MQ=255)

     ggTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGCATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCAt                                                                                                                                              <  2:352087/150‑1 (MQ=255)

        cGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGcc                                                                                                                                          >  1:1514575/1‑151 (MQ=255)

        cGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGcc                                                                                                                                          <  1:175760/151‑1 (MQ=255)

          ccTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCtt                                                                                                                                        <  1:70156/151‑1 (MQ=255)

             ggCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCttt                                                                                                                                       <  1:486208/149‑1 (MQ=255)

              gCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCtttttt                                                                                                                                    <  1:30040/151‑1 (MQ=255)

               ccTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTg                                                                                                                                   >  1:1258321/1‑151 (MQ=255)

                    aaCATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTgtcgtc                                                                                                                              <  2:576262/151‑1 (MQ=255)

                    aaCATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTgtcgtc                                                                                                                              <  2:1125517/151‑1 (MQ=255)

                    aaCATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTgtcgtc                                                                                                                              >  1:12739/1‑151 (MQ=255)

                      cATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTcgtcgt                                                                                                                            <  2:474008/151‑1 (MQ=255)

                      cATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTcgtcgt                                                                                                                            >  2:1655671/1‑151 (MQ=255)

                       atatTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGTACTTCATCGCCTTTTTTGTCGTCgtg                                                                                                                           >  2:318002/1‑151 (MQ=255)

                         atTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCgtgt                                                                                                                          <  2:852845/150‑1 (MQ=255)

                         atTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGtt                                                                                                                         >  2:683207/1‑151 (MQ=255)

                           tGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTg                                                                                                                        >  1:815383/1‑150 (MQ=255)

                             gCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGata                                                                                                                     <  1:1479216/151‑1 (MQ=255)

                                acaTAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTTCTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGatatat                                                                                                                  <  1:957875/151‑1 (MQ=255)

                                acaTAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGatatat                                                                                                                  <  1:507728/151‑1 (MQ=255)

                                    aaaGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGc                                                                                                                >  2:428621/1‑149 (MQ=255)

                                    aaaGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCa                                                                                                               <  2:1287348/150‑1 (MQ=255)

                                     aaGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGt                                                                                                             <  1:1336770/151‑1 (MQ=255)

                                      aGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTg                                                                                                            <  1:133756/151‑1 (MQ=255)

                                       gCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTgg                                                                                                           <  1:941150/151‑1 (MQ=255)

                                        ccTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGc                                                                                                          <  1:1696595/151‑1 (MQ=255)

                                          ttGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGTACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCa                                                                                                         >  2:620252/1‑150 (MQ=255)

                                           tgtgCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCaa                                                                                                        <  1:1312740/150‑1 (MQ=255)

                                           tgtgCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCa                                                                                                         <  2:1670623/149‑1 (MQ=255)

                                            gtgCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAAtt                                                                                                      >  1:1267616/1‑151 (MQ=255)

                                            gtgCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAAtt                                                                                                      <  2:129944/151‑1 (MQ=255)

                                              gCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTAc                                                                                                    <  1:9801/151‑1 (MQ=255)

                                              gCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTAc                                                                                                    >  1:882543/1‑151 (MQ=255)

                                                 aGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCa                                                                                                  <  2:1497079/150‑1 (MQ=255)

                                                       tCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACt                                                                                           >  2:346468/1‑151 (MQ=255)

                                                        ccTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACtt                                                                                          >  2:1545033/1‑151 (MQ=255)

                                                        ccTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACtt                                                                                          <  2:1287196/151‑1 (MQ=255)

                                                        ccTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACtt                                                                                          >  1:834762/1‑151 (MQ=255)

                                                          tctcATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTAc                                                                                        >  1:318441/1‑151 (MQ=255)

                                                            tcatcaAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACgg                                                                                      >  2:150948/1‑151 (MQ=255)

                                                                caAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAAc                                                                                   >  1:579275/1‑150 (MQ=255)

                                                                caAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACa                                                                                  <  2:813237/151‑1 (MQ=255)

                                                                  aaaaCATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGc                                                                                >  2:17281/1‑151 (MQ=255)

                                                                  aaaaCATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGc                                                                                <  2:978567/151‑1 (MQ=255)

                                                                   aaaCATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGTACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCt                                                                               >  1:1053650/1‑151 (MQ=255)

                                                                   aaaCATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCt                                                                               >  1:1663564/1‑151 (MQ=255)

                                                                    aaCATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTg                                                                              <  1:1375097/151‑1 (MQ=255)

                                                                     aCATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTgg                                                                             <  2:547134/151‑1 (MQ=255)

                                                                       aTCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGc                                                                            >  2:332634/1‑150 (MQ=255)

                                                                          agcagGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCaa                                                                          >  2:140309/1‑149 (MQ=255)

                                                                          agcagGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAg                                                                        <  2:1403588/151‑1 (MQ=255)

                                                                             agGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCat                                                                     >  2:1096513/1‑151 (MQ=255)

                                                                             agGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCat                                                                     >  2:911810/1‑151 (MQ=255)

                                                                               gTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCatat                                                                   <  1:780717/151‑1 (MQ=255)

                                                                                 gAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATgg                                                                 <  1:1027486/151‑1 (MQ=255)

                                                                                 gAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATgg                                                                 <  1:1418783/151‑1 (MQ=255)

                                                                                 gAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATg                                                                  <  2:105243/150‑1 (MQ=255)

                                                                                  aaGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGt                                                                >  1:1621340/1‑151 (MQ=255)

                                                                                  aaGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGt                                                                >  2:355416/1‑151 (MQ=255)

                                                                                  aaGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGt                                                                >  2:1174311/1‑151 (MQ=255)

                                                                                  aaGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGt                                                                >  1:556173/1‑151 (MQ=255)

                                                                                   aGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGtt                                                               >  2:1451969/1‑151 (MQ=255)

                                                                                      gATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTAtt                                                            <  2:1684608/151‑1 (MQ=255)

                                                                                             ggCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTc                                                     <  2:1329167/151‑1 (MQ=255)

                                                                                             ggCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTc                                                     <  2:1621439/151‑1 (MQ=255)

                                                                                               cgcgGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGc                                                   <  2:278149/151‑1 (MQ=255)

                                                                                                       cacacTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTAt                                           >  2:1693110/1‑151 (MQ=255)

                                                                                                        acacTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATg                                          >  2:130512/1‑151 (MQ=255)

                                                                                                        acacTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATg                                          >  1:972024/1‑151 (MQ=255)

                                                                                                             actaTGGCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTgg                                      <  1:1096428/150‑1 (MQ=255)

                                                                                                                  ggCAAGAAGATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTc                                <  2:178872/151‑1 (MQ=255)

                                                                                                                     aagaagATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTgg                             <  1:332608/151‑1 (MQ=255)

                                                                                                                     aagaagATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTgg                             <  1:1279286/151‑1 (MQ=255)

                                                                                                                     aagaagATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTgg                             >  1:787600/1‑151 (MQ=255)

                                                                                                                     aagaagATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTgg                             <  1:1479159/151‑1 (MQ=255)

                                                                                                                     aagaagATTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTgg                             <  2:1065668/151‑1 (MQ=255)

                                                                                                                      agaagaTTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTggg                            <  2:1332117/151‑1 (MQ=255)

                                                                                                                       gaagaTTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGgt                           >  2:601644/1‑151 (MQ=255)

                                                                                                                        aagaTTGGTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGgt                           <  1:1639534/150‑1 (MQ=255)

                                                                                                                              ggTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTAt                    <  2:469548/151‑1 (MQ=255)

                                                                                                                               gTAATCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTAtt                   >  1:822111/1‑151 (MQ=255)

                                                                                                                                 aaTCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTaa                 <  2:556223/151‑1 (MQ=255)

                                                                                                                                 aaTCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTa                  <  2:9449/150‑1 (MQ=255)

                                                                                                                                 aaTCTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTa                  >  1:1109879/1‑150 (MQ=255)

                                                                                                                                    cTGATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTc              >  1:591536/1‑151 (MQ=255)

                                                                                                                                      gATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTc              <  1:817419/149‑1 (MQ=255)

                                                                                                                                      gATTACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTc              <  2:1458736/149‑1 (MQ=255)

                                                                                                                                        ttACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGa           >  2:594753/1‑150 (MQ=255)

                                                                                                                                         tACCACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGaa          <  2:1207127/150‑1 (MQ=255)

                                                                                                                                           ccACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGAAtc        <  1:1655563/150‑1 (MQ=255)

                                                                                                                                            cACGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGAAtctc      >  1:628507/1‑151 (MQ=255)

                                                                                                                                              cGCTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGAATCTCAt    <  1:1530200/151‑1 (MQ=255)

                                                                                                                                                cTGCACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGAATCTCAttt  >  1:1248569/1‑151 (MQ=255)

GCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGGATAACGGGCGCGGTAACACACTACTATGGCAAGAAGATTGGTAATCTGATTACCACGCTGTACTTCATCGCCTTTTTTGTCGTCGTGTTGATATATGCAGTGGCAATTACCAACTCACTTACGGAACAGCTGGCAAAGCATATGGTTATTGATCTTCGCATCCGTATGTTGGTGAGTCTGGGTGTTGTATTAATTCTGAATCTCATTT  >  NC_000913/2985976‑2986270

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