BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	662,636	T→G	26.3%	T3P (ACC→CCC)	ybeD ←	UPF0250 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	662,636	0	T	G	26.3%	285.8 / 130.3	175	T3P (ACC→CCC)	ybeD	UPF0250 family protein
Reads supporting (aligned to +/- strand): ref base T (75/53); new base G (17/29); total (93/82)
Fisher's exact test for biased strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 183 total aligned reads displayed.

GTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACT  >  NC_000913/662537‑662735

  aCCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGgtt                                                                                                   >  1:3886991/1‑100 (MQ=255)

   ccACTTCGACCACCTGATCAACCAGCTCAGGTAACGCCTGTCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGtttt                                                                                                  >  1:2167253/1‑100 (MQ=255)

   ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGtttt                                                                                                  >  2:2224552/1‑100 (MQ=255)

   ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGtttt                                                                                                  >  1:1411072/1‑100 (MQ=255)

    cACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTc                                                                                                 >  2:3101430/1‑100 (MQ=255)

     aCTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTATGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCa                                                                                                <  2:2477925/100‑1 (MQ=255)

     aCTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCa                                                                                                <  2:3698107/100‑1 (MQ=255)

      cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCAt                                                                                               <  1:2254132/100‑1 (MQ=255)

      cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCAt                                                                                               >  2:366700/1‑100 (MQ=255)

      cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCAt                                                                                               <  1:2395372/100‑1 (MQ=255)

        tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGt                                                                                             <  1:475224/100‑1 (MQ=255)

        tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTAGTAAGTACAAGGCGTACGGACTTCAAACTGTTCGATAAGTATGGATTTCAt                                                                                               >  2:33115/1‑98 (MQ=255)

         cAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTc                                                                                            <  1:2226293/100‑1 (MQ=255)

           accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAg                                                                                          >  2:3405662/1‑100 (MQ=255)

           accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAg                                                                                          <  2:4014533/100‑1 (MQ=255)

            ccaccTTCTCAACAAGCTCAAGTAACGCCTGCTCAATGACTTTGTAAGTAAAAGGAGTAGGTAATTCAAGCAGTGCGTTAAGTTTGGTTTTCATGTc                                                                                            <  2:4038795/97‑1 (MQ=255)

            ccaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGc                                                                                         <  1:3174058/100‑1 (MQ=255)

            ccaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGc                                                                                         <  2:3335615/100‑1 (MQ=255)

              accTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTc                                                                                       >  2:1806345/1‑100 (MQ=255)

                 tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCgg                                                                                    >  1:3312870/1‑100 (MQ=255)

                 tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCgg                                                                                    >  1:3436754/1‑100 (MQ=255)

                 tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCg                                                                                     >  2:56395/1‑99 (MQ=255)

                   aTCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGcg                                                                                  <  1:2036061/100‑1 (MQ=255)

                    tCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGcgt                                                                                 >  1:1610795/1‑100 (MQ=255)

                    tCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgt                                                                                 <  1:2935044/100‑1 (MQ=255)

                      aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGcgtaa                                                                               <  1:3039362/100‑1 (MQ=255)

                       aCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGcgtaac                                                                              >  1:1451949/1‑100 (MQ=255)

                        ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGcgtaacg                                                                             <  1:3062924/100‑1 (MQ=255)

                        ccAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacg                                                                             <  1:3804550/100‑1 (MQ=255)

                         cAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacgt                                                                            <  2:947208/100‑1 (MQ=255)

                                ggTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTaaa                                                                     >  1:2686798/1‑100 (MQ=255)

                                ggTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTaaa                                                                     <  1:1167095/100‑1 (MQ=255)

                                 gTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAAt                                                                    >  1:1399714/1‑100 (MQ=255)

                                  tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                                   >  1:1025107/1‑100 (MQ=255)

                                  tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                                   >  2:1160962/1‑100 (MQ=255)

                                     cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCa                                                                >  2:506219/1‑100 (MQ=255)

                                      gCCTGCCCCATAACTTTGTAAGTAAAAGGAGTATGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGcn                                                                <  2:3521396/99‑2 (MQ=255)

                                       ccTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAAc                                                              >  2:2628066/1‑100 (MQ=255)

                                          gCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACt                                                             <  2:320714/98‑1 (MQ=255)

                                          gCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTcc                                                           >  1:4095633/1‑100 (MQ=255)

                                             ccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGc                                                        >  1:874337/1‑100 (MQ=255)

                                              cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGcc                                                       <  2:2732474/100‑1 (MQ=255)

                                                 aaCTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGa                                                    <  2:3255905/100‑1 (MQ=255)

                                                     ttGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGc                                                >  1:2456405/1‑100 (MQ=255)

                                                       gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                              <  1:767459/100‑1 (MQ=255)

                                                       gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                              >  1:4051592/1‑100 (MQ=255)

                                                        taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCggg                                             <  2:1644188/100‑1 (MQ=255)

                                                        taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCggg                                             <  2:3084810/100‑1 (MQ=255)

                                                        taagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGCAACGTAATTAAATAGCAACTCCCGccttctggcggg                                             <  1:2908309/100‑11 (MQ=255)

                                                         aagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGa                                            <  2:839442/100‑1 (MQ=255)

                                                          agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCTGCCAGAAGGCGGGAg                                           <  1:2361054/100‑1 (MQ=255)

                                                           gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                           >  2:2921224/1‑99 (MQ=255)

                                                           gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGt                                          <  1:3673412/100‑1 (MQ=255)

                                                             aaaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttt                                        >  1:4112926/1‑100 (MQ=255)

                                                             aaaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttt                                        <  2:2485896/100‑1 (MQ=255)

                                                             aaaaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttt                                        <  2:2677182/100‑1 (MQ=255)

                                                                aGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATTTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGTCAGAAGGCGGGAGtttttt                                     <  2:1978227/100‑1 (MQ=255)

                                                                aGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                     >  1:1957146/1‑100 (MQ=255)

                                                                aGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttttt                                     >  1:2163364/1‑100 (MQ=255)

                                                                 ggAGTAGGTAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttt                                       <  2:1266252/97‑1 (MQ=255)

                                                                 ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                                    >  2:830207/1‑100 (MQ=255)

                                                                 ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTGTCAGGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttgttt                                     >  2:2926466/1‑99 (MQ=255)

                                                                  gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                                   <  1:26938/100‑1 (MQ=255)

                                                                  gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                                   <  1:1233087/100‑1 (MQ=255)

                                                                  gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                                   <  2:1785795/100‑1 (MQ=255)

                                                                  gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                                   <  1:815705/100‑1 (MQ=255)

                                                                  gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTAc                                   <  1:20677/100‑1 (MQ=255)

                                                                      aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGTCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                               >  2:2284430/1‑100 (MQ=255)

                                                                      aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                               >  2:2233780/1‑100 (MQ=255)

                                                                      aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                               >  2:202447/1‑100 (MQ=255)

                                                                      aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                               >  2:1382964/1‑100 (MQ=255)

                                                                      aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                               >  2:1103073/1‑100 (MQ=255)

                                                                      aGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                               >  2:474852/1‑100 (MQ=255)

                                                                       gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                              >  1:1270711/1‑100 (MQ=255)

                                                                       gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                              >  1:1729717/1‑100 (MQ=255)

                                                                       gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                              >  2:3129528/1‑100 (MQ=255)

                                                                       gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                              >  1:360070/1‑100 (MQ=255)

                                                                       gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                              >  1:3577196/1‑100 (MQ=255)

                                                                        ggAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGc                             >  1:2429059/1‑100 (MQ=255)

                                                                         gAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCt                            >  1:697756/1‑100 (MQ=255)

                                                                             tCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                        <  2:2686782/100‑1 (MQ=255)

                                                                              cAAGCAGTTCGTTGAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                       <  2:2101819/100‑1 (MQ=255)

                                                                              cAAGCAGTTCGTTAAGTTTTGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCATAAGGCGGGAGTTTTTTACTGATGCTTAGt                        <  1:573921/99‑1 (MQ=255)

                                                                              cAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                       <  2:2304611/100‑1 (MQ=255)

                                                                              cAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                       >  1:4001693/1‑100 (MQ=255)

                                                                              cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                       >  1:1142300/1‑100 (MQ=255)

                                                                              cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                       >  1:2981913/1‑100 (MQ=255)

                                                                                aGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtata                     <  1:1103079/100‑1 (MQ=255)

                                                                                   aGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                  <  1:516856/100‑1 (MQ=255)

                                                                                   aGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                  >  1:2043950/1‑100 (MQ=255)

                                                                                   aGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                  <  1:2628082/100‑1 (MQ=255)

                                                                                   aGTTCGTTAAGTTTGGTTTTCATGTCAGCTACGGCGTAACGTAATTAAATAGCAACTCACGCAAGAAGGCGGGAGTTTTTTACTGAAGCTTAGTATCTgg                  >  1:3447492/1‑100 (MQ=255)

                                                                                     ttCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgggg                <  2:2380721/100‑1 (MQ=255)

                                                                                     ttCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                  <  2:3664342/98‑1 (MQ=255)

                                                                                       cGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc              >  1:3876384/1‑100 (MQ=255)

                                                                                       cGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc              >  1:1713893/1‑100 (MQ=255)

                                                                                         ttAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACgg            <  2:3460947/100‑1 (MQ=255)

                                                                                             gTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAt        >  2:3975250/1‑100 (MQ=255)

                                                                                                   ttttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACt  <  1:1442788/100‑1 (MQ=255)

                                                                                                   ttttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACt  >  2:2324078/1‑100 (MQ=255)

GTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACT  >  NC_000913/662537‑662735

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