Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,530,797 C→T Q285Q (CAG→CAA sgcX ← putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,7970CT100.0% 126.6 / NA 40Q285Q (CAG→CAAsgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/16);  total (24/16)

GCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGAT  >  NC_000913/4530557‑4531038
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
gCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAAc                                                                                                                                                                                                                                          >  1:191088/1‑250 (MQ=255)
gCCTGAATAAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAAc                                                                                                                                                                                                                                          >  2:103208/1‑250 (MQ=255)
          ggCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGTATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTTTGTTCGAGAGCCGCTTGTTCCAACATCCGTATCAg                                                                                                                                                                                                                               <  1:261371/251‑1 (MQ=255)
                    ggAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCCGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGc                                                                                                                                                                                                                     >  1:25085/1‑251 (MQ=255)
                    ggAAATGTGCTGCTGACATACCTGCCAGAGCGATCAATAGGCGGATGCAATAAGTCAAAACACGCAGGCTGGCAACATACGCCGGAGAAAGGGAATAGCGACAAGGAATAGAGAGAGTGGCGCAGGGAATACAATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCAATTCTCGCTGAACTGGAACGTTGTGTTCAAGAGCCGTTTGTTCCAACCTCAGAATCCGGCGAGGGGGc                                                                                                                                                                                                                     >  1:283186/1‑251 (MQ=255)
                                  gACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGc                                                                                                                                                                                                       <  1:16/251‑1 (MQ=255)
                                     aTACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTAAAGTCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGCGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGATCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGcc                                                                                                                                                                                                      >  1:294045/1‑249 (MQ=255)
                                            ccAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCCGGCTGGCAACATCCGCCGGAGAATGGGTATAGCGACCAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTccc                                                                                                                                                                                             >  1:286811/1‑251 (MQ=255)
                                            ccAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTccc                                                                                                                                                                                             >  1:138891/1‑251 (MQ=255)
                                                             cGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATTAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTc                                                                                                                                                                              >  1:216405/1‑249 (MQ=255)
                                                             cGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTc                                                                                                                                                                              >  1:237617/1‑249 (MQ=255)
                                                              agattcAATCAGTCGAATCACACAGTCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGGCTACCATCCTGCTCCATCTGGATATAACCGGTTTCTGTTATCACACCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGAGAGGAGGCGTGATTCATCCGGACAGCGGTCCCCGACCATGATAGTTcagg                                                                                                                                                                           >  2:133743/6‑251 (MQ=255)
                                                                  gCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTcaggcagg                                                                                                                                                                       >  2:265160/1‑251 (MQ=255)
                                                                         gTCAAATCACGCAGGCTGGCACCCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATccc                                                                                                                                                                <  2:118943/251‑1 (MQ=255)
                                                                           cAAATGACGCGGACTGGAAACCTCCGCCGGAGAATGTGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCCACCTGGATATACCCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCaa                                                                                                                                                              <  2:240628/251‑1 (MQ=255)
                                                                             aaTCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACg                                                                                                                                                            >  2:253696/1‑251 (MQ=255)
                                                                              aTCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGtt                                                                                                                                                                               <  2:72505/231‑1 (MQ=255)
                                                                              aTCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGtt                                                                                                                                                                               >  1:72505/1‑231 (MQ=255)
                                                                                             aaCCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTg                                                                                                                                             >  1:279417/1‑250 (MQ=255)
                                                                                               ccTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGAcc                                                                                                                                          >  1:173586/1‑251 (MQ=255)
                                                                                                cTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCt                                                                                                                                         >  1:262008/1‑251 (MQ=255)
                                                                                                           aaTGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCg                                                                                                                              >  2:217763/1‑251 (MQ=255)
                                                                                                           aaTGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCg                                                                                                                              <  1:265160/251‑1 (MQ=255)
                                                                                                              gggTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCgtg                                                                                                                            <  2:57524/250‑1 (MQ=255)
                                                                                                                  ataGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAAt                                                                                                                       >  1:21302/1‑251 (MQ=255)
                                                                                                                    aGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCGGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGAGTAATCCTGACCTCGGAATAATCGTGTAAATc                                                                                                                      <  1:14552/250‑1 (MQ=255)
                                                                                                                      cGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAgg                                                                                                                   <  1:131707/251‑1 (MQ=255)
                                                                                                                                                        cATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATc                                                                                 <  1:204658/251‑1 (MQ=255)
                                                                                                                                                                  aaCCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCcaggtcag                                                                       <  2:138891/251‑1 (MQ=255)
                                                                                                                                                                    ccTGGTTATAACCGGTTTCTGTTATCACGTCGGGAGCCTCTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTTATATCAATACCAATCGCcaggtcaggt                                                                     <  1:25474/251‑1 (MQ=255)
                                                                                                                                                                      tGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCAGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCACCGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATAAATACCAATCGCCAGGTCAGGttt                                                                   >  2:200590/1‑251 (MQ=255)
                                                                                                                                                                         atatAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTac                                                                >  2:76802/1‑251 (MQ=255)
                                                                                                                                                                                          tATCACGCCGGGAGCCACTTCTCGCTGCACCGGCATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAAt                                                                                                                                                 >  2:303416/1‑153 (MQ=255)
                                                                                                                                                                                          tATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGTCGAGGTGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGCAATCAAGACCAATCGCCAGGGCAGGTTTTACACGGCGGAATACGGGa                                                >  2:146450/1‑250 (MQ=255)
                                                                                                                                                                                          tATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTTATTAATCCGGCCAGCGTGCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAAt                                                                                                                                                 <  1:303416/153‑1 (MQ=255)
                                                                                                                                                                                             cACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACa                                             <  2:242437/250‑1 (MQ=255)
                                                                                                                                                                                             cACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACa                                             >  2:67699/1‑250 (MQ=255)
                                                                                                                                                                                                                              gttgtGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAGTCTTCCTGTACCGAAg           <  1:67699/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gtTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGccacc      <  2:106982/251‑1 (MQ=255)
                                                                                                                                                                                                                                        agagCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCATGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGAt  >  1:171935/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
GCCTGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGAT  >  NC_000913/4530557‑4531038

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