BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,859,354	T→C	35.6%	N255D (AAC→GAC)	yidK ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,859,354	0	T	C	35.6%	449.8 / inf	434	N255D (AAC→GAC)	yidK	putative transporter
Reads supporting (aligned to +/- strand): ref base T (144/135); new base C (91/63); total (236/198)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 437 total aligned reads displayed.

AACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATA  >  NC_000913/3859206‑3859500

aaCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATGCCCCTGATTTGTACACCAGTAAAAGGTGtt                                                                                                                                                    >  2:1230624/1‑150 (MQ=255)

aaCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTc                                                                                                                                                   >  2:851973/1‑150 (MQ=255)

 aCGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCa                                                                                                                                                  <  2:2948806/150‑1 (MQ=255)

  cGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCAc                                                                                                                                                 >  2:1065403/1‑150 (MQ=255)

   gCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCAcc                                                                                                                                                <  2:730907/150‑1 (MQ=255)

   gCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCAcc                                                                                                                                                >  2:763256/1‑150 (MQ=255)

   gCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCAcc                                                                                                                                                <  2:2257053/150‑1 (MQ=255)

    cAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCa                                                                                                                                               <  1:1145475/150‑1 (MQ=255)

     aaTCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAg                                                                                                                                              >  1:131160/1‑150 (MQ=255)

       tCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGaa                                                                                                                                            >  2:887277/1‑150 (MQ=255)

       tCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGaa                                                                                                                                            <  1:3079023/150‑1 (MQ=255)

       tCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGaa                                                                                                                                            >  2:2747694/1‑150 (MQ=255)

           cccTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCaa                                                                                                                                        >  1:2541700/1‑150 (MQ=255)

             cTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAAc                                                                                                                                      >  2:1132995/1‑150 (MQ=255)

                gcagcaCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCgg                                                                                                                                   <  2:1911901/150‑1 (MQ=255)

                  agcaCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTa                                                                                                                                 >  2:1638150/1‑150 (MQ=255)

                  agcaCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTa                                                                                                                                 >  2:581848/1‑150 (MQ=255)

                    caCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTaaa                                                                                                                               >  2:3087015/1‑150 (MQ=255)

                     aCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAAt                                                                                                                              >  1:2402303/1‑150 (MQ=255)

                     aCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAAt                                                                                                                              <  1:1733639/150‑1 (MQ=255)

                     aCCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAAt                                                                                                                              <  2:1712611/150‑1 (MQ=255)

                      ccAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATg                                                                                                                             >  1:1198531/1‑150 (MQ=255)

                       cAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGc                                                                                                                            >  1:564963/1‑150 (MQ=255)

                       cAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGc                                                                                                                            >  1:446785/1‑150 (MQ=255)

                       cAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGc                                                                                                                            >  1:1637182/1‑150 (MQ=255)

                         gTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCcg                                                                                                                          <  2:2624773/150‑1 (MQ=255)

                            ccAGCGGGTCGAGCATTTTCAGCACCGCCGTTATCAGCGCCCCCTTTTGCGCTTCCGCCAGGCTTTTTGACGCGAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAATGGTGTCCACCAGAATCAAACCGGTAAATGCCgcgc                                                                                                                       <  1:1900206/150‑1 (MQ=255)

                             cAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGcc                                                                                                                      >  1:2133820/1‑150 (MQ=255)

                               gCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCCGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCaa                                                                                                                    >  1:511445/1‑150 (MQ=255)

                               gCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCaa                                                                                                                    >  2:3016030/1‑150 (MQ=255)

                               gCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCaa                                                                                                                    <  2:2817284/150‑1 (MQ=255)

                                 gggTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATc                                                                                                                  >  1:2720113/1‑150 (MQ=255)

                                  ggTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTACATGCCGCGCCAATCg                                                                                                                 <  1:2787711/150‑1 (MQ=255)

                                   gTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCgg                                                                                                                >  2:1958573/1‑150 (MQ=255)

                                      gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                             <  2:151379/150‑1 (MQ=255)

                                      gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                             <  1:610630/150‑1 (MQ=255)

                                      gAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCaa                                                                                                             <  1:1077254/150‑1 (MQ=255)

                                           ttttCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAt                                                                                                        <  1:1246173/150‑1 (MQ=255)

                                            tttCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtc                                                                                                       >  1:209298/1‑150 (MQ=255)

                                            tttCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtc                                                                                                       >  2:1363062/1‑150 (MQ=255)

                                              tCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCGTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcgg                                                                                                     >  2:2074109/1‑150 (MQ=255)

                                                 gCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGAtcggtcg                                                                                                  >  1:2193221/1‑150 (MQ=255)

                                                     cgccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGcc                                                                                              >  1:2880502/1‑150 (MQ=255)

                                                      gccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACTCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAATAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCAAATAGGCAAGGGATCGGTCGGGCCa                                                                                             <  2:356641/150‑1 (MQ=255)

                                                      gccgTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCa                                                                                             <  2:1997808/150‑1 (MQ=255)

                                                            aaCAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGAtt                                                                                       <  2:1568869/150‑1 (MQ=255)

                                                              cAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGa                                                                                     <  1:1195929/150‑1 (MQ=255)

                                                               aGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAg                                                                                    >  1:66909/1‑150 (MQ=255)

                                                               aGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAg                                                                                    <  1:2857955/150‑1 (MQ=255)

                                                                   cccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTa                                                                                >  2:2118634/1‑150 (MQ=255)

                                                                    ccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTCAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGGCGGGCCACCGATTGAGTTTaa                                                                               >  2:1425357/1‑150 (MQ=255)

                                                                    ccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTaa                                                                               <  1:3111250/150‑1 (MQ=255)

                                                                    ccccTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTaa                                                                               >  1:32966/1‑150 (MQ=255)

                                                                         tttGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTtct                                                                          <  1:1603532/150‑1 (MQ=255)

                                                                           tGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCg                                                                        <  2:645382/150‑1 (MQ=255)

                                                                            gCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCgg                                                                       <  2:1679954/150‑1 (MQ=255)

                                                                            gCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCgg                                                                       >  2:536191/1‑150 (MQ=255)

                                                                               cTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGt                                                                    <  1:1827600/150‑1 (MQ=255)

                                                                               cTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGt                                                                    <  2:119220/150‑1 (MQ=255)

                                                                               cTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGt                                                                    >  1:842748/1‑150 (MQ=255)

                                                                                  ccGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGaa                                                                 <  1:1444489/150‑1 (MQ=255)

                                                                                  ccGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGaa                                                                 <  1:1910005/150‑1 (MQ=255)

                                                                                      cAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggt                                                             >  2:2923096/1‑150 (MQ=255)

                                                                                       aGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtg                                                            <  1:3018436/150‑1 (MQ=255)

                                                                                       aGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtg                                                            >  2:1869165/1‑150 (MQ=255)

                                                                                       aGGCTTGTTGCCGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACCCCAGTAACAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtg                                                            <  1:1353469/150‑1 (MQ=255)

                                                                                         gCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtggt                                                          >  2:2543300/1‑150 (MQ=255)

                                                                                         gCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACggtggt                                                          >  1:1993866/1‑150 (MQ=255)

                                                                                           tttttGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGa                                                        >  1:2335131/1‑150 (MQ=255)

                                                                                           tttttGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGa                                                        >  2:843299/1‑150 (MQ=255)

                                                                                           tttttGACGCCAGCGAGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACCGTGGTGa                                                        >  1:2368361/1‑150 (MQ=255)

                                                                                              ttGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGtt                                                     <  2:1164508/150‑1 (MQ=255)

                                                                                              ttGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGtt                                                     <  1:2951669/150‑1 (MQ=255)

                                                                                               tGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAACTGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTg                                                    <  2:500629/150‑1 (MQ=255)

                                                                                                   gCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCa                                                >  1:2468138/1‑150 (MQ=255)

                                                                                                      aGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATg                                             >  2:2984802/1‑150 (MQ=255)

                                                                                                         gtgcgtTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGccc                                          >  2:1022457/1‑150 (MQ=255)

                                                                                                         gtgcgtTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGccc                                          <  2:1280364/150‑1 (MQ=255)

                                                                                                             gtTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCa                                      >  2:3008323/1‑150 (MQ=255)

                                                                                                                gCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATaa                                   <  2:2374948/150‑1 (MQ=255)

                                                                                                                  aCGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTACATAAAg                                 <  1:913764/150‑1 (MQ=255)

                                                                                                                  aCGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAg                                 >  2:1831273/1‑150 (MQ=255)

                                                                                                                    gATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCt                               >  1:1107429/1‑150 (MQ=255)

                                                                                                                      tG‑CCCTGATTTGTACACCAGTAAAAGGTGCTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGc                             <  2:2689345/150‑1 (MQ=255)

                                                                                                                           cTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCtt                         >  1:1024919/1‑150 (MQ=255)

                                                                                                                           cTGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCtt                         >  1:1044153/1‑150 (MQ=255)

                                                                                                                            tGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTg                        <  1:2617923/150‑1 (MQ=255)

                                                                                                                            tGATTTGTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTg                        >  2:1395484/1‑150 (MQ=255)

                                                                                                                             gATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGc                       <  2:2579496/150‑1 (MQ=255)

                                                                                                                                  gTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATc                  <  1:794880/150‑1 (MQ=255)

                                                                                                                                  gTACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATc                  <  2:1731631/150‑1 (MQ=255)

                                                                                                                                   tACACCAGTAAAAGGTGTTCACCAGACTCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATcg                 <  1:15367/150‑1 (MQ=255)

                                                                                                                                   tACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATcg                 <  1:3005576/150‑1 (MQ=255)

                                                                                                                                   tACACCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATcg                 <  2:872532/150‑1 (MQ=255)

                                                                                                                                      acCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGa              <  2:401210/150‑1 (MQ=255)

                                                                                                                                      acCAGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGa              >  2:2400092/1‑150 (MQ=255)

                                                                                                                                         aGTAAAAGGTGTCCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATGTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTa           <  1:1160468/150‑1 (MQ=255)

                                                                                                                                          gTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAg          >  1:754942/1‑150 (MQ=255)

                                                                                                                                              aaGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCa      >  1:2429949/1‑150 (MQ=255)

                                                                                                                                                ggTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCaaa    <  1:860905/150‑1 (MQ=255)

AACGCAATCAACCCTGGCAGCACCAGTACCAGCGGGTCGAGCATTTTCAGCACCGCCGTTAACAGCGCCCCCTTTTGCCCTTCCGCCAGGCTTTTTGACGCCAGCGTGCGTTGCACGATG‑CCCTGATTTGTACACCAGTAAAAGGTGTTCACCAGAATCAAACCGGTAAATGCCGCGCCAATCGGCAAGGGATCGGTCGGGCCACCGATTGAGTTTAATTTCTCGGCGTGAACGGTGGTGAGTTGCTCAATGCCCTGCATAAAGCTGCCCTTGCCCATCGCGATTAGGCCAAATA  >  NC_000913/3859206‑3859500

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