Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,684,604 G→A 100% Q339* (CAG→TAG)  yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,6040GA100.0% 595.1 / NA 169Q339* (CAG→TAG) yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (89/80);  total (89/80)

Only 100 of 169 total aligned reads displayed.

AGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACA  >  NC_000913/3684507‑3684699
                                                                                                 |                                                                                               
agaTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAc                                                                                               <  1:2285191/100‑1 (MQ=255)
 gaTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAAcc                                                                                              <  2:624535/100‑1 (MQ=255)
  aTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCt                                                                                             >  2:784882/1‑100 (MQ=255)
   tCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTa                                                                                            <  2:2351004/100‑1 (MQ=255)
     gCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAg                                                                                          >  1:379483/1‑100 (MQ=255)
     gCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAg                                                                                          >  2:2427000/1‑100 (MQ=255)
       cGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGt                                                                                        <  2:3533309/100‑1 (MQ=255)
        gTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTg                                                                                       >  2:1887142/1‑100 (MQ=255)
           gAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGAtt                                                                                    >  1:668386/1‑100 (MQ=255)
           gAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGAtt                                                                                    >  1:1437664/1‑100 (MQ=255)
           gAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGAtt                                                                                    >  2:1933134/1‑100 (MQ=255)
           gAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGATCGCTCTTGATAGTGAGCACTTGCTAACCTAAGGTGAtt                                                                                    >  1:2709023/1‑100 (MQ=255)
            aaCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTg                                                                                   <  2:233900/100‑1 (MQ=255)
             aCATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGc                                                                                  <  2:1343663/100‑1 (MQ=255)
              cATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCg                                                                                 >  2:3368913/1‑100 (MQ=255)
              cATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGATCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCg                                                                                 >  1:2198191/1‑100 (MQ=255)
              cATCACCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCg                                                                                 >  1:3177909/1‑100 (MQ=255)
                tCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTg                                                                               >  2:2169414/1‑100 (MQ=255)
                tCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTg                                                                               <  1:9692/100‑1 (MQ=255)
                 cgccacgccaATGGTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGc                                                                              >  2:2614551/1‑100 (MQ=255)
                 cgccacgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGc                                                                              <  2:1500089/100‑1 (MQ=255)
                 cgccacgccaATGATACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGc                                                                              >  2:2614112/1‑100 (MQ=255)
                     acgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACg                                                                          <  1:3425968/100‑1 (MQ=255)
                      cgccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACgg                                                                         >  2:2891137/1‑100 (MQ=255)
                       gccaATGCTACAGTGCGGACGGAGTTGAATACGTTGAGTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGt                                                                        <  1:1605984/100‑1 (MQ=255)
                       gccaATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGt                                                                        <  1:678380/100‑1 (MQ=255)
                         caATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTc                                                                      <  2:1389164/100‑1 (MQ=255)
                            tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   <  2:2095145/100‑1 (MQ=255)
                            tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   <  1:2606748/100‑1 (MQ=255)
                            tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   >  2:2043063/1‑100 (MQ=255)
                            tGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTg                                                                   >  1:2801499/1‑100 (MQ=255)
                             gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  >  1:2566576/1‑100 (MQ=255)
                             gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  <  1:2517223/100‑1 (MQ=255)
                             gCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGt                                                                  <  2:663721/100‑1 (MQ=255)
                              cTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTa                                                                 >  1:2813163/1‑100 (MQ=255)
                               tACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTac                                                                <  2:2893049/100‑1 (MQ=255)
                               tACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTac                                                                >  2:2699129/1‑100 (MQ=255)
                                aCAGTTCGGACGGAGTGAAATTCGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGGGAGTGCGTGCCACGGTTCCTGTa                                                                 <  1:2713926/98‑1 (MQ=255)
                                aCAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTaca                                                               >  2:325210/1‑100 (MQ=255)
                                aCAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTaca                                                               >  1:98454/1‑100 (MQ=255)
                                     gCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGtt                                                          >  1:2772152/1‑100 (MQ=255)
                                      cggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTg                                                         <  2:2830379/100‑1 (MQ=255)
                                       ggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTgg                                                        <  2:664858/100‑1 (MQ=255)
                                       ggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTgg                                                        <  1:1887129/100‑1 (MQ=255)
                                       ggacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTgg                                                        >  1:3665626/1‑100 (MQ=255)
                                        gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTCAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                       >  1:1265402/1‑100 (MQ=255)
                                        gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                       >  2:1075391/1‑100 (MQ=255)
                                        gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                       <  1:3766433/100‑1 (MQ=255)
                                        gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                       <  1:3569316/100‑1 (MQ=255)
                                        gacggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGc                                                       <  1:3600646/100‑1 (MQ=255)
                                         acggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCa                                                      <  2:1446140/100‑1 (MQ=255)
                                          cggaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCATGGTTCCTGTACACCGTTGGCa                                                      >  1:2569780/1‑99 (MQ=255)
                                          cggaGTTGAATAAGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCaa                                                     >  2:2532328/1‑100 (MQ=255)
                                            gaGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                                   <  2:2926452/100‑1 (MQ=255)
                                              gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGTCAATGAc                                                 <  2:3430799/100‑1 (MQ=255)
                                              gTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGAc                                                 <  2:3704554/100‑1 (MQ=255)
                                               ttGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACa                                                <  2:2593254/100‑1 (MQ=255)
                                                tGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                               >  1:1240911/1‑100 (MQ=255)
                                                 gAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGc                                              <  2:3548739/100‑1 (MQ=255)
                                                 gAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGc                                              >  2:2477806/1‑100 (MQ=255)
                                                   aTACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCaa                                            >  2:1814544/1‑100 (MQ=255)
                                                     aCGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAg                                          <  1:1412520/100‑1 (MQ=255)
                                                       gTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                                        >  2:305437/1‑100 (MQ=255)
                                                          cAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATa                                     >  1:1045393/1‑100 (MQ=255)
                                                           aaTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATaa                                    >  2:1362638/1‑100 (MQ=255)
                                                            aTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                                   >  1:1737450/1‑100 (MQ=255)
                                                            aTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                                   <  2:1436198/100‑1 (MQ=255)
                                                            aTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                                   <  1:3735464/100‑1 (MQ=255)
                                                            aTCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                                   <  2:1524394/100‑1 (MQ=255)
                                                             tCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAcc                                  >  2:3372251/1‑100 (MQ=255)
                                                              cGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCg                                 <  1:533160/100‑1 (MQ=255)
                                                               ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                >  2:317175/1‑100 (MQ=255)
                                                               ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                <  2:3808351/100‑1 (MQ=255)
                                                               ggCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                >  2:3159999/1‑100 (MQ=255)
                                                                gCAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTGCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                                <  1:1126721/99‑1 (MQ=255)
                                                                 cAGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTa                              >  2:3811673/1‑100 (MQ=255)
                                                                  aGGCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTaa                             <  2:2815503/100‑1 (MQ=255)
                                                                   ggCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAt                            >  2:2264863/1‑100 (MQ=255)
                                                                   ggCGCTCGCTCATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAt                            >  2:2124309/1‑100 (MQ=255)
                                                                      gctcgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTg                         >  2:1578537/1‑100 (MQ=255)
                                                                        tcgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcg                       <  1:681953/100‑1 (MQ=255)
                                                                         cgctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcgc                      <  2:530433/100‑1 (MQ=255)
                                                                           ctcATGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGa                    >  2:1603738/1‑100 (MQ=255)
                                                                               tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt                >  1:1432629/1‑100 (MQ=255)
                                                                               tGATAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAAGCTGCGCGGGGa                 >  2:2100305/1‑99 (MQ=255)
                                                                                 aTAGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCa              <  2:2237143/100‑1 (MQ=255)
                                                                                   aGTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATa            >  1:426636/1‑100 (MQ=255)
                                                                                    gTGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACCGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATAc           >  2:2710794/1‑100 (MQ=255)
                                                                                     tGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACg          >  1:34779/1‑100 (MQ=255)
                                                                                     tGAGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACg          >  1:2994107/1‑100 (MQ=255)
                                                                                       aGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTg        <  1:305432/100‑1 (MQ=255)
                                                                                       aGCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTg        >  2:433030/1‑100 (MQ=255)
                                                                                        gCACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTgg       >  1:1740312/1‑100 (MQ=255)
                                                                                         cACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGc      <  2:487268/100‑1 (MQ=255)
                                                                                         cACTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGc      <  1:2132919/100‑1 (MQ=255)
                                                                                          ncTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCg     >  1:1349551/2‑100 (MQ=255)
                                                                                          aCTTGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCg     >  2:1231235/1‑100 (MQ=255)
                                                                                            ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  1:885819/100‑1 (MQ=255)
                                                                                            ttGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGAc   <  2:3300617/100‑1 (MQ=255)
                                                                                             tGCTAACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACa  >  2:851979/1‑100 (MQ=255)
                                                                                                 |                                                                                               
AGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGATCATACGTGGCGACA  >  NC_000913/3684507‑3684699

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