Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,013,607 T→A V121E (GTG→GAG)  SAUSA300_RS04960 → Serine protease HtrA‑like

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,013,6070TA98.9% 243.6 / NA 90V121E (GTG→GAG) SAUSA300_RS04960Serine protease HtrA‑like
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (89/0);  minor base G (1/0);  total (90/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGTGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAG  >  NC_007793/1013510‑1013702
                                                                                                 |                                                                                               
aaTCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAgaga                                                                                               >  1:1797174/1‑100 (MQ=255)
aaTCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAgaga                                                                                               >  1:611255/1‑100 (MQ=255)
 aTCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGaa                                                                                              >  1:4132692/1‑100 (MQ=255)
        tAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCaaaa                                                                                       >  1:850576/1‑100 (MQ=255)
         aGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAat                                                                                      >  1:18665/1‑100 (MQ=255)
         aGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAat                                                                                      >  1:2360301/1‑100 (MQ=255)
                cAGCATATAACAGAGATCATTATTTAACAGACGATGTCTCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGa                                                                               >  1:930501/1‑100 (MQ=255)
                 aGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAg                                                                              >  1:5328857/1‑100 (MQ=255)
                  gCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGc                                                                             >  1:6413750/1‑100 (MQ=255)
                   cctataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCa                                                                            >  1:1562545/3‑100 (MQ=255)
                    ctataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaa                                                                           >  1:5310859/2‑100 (MQ=255)
                      ataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaaga                                                                         >  1:3710444/1‑98 (MQ=255)
                      ataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaaca                                                                         >  1:6472864/1‑98 (MQ=255)
                      ataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaaaa                                                                         >  1:2122715/1‑100 (MQ=255)
                      ataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaaaa                                                                         >  1:1134949/1‑100 (MQ=255)
                      ataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaaaa                                                                         >  1:5460468/1‑100 (MQ=255)
                                        tAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATc                                                       >  1:2734928/1‑100 (MQ=255)
                                        tAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATc                                                       >  1:1922411/1‑100 (MQ=255)
                                        gaaCAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATc                                                       >  1:3040926/2‑100 (MQ=255)
                                         caCAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTAAAGCGACTTTATCa                                                      >  1:4870135/2‑100 (MQ=255)
                                          aCAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCaa                                                     >  1:4467100/1‑100 (MQ=255)
                                          aCAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCaa                                                     >  1:1644502/1‑100 (MQ=255)
                                           ccgACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAAc                                                    >  1:1938643/3‑100 (MQ=255)
                                              ccGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTaa                                                 >  1:6322536/2‑100 (MQ=255)
                                                         aaaaaaCAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATa                                      >  1:6038488/1‑100 (MQ=255)
                                                         aaaaaaCAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATa                                      >  1:6223000/1‑100 (MQ=255)
                                                         aaaaaaCAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATCATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATa                                      >  1:649928/1‑100 (MQ=255)
                                                                         tAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAGGTAGAATCAACTGa                      >  1:4917912/1‑100 (MQ=255)
                                                                              tCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATga                 >  1:2553572/1‑100 (MQ=255)
                                                                              ccAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATga                 >  1:5572725/2‑100 (MQ=255)
                                                                              ccAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATga                 >  1:1970110/2‑100 (MQ=255)
                                                                               cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:2070139/1‑100 (MQ=255)
                                                                               cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:4302125/1‑100 (MQ=255)
                                                                               cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:1329728/1‑100 (MQ=255)
                                                                               cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:2827469/1‑100 (MQ=255)
                                                                                tgTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:1821512/2‑100 (MQ=255)
                                                                                cctGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:900683/3‑100 (MQ=255)
                                                                                attGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:1062033/3‑100 (MQ=255)
                                                                                aGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:5738136/1‑100 (MQ=255)
                                                                                aGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:2550719/1‑100 (MQ=255)
                                                                                 ttGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaa              >  1:1106433/2‑100 (MQ=255)
                                                                                  tGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaaa             >  1:1040651/1‑100 (MQ=255)
                                                                                  tGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaaa             >  1:5797629/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAGCCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:5705264/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTGGAATCAACTGAAATGAGAAAg            >  1:2187282/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:3399440/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:609832/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:1917233/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:5432318/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:5058502/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:4557366/1‑100 (MQ=255)
                                                                                   ggACCAAGATACAGAGAAATCAAAATATTACGAGCAAAATTCTGAAGCGACTTTACCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:1442862/1‑100 (MQ=255)
                                                                                    ttccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGc           >  1:5430098/3‑100 (MQ=255)
                                                                                    taCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGc           >  1:3405481/2‑100 (MQ=255)
                                                                                    gACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGc           >  1:2992831/1‑100 (MQ=255)
                                                                                    gACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGc           >  1:3857397/1‑100 (MQ=255)
                                                                                     cccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:2701620/2‑100 (MQ=255)
                                                                                     atcAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:1178430/3‑100 (MQ=255)
                                                                                     aCCAAGATGCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:1959236/1‑100 (MQ=255)
                                                                                     aCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:890943/1‑100 (MQ=255)
                                                                                     aCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:1583530/1‑100 (MQ=255)
                                                                                      ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:3733758/1‑100 (MQ=255)
                                                                                      ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:3570973/1‑100 (MQ=255)
                                                                                      ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:4350765/1‑100 (MQ=255)
                                                                                      ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:395177/1‑100 (MQ=255)
                                                                                      ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:5010925/1‑100 (MQ=255)
                                                                                       ccaGATACAGGGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTaa        >  1:6385082/3‑100 (MQ=255)
                                                                                        taGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAg       >  1:5025725/2‑100 (MQ=255)
                                                                                        gaGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAg       >  1:1354524/2‑100 (MQ=255)
                                                                                        caGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAg       >  1:5592011/2‑100 (MQ=255)
                                                                                        caGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAg       >  1:4148446/2‑100 (MQ=255)
                                                                                        aaGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAg       >  1:2121653/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAACATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:6519210/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGTGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGTAAGCTAAGt      >  1:5878292/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaagatcggaa      >  1:3374328/1‑92 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:169389/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:763295/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:744784/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:1468999/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:28063/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:5208054/1‑100 (MQ=255)
                                                                                         aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:5689098/1‑100 (MQ=255)
                                                                                          gATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGtt     >  1:6210938/1‑100 (MQ=255)
                                                                                            tACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCa   >  1:407285/1‑100 (MQ=255)
                                                                                             aCAGAGGAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:5972143/1‑100 (MQ=255)
                                                                                             aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:1785515/1‑100 (MQ=255)
                                                                                             aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:6472784/1‑100 (MQ=255)
                                                                                             aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:5255290/1‑100 (MQ=255)
                                                                                             aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:2138358/1‑100 (MQ=255)
                                                                                             aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:490701/1‑100 (MQ=255)
                                                                                             aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:2444079/1‑100 (MQ=255)
                                                                                             aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGACAAGTTCAg  >  1:1969642/1‑100 (MQ=255)
                                                                                                 |                                                                                               
AATCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGTGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAG  >  NC_007793/1013510‑1013702

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: ‑