Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 612,932 2 bp→TG coding (1786‑1787/2844 nt) SAUSA300_RS02915 → serine‑aspartate repeat‑containing protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793612,9320AT98.1% 812.5 / 4.2 572M596L (ATG→TTG) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base A (10/1);  new base T (543/17);  total (553/19)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-03
*NC_007793612,9330TG97.9% 815.8 / 11.3 581M596R (ATG→AGG) SAUSA300_RS02915serine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base T (11/1);  new base G (551/18);  total (562/19)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-05

Only 100 of 582 total aligned reads displayed.

GCAAATGTAGGTACAGATGATGCTGTAGATTCTGATGGACTAACTACAACAGGTGTCATTAAAGACGCTGACAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGA  >  NC_007793/612857‑613029
                                                                           ||                                                                                                   
        aGGTACAGATGATGCTGTAGATTCTGATGGACTAACTACAACAGGTGTCATTAAAGACGCTGACAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑Ca                                                                   >  1:3565558/1‑100 (MQ=255)
                                    ggACTAACTACAACAGGTGTCATTAAAGACGCTGACAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAGTTTAGGTGATTATGTTTGGt                                       >  1:2885044/1‑100 (MQ=255)
                                                    gtacaaTTAAAGATGCTGACAACTGTACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCtaataa                          <  1:6617201/95‑1 (MQ=9)
                                                    gtacaaTTAAAGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCtaataa                          <  1:1167485/95‑1 (MQ=9)
                                                    gtacaaTTAAAGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCtaataa                          <  1:3964811/95‑1 (MQ=9)
                                                     cacaaTTAAAGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATaaa                         >  1:4847653/5‑100 (MQ=9)
                                                            taaGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGtt                  >  1:2275023/2‑99 (MQ=9)
                                                            aaaGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATgggt                  >  1:714512/1‑98 (MQ=9)
                                                            aaaGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGtt                  >  1:424201/1‑99 (MQ=9)
                                                             caGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGtta                 >  1:3401317/2‑98 (MQ=9)
                                                             aaGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGtta                 >  1:3653776/1‑98 (MQ=9)
                                                             aaGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGtta                 >  1:188923/1‑98 (MQ=9)
                                                             aaGACGCTGACAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACa              <  1:3255310/100‑1 (MQ=255)
                                                              ccatgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGttac                >  1:2452761/5‑97 (MQ=9)
                                                              aGATGCTGACAACTGGACTTTGGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGttac                >  1:687837/1‑97 (MQ=9)
                                                               ttgccTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGttacc               >  1:5134289/5‑96 (MQ=9)
                                                               tgtgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGAtcggaaga               >  1:1565014/4‑93 (MQ=2)
                                                               aatgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGttacc               >  1:4753185/4‑96 (MQ=9)
                                                                atgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGttacca              >  1:1520299/3‑95 (MQ=9)
                                                                atgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGttacca              >  1:6370946/3‑95 (MQ=9)
                                                                 tgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCaa             >  1:5882260/2‑100 (MQ=9)
                                                                 tgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCaa             >  1:5019137/2‑100 (MQ=9)
                                                                 tgCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCaa             >  1:5065058/2‑100 (MQ=9)
                                                                 atcTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCaa             >  1:2315523/3‑100 (MQ=9)
                                                                 accTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCaa             >  1:5387335/3‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:3182555/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:1814733/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:1224250/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:835820/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:2728135/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:5421531/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:2360185/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:4413722/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:2179713/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:2049185/1‑100 (MQ=9)
                                                                  gCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:1754200/1‑100 (MQ=9)
                                                                  ccTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:5104545/2‑100 (MQ=9)
                                                                  ccTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:1066215/2‑100 (MQ=9)
                                                                  ccTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:3064470/2‑100 (MQ=9)
                                                                  ccTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:6596905/2‑100 (MQ=9)
                                                                  ccTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAg            >  1:1427605/2‑100 (MQ=9)
                                                                   gtGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGa           >  1:1641541/2‑100 (MQ=9)
                                                                   cTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGa           >  1:3454881/1‑100 (MQ=9)
                                                                   cTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGa           >  1:6414719/1‑100 (MQ=9)
                                                                   cTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGa           >  1:5836615/1‑100 (MQ=9)
                                                                   cTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGa           >  1:4186463/1‑100 (MQ=9)
                                                                   cTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTAGGTATGGTACGATTCTAATAAAGATGGTTACCAAGa           >  1:3777145/1‑100 (MQ=9)
                                                                   actacaactGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGa           >  1:1370092/4‑100 (MQ=9)
                                                                    ttacaactGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGac          >  1:4651520/3‑99 (MQ=9)
                                                                    tGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGac          >  1:5553416/1‑99 (MQ=9)
                                                                    tGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGac          >  1:784101/1‑99 (MQ=9)
                                                                    tGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGac          >  1:4044319/1‑99 (MQ=9)
                                                                    tGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGac          >  1:5181085/1‑99 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACCGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2667451/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTCGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:3503066/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGGTGGTTACCAAGaca         >  1:4961695/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGCTGGTTACCAAGaca         >  1:4917752/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:6489332/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:5150918/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:5234056/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:6273934/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:6204373/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:560608/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:6049170/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:5932404/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:570513/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:1006628/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:938470/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:5789777/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2953948/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:1098938/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:1346451/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:1678061/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:1716047/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:1847066/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:1947129/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2098260/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2203334/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2417003/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:255648/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:259821/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2754402/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2816754/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:2855950/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:4696660/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:3213791/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:3295904/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:3619173/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:3723584/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:3830117/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:410265/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:4152603/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:4361516/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:4482381/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:4548943/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:4602354/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATAGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:4278391/1‑98 (MQ=9)
                                                                     gACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGGGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:300041/1‑98 (MQ=9)
                                                                     aacaactGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGaca         >  1:6085326/2‑98 (MQ=9)
                                                                      acaactGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGacaa        >  1:3895930/1‑97 (MQ=9)
                                                                           ||                                                                                                   
GCAAATGTAGGTACAGATGATGCTGTAGATTCTGATGGACTAACTACAACAGGTGTCATTAAAGACGCTGACAACATGACATTAGATAGTGGATTCTACAAAACAC‑‑‑CAAAATATAGTTTAGGTGATTATGTTTGGTACGACAGTAATAAAGATGGTAAACAAGATTCGACTGA  >  NC_007793/612857‑613029

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