BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	532,380	A→G	100%	intergenic (‑159/‑71)	allS ← / → allA	allD operon transcriptional activator/ureidoglycolate lyase, releasing urea

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	532,380	0	A	G	98.0%	356.2 / ‑2.7	112	intergenic (‑159/‑71)	allS/allA	allD operon transcriptional activator/ureidoglycolate lyase, releasing urea
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (51/55); total (55/57)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 112 total aligned reads displayed.

TATTTTTACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTAAG  >  NC_000913/532283‑532476

tATTTTTACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGtt                                                                                                >  2:681289/1‑100 (MQ=255)

   ttttACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGcc                                                                                             <  2:2436903/100‑1 (MQ=255)

   ttttACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGcc                                                                                             <  2:2405248/100‑1 (MQ=255)

     ttACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTg                                                                                           <  2:2595501/100‑1 (MQ=255)

     ttACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTg                                                                                           <  2:46374/100‑1 (MQ=255)

      tACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTg                                                                                           <  2:1166930/99‑1 (MQ=255)

      tACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGt                                                                                          <  1:314402/100‑1 (MQ=255)

      tACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGt                                                                                          >  1:557017/1‑100 (MQ=255)

          tGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTAtt                                                                                      <  2:2005078/100‑1 (MQ=255)

          tGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTAtt                                                                                      >  1:2831891/1‑100 (MQ=255)

          tGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTAtt                                                                                      <  2:1395812/100‑1 (MQ=255)

          tGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAAGTTTTAGGTTGCCTGTTAtt                                                                                      >  1:754686/1‑100 (MQ=255)

           gATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTg                                                                                     <  2:2275610/100‑1 (MQ=255)

           gATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTg                                                                                     >  1:727671/1‑100 (MQ=255)

           gATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTg                                                                                     <  2:788127/100‑1 (MQ=255)

             tAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGaa                                                                                   <  1:803938/100‑1 (MQ=255)

             tAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGaa                                                                                   >  1:2522899/1‑100 (MQ=255)

             tAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGaa                                                                                   >  2:2465712/1‑100 (MQ=255)

             tAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTAATGTTGTTCTCTTAAGGATTTAGATTGCCTTTTATTGaa                                                                                   >  1:128541/1‑100 (MQ=255)

              aaCCCTGTTCAATTTTTGGGCTAAATCTAGTTTTGGAGAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTACCGGTTATTGaa                                                                                   <  2:2399015/99‑1 (MQ=255)

              aaCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGaaa                                                                                  <  2:1624904/100‑1 (MQ=255)

              aaCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGaaa                                                                                  <  1:1777306/100‑1 (MQ=255)

               aCCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAAc                                                                                 <  2:521191/100‑1 (MQ=255)

                cccTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAAcc                                                                                >  2:308399/1‑100 (MQ=255)

                cccTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAAc                                                                                 <  2:1812900/99‑1 (MQ=255)

                   tGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAg                                                                             <  1:189894/100‑1 (MQ=255)

                    gTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGc                                                                            <  2:1108521/100‑1 (MQ=255)

                     ttCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCt                                                                           >  1:2641814/1‑100 (MQ=255)

                     ttCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCt                                                                           <  1:2988758/100‑1 (MQ=255)

                     ttCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTATTCTCTTAAGGATTTAGGTTGCCTGTTATTGAAACCAAGCt                                                                           >  1:1103776/1‑100 (MQ=255)

                      tCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGAATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTg                                                                          >  2:1725536/1‑100 (MQ=255)

                         aTTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGAcc                                                                       <  2:1816064/100‑1 (MQ=255)

                               ggACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTcgg                                                                 >  2:1950548/1‑100 (MQ=255)

                                gACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTcggc                                                                >  1:2273947/1‑100 (MQ=255)

                                 aCTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTcggcg                                                               >  2:2547777/1‑100 (MQ=255)

                                 aCTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTcggcg                                                               >  1:1981617/1‑100 (MQ=255)

                                    aaaTCTATTTTTGTAAAAATATTCCAACTTTTATATTAATGTTGTTCACTTAACGTTTTAGTTTGCCTGATATTGAAACCAAGATGAACGGTc                                                                   >  2:2469566/1‑93 (MQ=255)

                                    aaaTCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCggtg                                                            <  2:2352949/100‑1 (MQ=255)

                                     aaTCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCggtgg                                                           >  2:2063662/1‑100 (MQ=255)

                                      aTCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCggtggt                                                          >  2:1221179/1‑100 (MQ=255)

                                       tCTAGTTTTGGAAAAATATTCCAATTTGTGTATTTATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGtt                                                         <  2:63651/100‑1 (MQ=255)

                                        cTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTg                                                        <  1:1569659/100‑1 (MQ=255)

                                         tAGTTTTGGAAAAATATTCCAAATTTCGTATTGAAGTTGTGCTATTAAGGTTTTAGGTTGCATGGTATTGAAACCAAGCTGGCCGGTCGGCggtg                                                            >  2:1001139/1‑95 (MQ=255)

                                           gTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCCAGCTGACCGGTCGGCGGTGGTTGAAc                                                     >  2:902546/1‑100 (MQ=255)

                                           gTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAAc                                                     >  1:3039584/1‑100 (MQ=255)

                                           gTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAAc                                                     >  2:2790341/1‑100 (MQ=255)

                                            ttttGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACg                                                    <  1:394964/100‑1 (MQ=255)

                                            ttttGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACg                                                    <  1:2441526/100‑1 (MQ=255)

                                             tttGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACgg                                                   >  1:2551623/1‑100 (MQ=255)

                                             tttGGAAAAATATGCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACgg                                                   <  2:2336509/100‑1 (MQ=255)

                                               tGGAAAACTATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTTAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGaa                                                 >  2:699412/1‑100 (MQ=255)

                                               tGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTTACCGGTCGGCGGTGGTTGAACGGaa                                                 <  2:363989/100‑1 (MQ=255)

                                                ggAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGa                                                  >  2:2656996/1‑98 (MQ=255)

                                                ggAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAAt                                                >  2:2882259/1‑100 (MQ=255)

                                                ggAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAAt                                                <  1:1049875/100‑1 (MQ=255)

                                                ggAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAAt                                                >  1:1867548/1‑100 (MQ=255)

                                                 gAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAAt                                                >  1:441977/1‑99 (MQ=255)

                                                 gAAAAATATTCCAACTTTTGTATTGATGTCGTCCTCTTAAGGATCTAGCCAGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGaa                                                 >  2:1003326/1‑98 (MQ=255)

                                                  aaaaaTATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTa                                              <  1:772959/100‑1 (MQ=255)

                                                   atatttttCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTAt                                             >  2:931246/7‑100 (MQ=255)

                                                    aaaTATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATg                                            <  2:2258684/100‑1 (MQ=255)

                                                      atatTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGtt                                          <  1:1106452/100‑1 (MQ=255)

                                                      atatTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGtt                                          <  2:554547/100‑1 (MQ=255)

                                                       tatTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTa                                         <  1:460686/100‑1 (MQ=255)

                                                       tatTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTa                                         <  1:1972044/100‑1 (MQ=255)

                                                        atTCCAACTTTTGTATTGAGGTTGTTCTCTTAAGGTTTTAGCTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTTGTTTAACGGAATTATGTTAc                                        <  2:1458385/100‑1 (MQ=255)

                                                         ttCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACa                                       >  2:534556/1‑100 (MQ=255)

                                                          gatctcTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACaa                                      <  1:1186267/95‑1 (MQ=255)

                                                           ccAACTTTTGTATTGATGTTGTTTTCTTATGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAg                                     <  1:182520/100‑1 (MQ=255)

                                                           ccAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAg                                     <  1:1003332/100‑1 (MQ=255)

                                                           ccAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAg                                     <  2:2831859/100‑1 (MQ=255)

                                                           ccAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAg                                     <  1:2784887/100‑1 (MQ=255)

                                                           ccAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAg                                     >  1:371690/1‑100 (MQ=255)

                                                               cTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACaaagatc                                 >  2:1186256/1‑95 (MQ=255)

                                                               cTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACa                                 >  1:1031745/1‑100 (MQ=255)

                                                                  ttGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACaaaa                              <  2:665721/100‑1 (MQ=255)

                                                                   tGTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACaaaaa                             <  1:480386/100‑1 (MQ=255)

                                                                    gTATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCAACGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAg                            >  1:419962/1‑100 (MQ=255)

                                                                     tATTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGa                           <  1:1265091/100‑1 (MQ=255)

                                                                      aTTGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGAt                          >  1:2103604/1‑100 (MQ=255)

                                                                       ttGATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATg                         >  2:816158/1‑100 (MQ=255)

                                                                         gATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGaa                       <  1:219309/100‑1 (MQ=255)

                                                                         gATGTTGTTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGaa                       <  1:5153/100‑1 (MQ=255)

                                                                           tgttgtTCTTTTAAGGTTTTAGGTTGCTTGTTTTTGAAACCATGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAAc                     <  1:1712938/100‑1 (MQ=255)

                                                                           tgttgtTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAAc                     >  2:1080002/1‑100 (MQ=255)

                                                                           tgttgtTCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAAc                     <  2:2869219/100‑1 (MQ=255)

                                                                            gttgttCTCTTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACt                    >  1:1641879/1‑100 (MQ=255)

                                                                                ttCTCTTAAGGTTTTAGGTTGCCCGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAg                >  1:2779771/1‑100 (MQ=255)

                                                                                  ctctTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGt              >  2:396251/1‑100 (MQ=255)

                                                                                  ctctTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGt              <  1:2465741/100‑1 (MQ=255)

                                                                                   tctTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTa             <  1:2790373/100‑1 (MQ=255)

                                                                                   tctTAAGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTa             <  2:1981594/100‑1 (MQ=255)

                                                                                       aaGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTAc         >  1:182762/1‑100 (MQ=255)

                                                                                       aaGGTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTAc         >  2:2334275/1‑100 (MQ=255)

                                                                                         ggTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCg       >  1:2193722/1‑100 (MQ=255)

                                                                                          gTTTTAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGt      >  1:2249868/1‑100 (MQ=255)

                                                                                           ttttAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGtt     >  1:1070563/1‑100 (MQ=255)

                                                                                             ttAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTaa   <  1:1950571/100‑1 (MQ=255)

                                                                                              tAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTAAg  <  2:927616/100‑1 (MQ=255)

                                                                                              tAGGTTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTAAg  <  1:2977024/100‑1 (MQ=255)

TATTTTTACATGATAACCCTGTTCAATTTGTGGACTAAATCTAGTTTTGGAAAAATATTCCAACTTTTGTATTGATGTTGTTCTCTTAAGGTTTTAGATTGCCTGTTATTGAAACCAAGCTGACCGGTCGGCGGTGGTTGAACGGAATTATGTTACAAGGACAAAAAGATGAAACTTCAGGTATTACCGTTAAG  >  NC_000913/532283‑532476

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