Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007790 744 T→C D155D (GAT→GAC)  SAUSA300_RS14705 → replication protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077907440TC99.1% 321.8 / ‑3.5 110coding (465/717 nt)SAUSA300_RS14705replication protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (91/18);  minor base . (1/0);  total (92/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 110 total aligned reads displayed.

TCGGCAATCAAAGAGACCTCAAAATATTCGGTTAAATCATCTGATTTTTTAACTGATGATGATG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTG  >  NC_007790/700‑837
                                            |                                                                                              
gggccactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAg                                        >  1:5068987/11‑100 (MQ=255)
ctaacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAg                                        >  1:5503773/11‑100 (MQ=255)
cacccactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAg                                        >  1:520800/11‑100 (MQ=255)
ataacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAg                                        >  1:5889292/11‑100 (MQ=255)
ataacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAg                                        >  1:2923452/11‑100 (MQ=255)
ataacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAg                                        >  1:2630431/11‑100 (MQ=255)
 taacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAgg                                       >  1:5421800/10‑100 (MQ=255)
 atacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAgg                                       >  1:4765385/10‑100 (MQ=255)
  cacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:5829283/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCGTCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:4233513/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACGAAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAgg                                       >  1:3129256/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:2166547/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:5120625/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:2581656/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:4709477/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:4006692/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:3585548/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:5812252/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:1519846/9‑100 (MQ=255)
  aacactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGt                                      >  1:5851325/9‑100 (MQ=255)
   tcactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAggat                                     >  1:3609182/8‑98 (MQ=255)
   ccactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGtt                                     >  1:4888556/8‑100 (MQ=255)
   acactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGtt                                     >  1:823842/8‑100 (MQ=255)
   acactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGtt                                     <  1:2304890/93‑1 (MQ=255)
   acactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGtt                                     >  1:1226597/8‑100 (MQ=255)
   acactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTCAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGtt                                     >  1:5588745/8‑100 (MQ=255)
    ccctcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGttt                                    >  1:517141/7‑100 (MQ=255)
    cactcgaaGAGACATCAGAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGttt                                    <  1:1525180/94‑1 (MQ=255)
    cactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGttt                                    >  1:2783872/7‑100 (MQ=255)
    cactcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACA‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGttt                                    >  1:4067264/7‑100 (MQ=255)
     cctcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACCGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTg                                   >  1:5385068/6‑100 (MQ=255)
      ctcgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGt                                  >  1:5412572/5‑100 (MQ=255)
        cgaaGAGACATCAAAATATTCCGTTAAATCCTCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTAt                                >  1:933993/3‑100 (MQ=255)
        cgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTAt                                >  1:3707978/3‑100 (MQ=255)
        cgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTAt                                >  1:4962983/3‑100 (MQ=255)
        cgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTAt                                >  1:4493320/3‑100 (MQ=255)
        cgaaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTAt                                >  1:1823714/3‑100 (MQ=255)
         ccagagACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATc                               >  1:1633373/3‑100 (MQ=255)
           ctagaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGa                             >  1:3788997/3‑100 (MQ=255)
           cctgaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGa                             >  1:3013273/4‑100 (MQ=255)
           agagACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGa                             >  1:114282/1‑100 (MQ=255)
           acagaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGa                             >  1:5412526/3‑100 (MQ=255)
            gagaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaa                            >  1:5352579/1‑100 (MQ=255)
            ctgaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaa                            >  1:1781697/3‑100 (MQ=255)
            ctgaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaa                            >  1:5934962/3‑100 (MQ=255)
             tgaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:5399461/2‑100 (MQ=255)
             ggaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACGAAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaa                            >  1:4404528/2‑100 (MQ=255)
             ggaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:4594466/2‑100 (MQ=255)
             cgaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:569423/2‑100 (MQ=255)
             agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCaaaa                           >  1:221295/1‑100 (MQ=255)
             agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:849209/1‑100 (MQ=255)
             agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:2942558/1‑100 (MQ=255)
             agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:5301329/1‑100 (MQ=255)
             agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:5609060/1‑100 (MQ=255)
             agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:3894046/1‑100 (MQ=255)
              tgacatcaaaTATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGCAATTATAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:4056534/8‑100 (MQ=255)
              tacatCAAAATATTCCGTTAAATCATCTGA‑CTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAAc                         >  1:5033417/5‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACGAAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaa                           >  1:1985250/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTTAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:662995/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:4805610/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:4842300/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:2410384/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:758783/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:3735876/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:5077404/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGaaaa                          >  1:5256506/1‑100 (MQ=255)
              gaCATCAAAATATTCCGTTAAATCATCTGACCTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGATTGTATCGaaaa                          >  1:3211523/1‑100 (MQ=255)
               tcatCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAAc                         >  1:2756275/4‑100 (MQ=255)
                catCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACg                        >  1:1804114/3‑100 (MQ=255)
                 atctaaaTATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGt                       >  1:4518206/5‑100 (MQ=255)
                   cccggcatTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTAt                     >  1:5485925/7‑100 (MQ=255)
                   cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACGAAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTa                      >  1:2776859/1‑100 (MQ=255)
                   cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTAt                     >  1:1185418/1‑100 (MQ=255)
                   cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTAt                     >  1:4515747/1‑100 (MQ=255)
                   cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTAt                     >  1:5592680/1‑100 (MQ=255)
                   cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTAt                     >  1:1053008/1‑100 (MQ=255)
                      aaTATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGtt                  <  1:5629997/100‑1 (MQ=38)
                       atatTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTa                 <  1:2243138/100‑1 (MQ=38)
                        tgtcccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTaa                >  1:503544/7‑100 (MQ=37)
                        caggccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTaa                >  1:1430897/7‑100 (MQ=37)
                         gttCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAg               >  1:5079888/2‑100 (MQ=37)
                         atTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAg               >  1:1847466/1‑100 (MQ=38)
                          ctccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGt              >  1:5052367/5‑100 (MQ=37)
                            ccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTa            >  1:829227/3‑100 (MQ=37)
                            ccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTa            >  1:2937957/3‑100 (MQ=37)
                             ccgatctacatcCGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAt           <  1:1689558/92‑1 (MQ=25)
                                tAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGGAATTGTAAGTGATTTAGAAGAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt        <  1:554601/100‑1 (MQ=37)
                                tAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt        <  1:5922058/100‑1 (MQ=255)
                                tAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt        >  1:342076/1‑100 (MQ=38)
                                tAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt        <  1:2554331/100‑1 (MQ=255)
                                tAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt        <  1:3669135/100‑1 (MQ=255)
                                tAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt        <  1:3307914/100‑1 (MQ=255)
                                 ttccgatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggtg       <  1:618255/95‑1 (MQ=38)
                                  aaTCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggtgg      <  1:2884890/100‑1 (MQ=255)
                                  aaTCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggtgc      <  1:1968946/100‑2 (MQ=255)
                                  aaTCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggtgc      <  1:5175530/100‑2 (MQ=255)
                                   atcatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGa     <  1:1534908/100‑1 (MQ=255)
                                    tcatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGAt    <  1:1437035/100‑1 (MQ=255)
                                    tcatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGAt    <  1:2300193/100‑1 (MQ=255)
                                      atcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGAttt  <  1:731816/100‑2 (MQ=255)
                                            |                                                                                              
TCGGCAATCAAAGAGACCTCAAAATATTCGGTTAAATCATCTGATTTTTTAACTGATGATGATG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTG  >  NC_007790/700‑837

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑