BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	327,213	G→C	5.9%	S12* (TCA→TGA)	betA ←	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	327,213	0	G	C	5.9%	581.3 / 19.0	170	S12* (TCA→TGA)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (89/70); new base C (4/6); total (93/77)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 174 total aligned reads displayed.

AGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAGTG  >  NC_000913/327064‑327362

aGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTg                                                                                                                                                       <  2:1978345/150‑1 (MQ=255)

   tGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGc                                                                                                                                                    >  1:1007424/1‑150 (MQ=255)

    gTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGcc                                                                                                                                                   <  2:1691606/150‑1 (MQ=255)

    gTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGcc                                                                                                                                                   <  1:752206/150‑1 (MQ=255)

      aGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGATAACGTTGCCGGCTGAGCCgg                                                                                                                                                 <  2:2277965/150‑1 (MQ=255)

      aGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCgg                                                                                                                                                 <  2:132674/150‑1 (MQ=255)

        ggTTTACCCTGTAGCGGGAATGCCAGTGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACTGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCa                                                                                                                                               <  2:1595103/149‑1 (MQ=255)

        cGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCa                                                                                                                                               <  1:544759/150‑1 (MQ=255)

          tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAcc                                                                                                                                             <  1:414848/150‑1 (MQ=255)

           ttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaacg                                                                                                                                            <  2:1104772/150‑4 (MQ=255)

            tACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaa                                                                                                                                           >  1:2196551/1‑150 (MQ=255)

            tACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaacgt                                                                                                                                           <  2:1819301/150‑5 (MQ=255)

            tACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaacgt                                                                                                                                           <  1:2183011/150‑5 (MQ=255)

              cccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaata                                                                                                                                         >  2:2175797/1‑150 (MQ=255)

               cctcatccGGGAATGCCCGGGCAGCGGTGATCTGCGGTCGGCTCTCACCTCGACCGGCCGGGACGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataa                                                                                                                                        <  1:2336012/143‑1 (MQ=255)

                    aGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAatgatg                                                                                                                                   >  1:933178/1‑150 (MQ=255)

                    aGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAatgatg                                                                                                                                   <  1:538248/150‑1 (MQ=255)

                     gCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGt                                                                                                                                  <  1:1958600/150‑1 (MQ=255)

                       gggAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAg                                                                                                                                <  2:20426/150‑1 (MQ=255)

                       gggAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAg                                                                                                                                >  1:2388304/1‑150 (MQ=255)

                          aaTGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCa                                                                                                                             <  2:1685445/150‑1 (MQ=255)

                            tGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCaaa                                                                                                                           >  1:295305/1‑150 (MQ=255)

                             gCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAt                                                                                                                          >  1:478656/1‑150 (MQ=255)

                               cAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTg                                                                                                                        >  1:852526/1‑150 (MQ=255)

                               cAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTg                                                                                                                        <  1:1466802/150‑1 (MQ=255)

                               cAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGAAAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTg                                                                                                                        >  1:1425608/1‑150 (MQ=255)

                                aGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGc                                                                                                                       >  2:451324/1‑150 (MQ=255)

                                 gggCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCa                                                                                                                      <  1:2374394/150‑1 (MQ=255)

                                 gggCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCa                                                                                                                      <  1:136992/150‑1 (MQ=255)

                                 gggCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCa                                                                                                                      <  2:505323/150‑1 (MQ=255)

                                      gCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCAGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCAGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATaa                                                                                                                 >  2:1926104/1‑150 (MQ=255)

                                       cGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATaaa                                                                                                                <  1:1602568/150‑1 (MQ=255)

                                            aTCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAAcctcc                                                                                                           <  2:613717/150‑1 (MQ=255)

                                              cTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTg                                                                                                         <  2:1790712/150‑1 (MQ=255)

                                                 ggtgttgAAGTCAAAGAGTTAGTCCGGGCCGCCCGCTTCAAGCAGCAGGACGGAGGTATTTGTATTTTGAGTGATACGGGTTGTGTTTATGTTTTCGGGTTAGCTTGCACCAATAATGATGTACTCTTTTTGCAATTATAGCTGCTGGtt                                                                                                      <  2:2236908/144‑1 (MQ=255)

                                                   tGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAg                                                                                                    <  1:2316158/150‑1 (MQ=255)

                                                    gCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGa                                                                                                   >  2:1012077/1‑150 (MQ=255)

                                                    gCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGa                                                                                                   >  1:1827158/1‑150 (MQ=255)

                                                     cGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGaa                                                                                                  >  1:1697067/1‑150 (MQ=255)

                                                       gAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGTGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAata                                                                                                >  2:496170/1‑150 (MQ=255)

                                                       gAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGCGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAata                                                                                                >  2:2132693/1‑150 (MQ=255)

                                                        aaGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAatat                                                                                               <  2:1166184/150‑1 (MQ=255)

                                                         aGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATg                                                                                              <  1:2084040/150‑1 (MQ=255)

                                                            cAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGTGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGAc                                                                                           >  2:1671769/1‑150 (MQ=255)

                                                            cAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGATGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGTCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGAc                                                                                           <  2:864995/150‑1 (MQ=255)

                                                            cAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGAc                                                                                           <  1:1920452/150‑1 (MQ=255)

                                                            cAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGAc                                                                                           >  2:466308/1‑150 (MQ=255)

                                                             aaaGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGTGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACt                                                                                          >  1:907632/1‑150 (MQ=255)

                                                               aGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTgg                                                                                        >  2:783400/1‑150 (MQ=255)

                                                                 cGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGt                                                                                                       >  2:2084768/1‑133 (MQ=255)

                                                                    tAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAAttt                                                                                   <  2:2189007/150‑1 (MQ=255)

                                                                    tAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAAttt                                                                                   >  1:1774974/1‑150 (MQ=255)

                                                                      gTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAg                                                                                 <  1:734092/150‑1 (MQ=255)

                                                                      gTCCGGGCCGACCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGCACGTTGCCGGCTGCGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAg                                                                                 >  2:539627/1‑150 (MQ=255)

                                                                       tCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGc                                                                                <  1:613152/150‑1 (MQ=255)

                                                                        ccGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGcc                                                                               <  2:2196551/150‑1 (MQ=255)

                                                                        ccGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGATAGAATATGGACTGGAATTTAGcc                                                                               <  1:1321355/150‑1 (MQ=255)

                                                                         cGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCa                                                                              <  1:2175797/150‑1 (MQ=255)

                                                                              cgcccgcTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCa                                                                         >  2:762708/1‑150 (MQ=255)

                                                                                   gcTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtg                                                                    >  1:1939685/1‑150 (MQ=255)

                                                                                     ttCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtgga                                                                  <  2:298243/150‑1 (MQ=255)

                                                                                     ttCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtgga                                                                  <  2:569429/150‑1 (MQ=255)

                                                                                      tCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggat                                                                 >  2:498729/1‑150 (MQ=255)

                                                                                      tCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggat                                                                 >  1:1120846/1‑150 (MQ=255)

                                                                                      tCAAGCAGCAGCACGGAGGTATTCGGAGCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggat                                                                 <  2:2172999/150‑1 (MQ=255)

                                                                                       cAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatg                                                                <  2:2388304/150‑1 (MQ=255)

                                                                                       cAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatg                                                                <  1:1102993/150‑1 (MQ=255)

                                                                                        aaGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGACTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatgg                                                               >  1:890389/1‑150 (MQ=255)

                                                                                         agcagcagcACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatgga                                                              >  1:120311/1‑150 (MQ=255)

                                                                                         agcagcagcACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatgga                                                              <  2:1156869/150‑1 (MQ=255)

                                                                                         agcagcagcACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatgga                                                              >  1:2088170/1‑150 (MQ=255)

                                                                                          gcagcagcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGAc                                                             >  1:873224/1‑150 (MQ=255)

                                                                                               agcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCAc                                                        >  2:594971/1‑150 (MQ=255)

                                                                                                gcaCGGAGGTATTCGGATATTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAGTGCAAATAAACCTCATGGTTAGAATATGGACTGGAATTTAGCCATCTAAACCTGGATGGACTTAAcc                                                       >  1:816091/1‑150 (MQ=255)

                                                                                                  aCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTg                                                     <  2:1425608/150‑1 (MQ=255)

                                                                                                    ggAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTggg                                                   <  1:1012077/150‑1 (MQ=255)

                                                                                                    ggAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTggg                                                   <  1:2462430/150‑1 (MQ=255)

                                                                                                        gTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTa                                               <  1:1037666/150‑1 (MQ=255)

                                                                                                        gTATTCGGATCTTCAGTCAGACGGGGAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATTAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTa                                               <  2:2348729/150‑1 (MQ=255)

                                                                                                         tATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTaa                                              >  1:521773/1‑150 (MQ=255)

                                                                                                            tCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAActc                                           >  2:351782/1‑150 (MQ=255)

                                                                                                               gATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTgg                                        <  1:1345547/150‑1 (MQ=255)

                                                                                                                  cTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGc                                     >  2:765317/1‑150 (MQ=255)

                                                                                                                  cTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGc                                     >  1:98442/1‑150 (MQ=255)

                                                                                                                  cTTCAGTCAGACGCGTAGCGAGAACGTTGCCGGCAGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATCGACTGGAATTTAGCCATCTCAACCTCGATGGACTTCACCTGGCAGTAACTCTAGAGc                                     >  1:1631350/1‑150 (MQ=255)

                                                                                                                       gtcagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtcat                                >  2:1746845/1‑150 (MQ=255)

                                                                                                                             cGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGcc                          >  2:1443493/1‑150 (MQ=255)

                                                                                                                                  aGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTtct                     >  2:886124/1‑150 (MQ=255)

                                                                                                                                    cGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATaaa                                                                                                                <  2:1033411/57‑1 (MQ=255)

                                                                                                                                    cGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTcg                   <  1:1800002/150‑1 (MQ=255)

                                                                                                                                        aaCGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGAc               <  2:242061/150‑1 (MQ=255)

                                                                                                                                         aCGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGAcc              <  2:369491/150‑1 (MQ=255)

                                                                                                                                         aCGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGAcc              <  2:2166032/150‑1 (MQ=255)

                                                                                                                                          cGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCa             >  2:1901168/1‑150 (MQ=255)

                                                                                                                                            ttGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAAt           >  1:335132/1‑150 (MQ=255)

                                                                                                                                            gtGCCGGCTCAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAAt           >  2:1820259/2‑150 (MQ=255)

                                                                                                                                               ccGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGcc        >  2:1537206/1‑150 (MQ=255)

                                                                                                                                                  gCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGa     >  1:1068814/1‑150 (MQ=255)

                                                                                                                                                     gAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAgtg  <  1:451324/150‑1 (MQ=255)

                                                                                                                                                     gAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAgtg  <  1:451263/150‑1 (MQ=255)

AGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAGTG  >  NC_000913/327064‑327362

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