Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,847,298 A→G K146K (AAA→AAG sppA → protease IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,847,2980AG100.0% 737.8 / NA 229K146K (AAA→AAGsppAprotease IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (119/110);  total (119/110)

Only 100 of 235 total aligned reads displayed.

TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387
                                                                                          |                                                                                         
tgtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGcc                                                                                         <  1:6287453/93‑1 (MQ=255)
 gtgATGGATTTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCg                                                                                        >  2:4864403/1‑93 (MQ=255)
  tgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCgg                                                                                       <  2:4941456/93‑1 (MQ=255)
  tgAGGGATCTGAAAAGCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCgg                                                                                       <  2:1204164/93‑1 (MQ=255)
   gATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGt                                                                                      <  2:5220106/93‑1 (MQ=255)
    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     <  2:5892378/93‑1 (MQ=255)
     tGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGttt                                                                                    >  1:4363890/1‑93 (MQ=255)
      ggATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   >  2:5323601/1‑93 (MQ=255)
       gATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAATACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGGTTAt                                                                                  >  1:878310/1‑93 (MQ=255)
        aTCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATg                                                                                 <  1:5773362/93‑1 (MQ=255)
         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                >  2:1308777/1‑93 (MQ=255)
         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                <  2:1075941/93‑1 (MQ=255)
         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                >  2:5144106/1‑93 (MQ=255)
         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGGTTCGTGACAGCGGGAAGACGGTTTATGc                                                                                >  2:4009443/1‑93 (MQ=255)
           tGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCg                                                                              >  2:5462092/1‑93 (MQ=255)
            gAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGt                                                                             <  2:1181025/93‑1 (MQ=255)
            gAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGt                                                                             <  2:1855003/93‑1 (MQ=255)
            gAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGt                                                                             <  2:5936768/93‑1 (MQ=255)
              aaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTg                                                                           <  1:5533633/93‑1 (MQ=255)
                aaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGc                                                                         >  2:2407110/1‑93 (MQ=255)
                aaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTAACAGCGGGAAGCCGGTTTATGCCGTTGGc                                                                         >  2:3808186/1‑93 (MQ=255)
                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCgt                                                                       <  2:6317726/93‑2 (MQ=255)
                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCga                                                                       >  1:308932/1‑93 (MQ=255)
                   ttCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCgag                                                                      >  1:375597/1‑93 (MQ=255)
                     cGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGaa                                                                    <  2:35219/93‑1 (MQ=255)
                     cGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGaa                                                                    >  1:5077616/1‑93 (MQ=255)
                      gCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   <  2:4795007/93‑1 (MQ=255)
                      gCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   >  2:2437324/1‑93 (MQ=255)
                      gCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   <  1:137828/93‑1 (MQ=255)
                      gCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   >  2:6069354/1‑93 (MQ=255)
                       cAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACt                                                                  <  2:3707253/93‑1 (MQ=255)
                        aGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTa                                                                 >  2:1728353/1‑93 (MQ=255)
                        aGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTa                                                                 >  2:3187973/1‑93 (MQ=255)
                        aGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTa                                                                 >  1:1728159/1‑93 (MQ=255)
                         ggcggcGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                >  2:2548932/1‑93 (MQ=255)
                         ggcggcGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                >  1:2583319/1‑93 (MQ=255)
                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  2:1682984/93‑1 (MQ=255)
                             gcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGcc                                                            >  2:1522014/1‑93 (MQ=255)
                             gcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGcc                                                            >  1:3821749/1‑93 (MQ=255)
                             gcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCAGGGTTATGCCGTTGGCGAGAACTACAGcc                                                            >  1:1792/1‑93 (MQ=255)
                               gACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAg                                                          <  2:3621080/93‑1 (MQ=255)
                               gACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAg                                                          <  1:3540605/93‑1 (MQ=255)
                                aCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                         <  1:1143034/93‑1 (MQ=255)
                                aCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                         >  1:557526/1‑93 (MQ=255)
                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        <  2:6130159/93‑1 (MQ=255)
                                  cAACCGTCTATGCAGTACATCGGCATAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgggg                                                       <  1:2790833/93‑1 (MQ=255)
                                  cAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgggg                                                       >  2:4768330/1‑93 (MQ=255)
                                  cAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgggg                                                       <  2:250854/93‑1 (MQ=255)
                                   aaCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                      >  2:1030718/1‑93 (MQ=255)
                                    aCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCa                                                     >  2:645083/1‑93 (MQ=255)
                                     ccGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCaa                                                    >  2:6236128/1‑93 (MQ=255)
                                      cGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAat                                                   >  2:4587607/1‑93 (MQ=255)
                                       gTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAata                                                  <  2:5610416/93‑1 (MQ=255)
                                       gTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAata                                                  >  2:2194806/1‑93 (MQ=255)
                                       gTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAata                                                  <  1:6552987/93‑1 (MQ=255)
                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  1:6077533/1‑93 (MQ=255)
                                          tATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAAtatta                                               >  2:597759/1‑93 (MQ=255)
                                                gTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         <  2:5680112/93‑1 (MQ=255)
                                                gTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         >  1:1864548/1‑93 (MQ=255)
                                                 tACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGcc                                        <  1:4137298/93‑1 (MQ=255)
                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGtt                                    <  1:3475378/93‑1 (MQ=255)
                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGtt                                    >  2:5156273/1‑93 (MQ=255)
                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGtt                                    <  1:2206610/93‑1 (MQ=255)
                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGtt                                    <  1:2267590/93‑1 (MQ=255)
                                                      cGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGttt                                   >  2:3426676/1‑93 (MQ=255)
                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTc                                  >  1:6115906/1‑93 (MQ=255)
                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTc                                  >  2:3225657/1‑93 (MQ=255)
                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTc                                  <  2:3759178/93‑1 (MQ=255)
                                                        gCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCg                                 >  1:2404275/1‑93 (MQ=255)
                                                         cAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGc                                >  2:1953310/1‑93 (MQ=255)
                                                            aGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCaa                             >  1:674572/1‑93 (MQ=255)
                                                            aGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCaa                             >  1:4921674/1‑93 (MQ=255)
                                                               tctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaa                          >  1:5376763/1‑93 (MQ=255)
                                                               tctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaa                          >  1:4295303/1‑93 (MQ=255)
                                                                ctGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaaa                         >  2:6702202/1‑93 (MQ=255)
                                                                  gAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAt                       >  2:4351265/1‑93 (MQ=255)
                                                                  gAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAt                       <  1:2073753/93‑1 (MQ=255)
                                                                   aaaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAtt                      <  1:4640394/93‑1 (MQ=255)
                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     <  2:4256203/93‑1 (MQ=255)
                                                                     agagTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGCCGAGACCTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTg                    <  1:5714713/93‑1 (MQ=255)
                                                                      gagTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTgg                   >  2:1665054/1‑93 (MQ=255)
                                                                        gTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCt                 >  1:2992101/1‑93 (MQ=255)
                                                                           tCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtc              <  1:1390114/93‑1 (MQ=255)
                                                                            cGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtct             >  2:3139471/1‑93 (MQ=255)
                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            <  1:5209095/93‑1 (MQ=255)
                                                                              tGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTcc           <  2:53576/93‑1 (MQ=255)
                                                                              tGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTcc           <  2:2891733/93‑1 (MQ=255)
                                                                                aCAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGc         >  2:5509998/1‑93 (MQ=255)
                                                                                aCAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGc         >  1:6425817/1‑93 (MQ=255)
                                                                                  aGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCaa       <  2:1379096/93‑1 (MQ=255)
                                                                                   gCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAg      <  2:5853154/93‑1 (MQ=255)
                                                                                   gCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAg      >  1:3593035/1‑93 (MQ=255)
                                                                                    cGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAgg     >  2:1915874/1‑93 (MQ=255)
                                                                                    cGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAgg     >  1:944974/1‑93 (MQ=255)
                                                                                     gggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGc    <  1:5752830/93‑1 (MQ=255)
                                                                                     gggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGc    >  2:1604344/1‑93 (MQ=255)
                                                                                       gAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGt  <  1:4864403/93‑1 (MQ=255)
                                                                                       gAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGt  <  1:395663/93‑1 (MQ=255)
                                                                                       gAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGt  <  2:6164744/93‑1 (MQ=255)
                                                                                       gAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGt  >  1:3370244/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: