BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,307,294	A→T	12.6%	L169Q (CTG→CAG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,294	0	A	T	12.6%	332.5 / 36.1	111	L169Q (CTG→CAG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand): ref base A (46/51); new base T (6/8); total (52/59)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

Only 100 of 113 total aligned reads displayed.

GGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAG  >  NC_000913/4307150‑4307442

ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                  <  4:489461/149‑1 (MQ=255)

ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                  <  4:539735/149‑1 (MQ=255)

  ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                <  2:31796/149‑1 (MQ=255)

  ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGACCAATTGGTATCCAGGCAGCCAGGATTGCCACACGCGCAGTGTTGGGTAATATCTCCCAGGGGGa                                                                                                                                                >  4:551612/1‑149 (MQ=255)

    cGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGGAGTGTTGGGTCATATCTCCCAGGGGGata                                                                                                                                              <  3:551612/149‑1 (MQ=255)

     gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGatat                                                                                                                                             <  4:577593/149‑1 (MQ=255)

     gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGatat                                                                                                                                             <  3:108796/149‑1 (MQ=255)

      ttttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATg                                                                                                                                            >  3:223942/1‑149 (MQ=255)

       tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGa                                                                                                                                           <  1:869103/149‑1 (MQ=255)

       tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                           <  4:724053/149‑1 (MQ=255)

       tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                           <  1:306747/149‑1 (MQ=255)

       tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCAAATCCCCCTGGGAGATATGa                                                                                                                                           <  3:765229/149‑1 (MQ=255)

         tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                         <  4:66395/149‑1 (MQ=255)

          ccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAccc                                                                                                                                        <  1:714583/149‑1 (MQ=255)

             cATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                     <  4:794687/149‑1 (MQ=255)

             cATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                                                                                     >  2:5182/1‑149 (MQ=255)

                  aCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                <  1:5182/149‑1 (MQ=255)

                   cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCt                                                                                                                               <  3:459411/149‑1 (MQ=255)

                    gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                              >  4:564297/1‑149 (MQ=255)

                    gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                              >  4:388467/1‑149 (MQ=255)

                    gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAccc                                                                                                                                        <  2:477807/139‑1 (MQ=255)

                    gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAccc                                                                                                                                        >  1:477807/1‑139 (MQ=255)

                     aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGc                                                                                                                             >  2:529653/1‑149 (MQ=255)

                      aaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                            >  4:645450/1‑149 (MQ=255)

                         gATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCcaca                                                                                                                         >  4:500644/1‑149 (MQ=255)

                             gcgcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgt                                                                                                                     >  1:242729/1‑149 (MQ=255)

                               gcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtgt                                                                                                                   <  2:606178/149‑1 (MQ=255)

                               gcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtgt                                                                                                                   <  4:88294/149‑1 (MQ=255)

                                cAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGtgtg                                                                                                                  >  3:470692/1‑149 (MQ=255)

                                     ggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAccc                                                                                                             <  3:23534/149‑1 (MQ=255)

                                     ggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAccc                                                                                                             >  2:243272/1‑149 (MQ=255)

                                              ggggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgc                                                                                                    <  4:883848/149‑1 (MQ=255)

                                                 gCTGTTGTTCGTATCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCa                                                                                                 >  1:58616/1‑149 (MQ=255)

                                                 gCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCa                                                                                                 >  4:655329/1‑149 (MQ=255)

                                                 gCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCa                                                                                                 >  1:205295/1‑149 (MQ=255)

                                                     ttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtgt                                                                                                                   >  2:625795/1‑127 (MQ=255)

                                                     ttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtgt                                                                                                                   <  1:625795/127‑1 (MQ=255)

                                                     ttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATg                                                                                                                                            <  4:254541/102‑1 (MQ=255)

                                                     tngttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATg                                                                                                                                            >  3:254541/1‑102 (MQ=255)

                                                       gttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                           <  3:655329/149‑1 (MQ=255)

                                                        ttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                          <  1:584837/149‑1 (MQ=255)

                                                          cGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACTGTGTGCACCCGTCCACGGCGCACCGTTCATc                                                                                        <  2:595922/149‑1 (MQ=255)

                                                           gTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc                                                                                       <  2:289541/149‑1 (MQ=255)

                                                              ccAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcaca                                                                                    <  2:622713/149‑1 (MQ=255)

                                                               cAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                   >  1:790390/1‑149 (MQ=255)

                                                                     cTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGaa                                                                             >  2:665914/1‑149 (MQ=255)

                                                                         gcgcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAcccc                                                                         >  4:500153/1‑149 (MQ=255)

                                                                            ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                                                                      <  4:841794/149‑1 (MQ=255)

                                                                            ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc                                                                      >  1:125159/1‑149 (MQ=255)

                                                                                tCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCg                                                                  >  4:41478/1‑149 (MQ=255)

                                                                                tCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCg                                                                  >  3:790387/1‑149 (MQ=255)

                                                                                   agagCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGTTCATATCTCCCATGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAc                                                               <  2:779117/149‑1 (MQ=255)

                                                                                      gCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCga                                                            <  2:205295/149‑1 (MQ=255)

                                                                                        aaTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                          >  3:887559/1‑149 (MQ=255)

                                                                                         aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                              >  1:680901/1‑112 (MQ=255)

                                                                                         aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                              <  2:680901/112‑1 (MQ=255)

                                                                                         aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCACCCACACTGCGAACCCCATCCCCGTACCGAGAt                                                         <  2:928114/149‑1 (MQ=255)

                                                                                           tGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCGACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAg                                                       <  4:621696/149‑1 (MQ=255)

                                                                                             gTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgc                                                     <  3:810009/149‑1 (MQ=255)

                                                                                              tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCga                                                            >  4:476944/1‑141 (MQ=255)

                                                                                              tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCga                                                            <  3:476944/141‑1 (MQ=255)

                                                                                              tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgcc                                                    >  2:26089/1‑149 (MQ=255)

                                                                                              tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCGCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACAGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgcc                                                    >  4:568925/1‑149 (MQ=255)

                                                                                                   aGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                     >  4:879150/1‑127 (MQ=255)

                                                                                                   aGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAg                                               >  1:394681/1‑149 (MQ=255)

                                                                                                     gCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                             >  1:888176/1‑149 (MQ=255)

                                                                                                     gCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                             <  1:870770/149‑1 (MQ=255)

                                                                                                        cccAGGATTGCCACACGCGCAGTGTTGGGTCATAGCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCa                                          <  1:243272/149‑1 (MQ=255)

                                                                                                         ccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGTGAACCCTATCCCCGTACCGAGACAGGCCGCCAGAACCAg                                         >  3:527416/1‑149 (MQ=255)

                                                                                                          cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                        >  2:702125/1‑149 (MQ=255)

                                                                                                          cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                        <  1:858030/149‑1 (MQ=255)

                                                                                                           aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgtt                                       <  4:470692/149‑1 (MQ=255)

                                                                                                             gATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgt                                     >  2:163467/1‑149 (MQ=255)

                                                                                                              aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                    <  3:902614/149‑1 (MQ=255)

                                                                                                                tGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGc                                  <  4:424171/149‑1 (MQ=255)

                                                                                                                tGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGc                                  <  3:500153/149‑1 (MQ=255)

                                                                                                                  ccACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                >  2:898171/1‑149 (MQ=255)

                                                                                                                    agacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTaa                              >  2:743147/3‑149 (MQ=255)

                                                                                                                       cgcgCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGc                           <  1:415595/149‑1 (MQ=255)

                                                                                                                           cAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGtt                       <  2:58616/149‑1 (MQ=255)

                                                                                                                               gtTGGGTCATAGCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTct                   >  3:598589/1‑149 (MQ=255)

                                                                                                                                 tGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                 <  4:486843/149‑1 (MQ=255)

                                                                                                                                 tGGGTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                 >  3:62746/1‑149 (MQ=255)

                                                                                                                                   ggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGt                        <  1:388548/140‑1 (MQ=255)

                                                                                                                                     tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg             >  2:175999/1‑149 (MQ=255)

                                                                                                                                     tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg             >  2:305872/1‑149 (MQ=255)

                                                                                                                                     tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg             >  3:712736/1‑149 (MQ=255)

                                                                                                                                     tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg             >  2:379534/1‑149 (MQ=255)

                                                                                                                                        tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCCCACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTcc          >  3:489441/1‑149 (MQ=255)

                                                                                                                                        tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTcc          <  3:100023/149‑1 (MQ=255)

                                                                                                                                        tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCACACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTcc          <  3:500644/149‑1 (MQ=255)

                                                                                                                                         atCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc         <  3:710988/149‑1 (MQ=255)

                                                                                                                                         atCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc         <  2:888176/149‑1 (MQ=255)

                                                                                                                                         atCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc         >  2:834387/1‑149 (MQ=255)

                                                                                                                                           ctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAg       <  4:790387/149‑1 (MQ=255)

                                                                                                                                            tcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAgg      >  2:768389/1‑149 (MQ=255)

                                                                                                                                             cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGga     <  4:634240/149‑1 (MQ=255)

                                                                                                                                              ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgag    >  4:127219/1‑149 (MQ=255)

                                                                                                                                                aGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag  >  1:5545/1‑149 (MQ=255)

                                                                                                                                                aGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag  >  2:359525/1‑149 (MQ=255)

GGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAG  >  NC_000913/4307150‑4307442

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