Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 277,817 G→A S179S (TCG→TCA)  SAUSA300_RS01220 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793277,8170GA100.0% 261.9 / NA 91S179S (TCG→TCA) SAUSA300_RS01220membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (80/11);  total (80/11)

TGGATTATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCGGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGTATTTTCATTACTA  >  NC_007793/277730‑277913
                                                                                       |                                                                                                
tGGATTATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGttcttc                                                                                      >  1:4646931/1‑100 (MQ=255)
 ggATTATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTtcttct                                                                                     >  1:314650/1‑100 (MQ=255)
 ggATTATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTtcttct                                                                                     >  1:4973203/1‑100 (MQ=255)
     tATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTa                                                                                 >  1:3325148/1‑100 (MQ=255)
     tATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTa                                                                                 >  1:844395/1‑100 (MQ=255)
     tATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTa                                                                                 >  1:3589311/1‑100 (MQ=255)
     tATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTa                                                                                 >  1:5497633/1‑100 (MQ=255)
      ttGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTAt                                                                                >  1:4725435/2‑100 (MQ=255)
      aTGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTAt                                                                                >  1:3937383/1‑100 (MQ=255)
      aTGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTAt                                                                                >  1:4319056/1‑100 (MQ=255)
       tGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATg                                                                               >  1:4915280/1‑100 (MQ=255)
        gggTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGa                                                                              >  1:1201352/1‑100 (MQ=255)
        gggTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGa                                                                              >  1:2478192/1‑100 (MQ=255)
           gcATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATc                                                                           >  1:4215964/2‑100 (MQ=255)
           atctGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATc                                                                           >  1:1703510/4‑100 (MQ=255)
           acATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATc                                                                           >  1:4543037/2‑100 (MQ=255)
           acATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATc                                                                           >  1:3949209/2‑100 (MQ=255)
            cctGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCa                                                                          >  1:739882/3‑100 (MQ=255)
            cATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCa                                                                          >  1:148833/1‑100 (MQ=255)
            cATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCa                                                                          >  1:4587333/1‑100 (MQ=255)
            cATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCa                                                                          >  1:5489045/1‑100 (MQ=255)
            cATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCa                                                                          >  1:4486978/1‑100 (MQ=255)
            cATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATAAATCa                                                                          >  1:56110/1‑100 (MQ=255)
            agtGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCa                                                                          >  1:5053775/3‑100 (MQ=255)
             ctGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCaa                                                                         >  1:1004042/2‑100 (MQ=255)
             aTGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCaa                                                                         >  1:2201880/1‑100 (MQ=255)
              tGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCaaa                                                                        >  1:343951/1‑100 (MQ=255)
              tGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCaaa                                                                        >  1:3542530/1‑100 (MQ=255)
              tGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCaaa                                                                        >  1:4722652/1‑100 (MQ=255)
              tGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCaaa                                                                        >  1:4743648/1‑100 (MQ=255)
               tgcagcaATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAAt                                                                       >  1:2509512/2‑100 (MQ=255)
               ggCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAAt                                                                       >  1:5144046/1‑100 (MQ=255)
               ggCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAAt                                                                       >  1:5321073/1‑100 (MQ=255)
                 cagcaATTTTTAATAACACTAACGCGATTCCGGTTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGACTCAAATGc                                                                     >  1:1261549/1‑100 (MQ=255)
                  ggcaATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCg                                                                    >  1:918523/2‑100 (MQ=255)
                  agcaATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCg                                                                    >  1:2693087/1‑100 (MQ=255)
                   gcaATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGt                                                                   >  1:3487775/1‑100 (MQ=255)
                             tcaaCACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAatat                                                         >  1:4374812/3‑100 (MQ=255)
                             ataaCACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAatat                                                         >  1:558757/1‑100 (MQ=255)
                              taaCACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATAtt                                                        <  1:3756693/100‑1 (MQ=255)
                               aaCACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTa                                                       >  1:1752263/1‑100 (MQ=255)
                                  ccTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCa                                                    >  1:2895636/2‑100 (MQ=255)
                                  acTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCa                                                    >  1:4610580/1‑100 (MQ=255)
                                   ccgaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAt                                                   >  1:86194/4‑100 (MQ=255)
                                   cTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAt                                                   >  1:4371902/1‑100 (MQ=255)
                                   cTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAt                                                   >  1:2285832/1‑100 (MQ=255)
                                   acaaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAt                                                   >  1:1823279/3‑100 (MQ=255)
                                    cgaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAtt                                                  >  1:2221903/3‑100 (MQ=255)
                                     aaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAttt                                                 >  1:4536541/1‑100 (MQ=255)
                                     aaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAttt                                                 >  1:674108/1‑100 (MQ=255)
                                     aaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAttt                                                 >  1:2066467/1‑100 (MQ=255)
                                     aaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAttt                                                 >  1:3829309/1‑100 (MQ=255)
                                     aaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCAttt                                                 >  1:2838041/1‑100 (MQ=255)
                                     aaCGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGAGTTTAAAAATATTAGCAttt                                                 >  1:3344947/1‑100 (MQ=255)
                                               cGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGttat                                       >  1:1823735/1‑100 (MQ=255)
                                               cGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGttat                                       >  1:4743220/1‑100 (MQ=255)
                                               cGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGttat                                       >  1:324520/1‑100 (MQ=255)
                                                    ctcATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCa                                  >  1:2424404/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:261759/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:368563/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:576728/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:4923611/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:3090958/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:1503937/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:894660/1‑100 (MQ=255)
                                                      cATTTTTATAGGGGCTTTAAAAGTACCATTATCAGGTTTGTTCTTTTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATg                                >  1:3172692/1‑100 (MQ=255)
                                                             gtaGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTa                         >  1:4900418/2‑100 (MQ=255)
                                                              taGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTAt                        >  1:2348361/1‑100 (MQ=255)
                                                              taGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTAt                        >  1:3295087/1‑100 (MQ=255)
                                                              taGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTAt                        >  1:1985586/1‑100 (MQ=255)
                                                              taGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCCTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTAt                        >  1:3242589/1‑100 (MQ=255)
                                                               aGGGGCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTAtc                       <  1:646443/100‑2 (MQ=255)
                                                                 tggCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTgg                     >  1:2278951/2‑100 (MQ=255)
                                                                 gggCTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTgg                     >  1:3025609/1‑100 (MQ=255)
                                                                    cTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTggcgg                  <  1:2457298/100‑1 (MQ=255)
                                                                    cTTTAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTggcga                  <  1:5292851/100‑2 (MQ=255)
                                                                     tttAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTggcggg                 <  1:201453/100‑2 (MQ=255)
                                                                     tttAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTggcggc                 <  1:5510547/100‑1 (MQ=255)
                                                                      ttAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTggcggca                <  1:4593861/100‑2 (MQ=255)
                                                                      ttAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGgcggcg                <  1:2948923/100‑1 (MQ=255)
                                                                       tAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGt               >  1:3141127/1‑100 (MQ=255)
                                                                       tAACAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGt               <  1:2856616/100‑1 (MQ=255)
                                                                         catgccccATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTggcggccaga             >  1:613763/7‑96 (MQ=255)
                                                                         aCAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTggcgtagatc             >  1:2653060/1‑94 (MQ=255)
                                                                         aCAGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGTag             <  1:1662299/100‑2 (MQ=255)
                                                                           tgTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGgcggcgagatc           >  1:1158604/2‑95 (MQ=255)
                                                                           aGTACCATTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGgcggcgagatc           >  1:2560991/1‑95 (MQ=255)
                                                                                 aTTATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGTATTTTCATTa     >  1:1559616/1‑100 (MQ=255)
                                                                                   tATCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGTATTTTCATTACt   <  1:125155/100‑1 (MQ=255)
                                                                                    aTCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGTATTTTCATTACTa  >  1:2692925/1‑100 (MQ=255)
                                                                                    aTCAGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGTATTTTCATTACTa  >  1:4937245/1‑100 (MQ=255)
                                                                                       |                                                                                                
TGGATTATGGGTCATGGCAGCAATTTTTAATAACACTAACGCGATTCCGGGTCTCATTTTTATAGGGGCTTTAACAGTACCATTATCGGGTTTGTTCTTCTTTTATGAATCAAATGCGTTTAAAAATATTAGCATTTTTGAAGTTATTATCATGTTCTTTATTGGCGGCGTATTTTCATTACTA  >  NC_007793/277730‑277913

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