Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 304,266 G→A E317E (GAG→GAA)  SAUSA300_RS01345 → glycosyl transferase family 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793304,2660GA99.5% 1807.2 / ‑4.8 626E317E (GAG→GAA) SAUSA300_RS01345glycosyl transferase family 2
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (619/4);  total (622/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 638 total aligned reads displayed.

AGTTCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAGGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAGGCATT  >  NC_007793/304170‑304362
                                                                                                |                                                                                                
aGTTCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGgaagaa                                                                                               >  1:431325/1‑100 (MQ=255)
   tCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGt                                                                                            >  1:2433626/1‑100 (MQ=255)
   tCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGt                                                                                            >  1:535072/1‑100 (MQ=255)
   tCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGt                                                                                            >  1:3184686/1‑100 (MQ=255)
   tCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGt                                                                                            >  1:1507634/1‑100 (MQ=255)
   tCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACCGTggaagatcgg                                                                                            >  1:4058965/1‑96 (MQ=255)
    ccAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGtt                                                                                           >  1:1535449/1‑100 (MQ=255)
    ccAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGtt                                                                                           >  1:2680135/1‑100 (MQ=255)
    ccAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGtt                                                                                           >  1:2588357/1‑100 (MQ=255)
    ccAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGtt                                                                                           >  1:5124280/1‑100 (MQ=255)
    ccAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGtt                                                                                           >  1:1269350/1‑100 (MQ=255)
    ccAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGtt                                                                                           >  1:2330302/1‑100 (MQ=255)
     tccgaCGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACCGTggaagatcggaa                                                                                          >  1:1647914/5‑94 (MQ=255)
     cAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTa                                                                                          >  1:3909751/1‑100 (MQ=255)
      aGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTAc                                                                                         >  1:5005515/1‑100 (MQ=255)
      aGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGAAAGAAAGTTAc                                                                                         >  1:525889/1‑100 (MQ=255)
       gAACGTGTAGATGAATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTggaagatcggaaga                                                                                        >  1:1237941/1‑92 (MQ=255)
                gATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTc                                                                               >  1:964684/1‑100 (MQ=255)
                caTGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTc                                                                               >  1:5623411/2‑100 (MQ=255)
                 aTGCATTGGTGATGAGTAAATTACGTCCATTGTTGTACTACGCGAGTGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCa                                                                              >  1:2740742/1‑100 (MQ=255)
                 aTGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAGGAAAGTTACCGTCAAGGTCa                                                                              >  1:2305071/1‑100 (MQ=255)
                   gCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAAt                                                                            >  1:2560853/1‑100 (MQ=255)
                   gCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAAt                                                                            >  1:106184/1‑100 (MQ=255)
                   gCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAAt                                                                            >  1:2102073/1‑100 (MQ=255)
                      ttGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAAtact                                                                         >  1:43553/1‑100 (MQ=255)
                       tGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAAtacta                                                                        >  1:3156139/1‑100 (MQ=255)
                           gatgaGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGt                                                                    >  1:5178881/1‑100 (MQ=255)
                                       aCGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTa                                                        >  1:1595429/1‑100 (MQ=255)
                                        tgACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTATCGTCAAGGTCAATACTACCGTTTTGATATTGTAg                                                       >  1:304591/2‑100 (MQ=255)
                                        accccATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAg                                                       >  1:334807/4‑100 (MQ=255)
                                           gcATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATg                                                    >  1:2464992/2‑100 (MQ=255)
                                                cttgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGCTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATTGTaaa                                               >  1:4925708/2‑100 (MQ=255)
                                                attgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTaaa                                               >  1:4906010/2‑100 (MQ=255)
                                                 ttgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAt                                              >  1:1972737/1‑100 (MQ=255)
                                                 ttgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTAACGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAt                                              >  1:567165/1‑100 (MQ=255)
                                                 tggCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAt                                              >  1:6567/3‑100 (MQ=255)
                                                 tcgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAt                                              >  1:1900541/3‑100 (MQ=255)
                                                 ctgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAt                                              >  1:291962/2‑100 (MQ=255)
                                                  tgCGCTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:2147434/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGACAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:4158625/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:3543487/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:5057242/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:2052518/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:482510/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:2251777/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:371811/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:5505434/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:5545657/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:5758268/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:3108205/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:1808969/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGACAATACTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:5424633/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGAGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGCAAAtt                                             >  1:1122652/1‑100 (MQ=255)
                                                  tgCACTACGCGAGAGCGAAAGATATAGACAACTATAACACTGTGGCAGAAAGTTACCGTCAAGGTCAATCCTACCGTTTTGATATTGTAGATGGTAAAtt                                             >  1:5583054/1‑100 (MQ=255)
                                                   gtaCTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTa                                            >  1:4520810/3‑100 (MQ=255)
                                                   gCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGTTGGTAAATTa                                            >  1:5401607/1‑100 (MQ=255)
                                                   gCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTTAATTa                                            >  1:785892/1‑100 (MQ=255)
                                                   gCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTa                                            >  1:1016848/1‑100 (MQ=255)
                                                   gCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTa                                            >  1:4485939/1‑100 (MQ=255)
                                                   gCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTa                                            >  1:3827790/1‑100 (MQ=255)
                                                   gCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTa                                            >  1:69122/1‑100 (MQ=255)
                                                   gCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGATACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTa                                            >  1:442009/1‑100 (MQ=255)
                                                       caCGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACa                                        >  1:4752085/2‑100 (MQ=255)
                                                        aCGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACat                                       >  1:141425/1‑100 (MQ=255)
                                                        aCGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACat                                       >  1:5223998/1‑100 (MQ=255)
                                                        aCGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACat                                       >  1:754953/1‑100 (MQ=255)
                                                        aCGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACat                                       >  1:4614802/1‑100 (MQ=255)
                                                        aCGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACat                                       >  1:4100322/1‑100 (MQ=255)
                                                                gCGAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACattcaattca                               >  1:3996789/1‑100 (MQ=255)
                                                                  gAAAGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATtcaattcaat                             >  1:3946421/1‑100 (MQ=255)
                                                                   aaaGATATAGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATg                            <  1:385725/100‑1 (MQ=255)
                                                                         ctaGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCg                      >  1:5798151/2‑100 (MQ=255)
                                                                          atgccAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGa                     >  1:5716690/5‑100 (MQ=255)
                                                                           aGACAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGaa                    >  1:3598047/1‑100 (MQ=255)
                                                                              cAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAgatcggaag                 >  1:1164162/1‑92 (MQ=255)
                                                                              cAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCa                 >  1:2837956/1‑100 (MQ=255)
                                                                              cAACTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCa                 >  1:3376278/1‑100 (MQ=255)
                                                                               gccTATAGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGacagat                >  1:4575264/3‑95 (MQ=255)
                                                                                   ataGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACtt            >  1:3002897/1‑100 (MQ=255)
                                                                                   ataGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACtt            >  1:876683/1‑100 (MQ=255)
                                                                                   ataGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACtt            >  1:2879140/1‑100 (MQ=255)
                                                                                   ataGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACtt            >  1:1694687/1‑100 (MQ=255)
                                                                                   ataGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACtt            >  1:1724002/1‑100 (MQ=255)
                                                                                   ataGAACTGTGGAAGAAAGTTACCGCCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACtt            >  1:932770/1‑100 (MQ=255)
                                                                                     aGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGaagatcggaaga          >  1:1864953/1‑90 (MQ=255)
                                                                                     aGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTa          >  1:2531225/1‑100 (MQ=255)
                                                                                     aGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTa          >  1:3283685/1‑100 (MQ=255)
                                                                                     aGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTa          >  1:3250081/1‑100 (MQ=255)
                                                                                     aGAACTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTAATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTa          >  1:239071/1‑100 (MQ=255)
                                                                                       cccTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTaaa        >  1:3478452/3‑100 (MQ=255)
                                                                                        gcTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAg       >  1:3663052/2‑100 (MQ=255)
                                                                                        aCTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAg       >  1:2636683/1‑100 (MQ=255)
                                                                                        aCTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAg       >  1:1370390/1‑100 (MQ=255)
                                                                                         cTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAgg      >  1:4264853/1‑100 (MQ=255)
                                                                                         cTGTGGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAgg      >  1:3780732/1‑100 (MQ=255)
                                                                                           gtgGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAGGCa    >  1:3400782/1‑100 (MQ=255)
                                                                                           gtgGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAGGCa    >  1:2790113/1‑100 (MQ=255)
                                                                                           gtgGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAGGCa    >  1:4690230/1‑100 (MQ=255)
                                                                                           gtgGAAGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAGGCa    >  1:1772030/1‑100 (MQ=255)
                                                                                           gtgGAAGAAAGTTAACGTCAAGGTCAATACTACCG‑TTTGATTTTGGAGATGGGAAATTAAACATTCAATTCAATGAAGGCGAACCATCCTTTAAAGGCAt   >  1:4186771/1‑100 (MQ=255)
                                                                                                |                                                                                                
AGTTCCAGAACGTGTAGATGCATTGGTGATGAGTAAATTACGACCATTGTTGCACTACGCGAGAGCGAAAGATATAGACAACTATAGAACTGTGGAGGAAAGTTACCGTCAAGGTCAATACTACCGTTTTGATATTGTAGATGGTAAATTAAACATTCAATTCAATGAAGGCGAACCATACTTTAAAGGCATT  >  NC_007793/304170‑304362

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